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Механические и конечноэлементные модели кератопротезов 
роговицы глаза 


А. Н. Соловьев, Н. И. Глушко", А.Н. Епихин”, М. Свейн`, О. Н. Лесняк,, А. Е. Иванов' 
'ФГБОУ ВО «Донской государственный технический университет» (г. Ростов-на-Дону, Российская Федерация) 
`ФГБОУ ВО «Ростовский государственный медицинский университет» (г. Ростов-на-Дону, Российская Федерация) 
"Сиднейский университет (г. Сидней, Австралия) 





Введение. При разработке глазных протезов возникает ряд проблем, одна из которых — конструирование со- 
единения между жесткой оптической частью и мягкой тканью роговицы. Их модули Юнга могут отличаться на 
три порядка. При этом возникает задача создания промежуточного слоя, возможно с градиентными свойствами, 
назначение которого — исключить травмирование мягких биологических тканей. Рассматривается два типа 
кератопротеза: с опорной пластиной и с промежуточным функционально-градиентным слоем. Для первого типа 
рассчитано напряженно-деформированное состояние. Для второго проведено аналитическое и конечноэлемент- 
ное моделирование взаимодействия цилиндрического оптического протеза, промежуточного неоднородного 
слоя и роговицы в упругой среде. Рассмотрены два варианта: без учета кривизны (круглая плита или пластина) 
и с учетом кривизны (сферический купол или оболочка). Цель работы — исследование напряженно- 
деформированного состояния керапротеза и роговицы в области контакта. 

Материалы и методы. Математические модели рассматриваемых конструкций — это краевые задачи линейной 
теории упругости. Аналитическое решение построено для упрощенной модели в виде составной круговой пла- 
стины. Пространственные трехмерные задачи и задачи в осесимметричной постановке решаются методом ко- 
нечных элементов. Конечноэлементное моделирование рассматриваемых конструкций проводилось в САЕ- 
пакете АМЗУЗ и АСЕГАМ. 

Результаты исследования. Построены САО-модели кератопротезов с условиями закрепления и нагружения. 
Установлена нагрузка на керапротез под действием внутреннего глазного давления. Рассчитано напряженно- 
деформированное состояние элементов кератопротеза и роговицы. Особое внимание уделено окрестности ее 
контакта с кератопротезом. 

Обсуждение и заключения. Результаты расчета осевых смещений и механических напряжений в кератопротезе 
первого типа показывают, что выбранные для него геометрические параметры удовлетворяют кинематическим 
и прочностным требованиям. Предложенные в работе модели деформированного состояния мягких биологиче- 
ских тканей позволяют оценить их травмирование при использовании кератопротеза второго типа, а также вы- 
брать геометрические параметры и градиентные свойства промежуточного слоя. 


Ключевые слова: глазной протез, неоднородные упругие свойства, пластина, оболочка, метод конечных эле- 
ментов. 
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Введение. Кератопротез — это бесклеточный искусственный имплантат, сконструированный таким 
образом, что цилиндрическая оправа удерживает оптику. Кератопротез заменяет удаленную часть роговицы. 
Раньше трансплантации роговицы имели высокие показатели инфицирования и отторжения. В конце 1980-х 
годов наибольшее распространение получила конструкция «ядро и юбка» (биоинтегрируемая «юбка» окружает 
оптику). Важны были не только размеры, но и расположение пор в пористой юбке. Такие кератопротезы, как 
А1рваСог (раннее название Сви!а) представляли собой полиметилметакрилатное устройство с центральной 
оптической областью, слитой с окружающей губчатой юбкой. 

Современные кератопротезы состоят из оптического элемента и опорной пластины. Оптический про- 
зрачный элемент имеет форму цилиндра или линзы. Опорная пластина, соединенная с оптическим элементом, 
может иметь различные формы: колечка с отверстиями, колесика со спицами, ушек или амебовидных ножек. В 
книге С. Н. Федорова [1] приведены описания различных видов и форм кератопротезов и креплений. Основное 
осложнение после кератопротезирования — асептический некроз роговицы, развивающийся спереди от опор- 
ного крепления имплантата. Именно это осложнение довольно часто приводит к отторжению кератопротеза [1 |. 
Основная причина асептического некроза — блокирование поступления жизненно необходимых веществ в 
слои роговицы, которые располагаются над опорой кератопротеза. В связи с этим для изготовления опорной 
пластины необходимо использовать биосовместимый материал, который позволит тканям роговицы прорасти 
сквозь опорное крепление. 

Американские исследователи описали микропористую опорную пластину" для кератопротеза из растя- 
нутого политетрафторэтилена (ПТФЭ). Его структура имеет конфигурацию в виде узлов полимера, которые 


' Патент США № 5713956, М. кл. 7 А 61 Е 2/14, опубл. 03.02.1998. 
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соединяются фибриллами длиной 7-8 мкм. Слишком маленькие поры препятствуют прорастанию тканей рого- 
вицы. Поэтому необходимо дополнительно пронизывать протез порами, перпендикулярными поверхностям 
опорной пластины диаметром 20-150 мкм (лучше 50 мкм). Толщина опорной пластины должна составлять око- 
ло 0,2 мкм, но не более 0,3 мкм. 

В модели с опорным креплением в форме кольца" оптический элемент сделан из прозрачного вещества 
типа полиметилметакрилата (ПММА), крепление — из гидрофильного пористого материала с волокнистой 
структурой, через который могут прорастать ткани роговицы. Из этого волокна сформирована опорная пласти- 
на толщиной 0,15-—0,30 мкм. Кольцо с внешним диаметром около 9,5 мм надето на оптический элемент. Его 
передняя часть — цилиндрическая, задняя имеет форму усеченного конуса, большее сечение которого направ- 
лено внутрь. Опорная пластина имеет упругие характеристики, необходимые для предотвращения отторжения 
кератопротеза и некроза из-за давления окружающих тканей глаза. 

В другой конструкции” опорная пластина представляет собой две параллельные оболочки, отстоящие 
друг от друга на 0,4-0,7 мкм. Оптимальное расстояние 0,3=0,02 мм. Оболочки конической формы прикреплены 
к цилиндрическому оптическому элементу. Материал опорных пластин должен быть биозаселяемым, следова- 
тельно, предполагается, что его пористость — 50% и выше, поры должны быть открыты, их диаметр — около 
20—100 мкм. 

Элементы кератопротеза”, предложенного российскими учеными (рис. 1): 

— опорная пластина — 1 (радиус кривизны 7-10 мм, внешний диаметр 7-12 мм); 

— оптическая часть — 2 (оптический элемент в форме цилиндра с закругленными торцами, которые служат 
для выбора диоптрий); 

— цилиндрическая поверхность оптического элемента с винтовой нарезкой — 3; 

— соединительная шайба — 4 (жестко скреплена с опорной пластиной). 





Рис. 1. Схема кератопротеза 


Кератопротез имплантируется в интраламеллярный карман роговицы. Опорная пластина располагается 
внутри роговицы глаза, а оптический элемент пронизывает всю роговицу. 

Имплантат для укрепления роговицы” представляет собой пластину — круглую, овальную, в форме 
трилистника, ромашки или выпукло-вогнутой линзы с радиусом кривизны 7-10 мм. В центре линзы может 
быть цилиндрическое отверстие. В структуре материала имплантата объемная доля пустот — 15—40 %, их 
удельная поверхность — 0,25-—0,55 мкм”/мкм°, среднее расстояние между пустотами — 25—50 мкм, средняя 
объемная хорда — 8—25 мкм. Имплантат представляет собой выпукло-вогнутую линзу, радиус кривизны кото- 
рой соответствует радиусу кривизны роговицы пациента. В центральной части толщина имплантата — 0,3— 
0,7 мм, на краях — 0,01 мм. 

В другой модификации кератопротеза”, в отличие от описанного выше, опорный блок выполнен из по- 
ристого никелида титана, кривизна пластины совпадает с кривизной роговицы. Опорная пластина устанавлива- 
ется на поверхности роговицы и ретенционно фиксируется аллотрансплантатом склеры. Цилиндрический опти- 
ческий элемент соединяется с опорной пластиной методом тугой посадки. Такой подход позволяет более эф- 
фективно лечить бельма, благодаря быстрой и успешной интеграции пористого имплантата (базиса) с окружа- 


' Патент США № 5489301, М. кл. 6 А 61 Е 2/14, опубл. 05.02.1996. 

` Патент США № 6106552, М. кл7 А 61 Е 2/14, опубл. 22.08.2000. 

° Кератопротез. Пат. 2270643 Рос. Федерация. № 2004119480/04; заявл. 17.06.2004; опубл. 27.02.2006. Бюл. № 6. 8 с. 

*Имплантат для укрепления роговицы. Пат. 2270642 Рос. Федерация. № 2004119479/04; заявл. 17.06.2004: опубл. 27.02.2006. Бюл. № 6. 9 с. 
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ющими тканями и механически прочной фиксации оптического элемента, которая снижает опасность репози- 
ции кератопротеза. 

Разработчик комбинированного кератопротеза' отмечает, что опорная часть имплантата должна мак- 
симально сопротивляться выталкивающему воздействию давления внутриглазной жидкости на оптическую 
часть кератопротеза и при этом минимально деформировать надлежащие ткани. Материал опорной части дол- 
жен быть не только биосовместим, но и структурирован таким образом, чтобы обеспечить сращение над- и 
подлежащих слоев роговицы, разделяющихся в процессе имплантации кератопротеза. Ранее упоминались кера- 
топротезы с оптическим элементом из прозрачного вещества (например, из ПИММА) и опорной пластиной в 
форме кольца из гидрофильного пористого материала ПТФЭ или полиэтилена”. Их использование обеспечива- 
ет сращения над- и подопорных слоев роговицы за счет высокой пористости материала (от 50 %) и диаметра 
пор до 100 мкм. Главные недостатки таких конструкций — избыточная гибкость и низкая механическая проч- 
ность. Эти характеристики не позволяют продолжительное время надежно удерживать и фиксировать установ- 
ленный кератопротез с большой оптикой. 

Оптическая часть комбинированного кератопротеза выполнена из прозрачного упругого полимера, 
насыщенного УФ-адсорбентом. По форме это съемный болт диаметром 5—6 мм. Передний торец — грибовид- 
ный, другие — сферические (асферические). Гаптическая часть — это опорный диск наружным диаметром 8— 
12 мм и толщиной 0,3-0,9 мм. На ее боковую поверхность нанесена резьба. Для этой части используют пори- 
стый (перфорированный) твердый (упругий) ареактивный полимер (металл) с объемной долей пустот не более 
50 %. Переднюю часть гаптического элемента делают из биологических материалов, заднюю — из искусствен- 
ных биосовместимых. Опорный диск задней части — плоский и перпендикулярно соединяется с размещенным 
в центре опорным полым цилиндром. Соединение обеспечивает канавка и (или) фланец на внешней задней ча- 
сти цилиндра. Профиль и диаметр внутренней резьбы цилиндра соответствуют резьбе оптического болта. 
Опорный диск может переходить во втулку, которая обхватывает опорный цилиндр снаружи. Опорная пласти- 
на спереди соединена с круглым биодиском, который вырезается из свежеизъятых или консервированных сте- 
рильных биологических материалов, которые представляют собой разновариантные послойные сочетания био- 
логических тканей. Перед соединением части кератопротеза стерилизуют. Соединение проводят в стерильных 
условиях. Сборку выполняют непосредственно перед протезированием или заранее (в этом случае имплантат 
хранится до операции). Внешний диаметр биодиска должен быть больше внешнего диаметра опорного диска на 
0,3—1 мм, в центре биодиска делается соосное отверстие с диаметром, соответствующим наружному диаметру 
покровной втулки или опорного цилиндра. 

Для сквозной кератопластики разработаны модели с использованием различных материалов, конструк- 
ций и хирургических методов [2]. Однако проблема окончательно не решена. Это подтверждают случаи слепо- 
ты роговицы у пациентов с повторной недостаточностью трансплантата или с тяжелым поражением глазной 
поверхности [3—6]. Вопросы кератопротезирования при заболеваниях роговицы и глазной поверхности рас- 
смотрены в [7]. Обзор по применению бостонского кератопротеза представлен в [8]. В [9] обсуждаются ослож- 
нения с широкими последствиями при применении кератопротезов 

В настоящей работе рассматриваются два типа кератопротезов с оптическим элементом цилиндриче- 
ской формы. Схема первого типа представлена на рис. 1. Во втором типе на цилиндрической поверхности при- 
сутствует некоторый промежуточный слой с неоднородными механическими свойствами. На рис. 2 представ- 
лена схема половины осевого сечения кератопротеза второго типа без учета кривизны (а) и с учетом кривизны 
роговицы (6). Серым выделен слой с функционально градиентными свойствами. Он нужен для нетравматиче- 
ского контакта оптики с мягкими тканями роговицы. 





а) 6) 


Рис. 2. Половина осевого сечения кератопротеза: модель пластины (а), модель купола (6) 


' Кератопротез комбинированный. Пат. 2707646 Рос. Федерация. № 2018146759; заявл. 27.12.2018; опубл. 28.11.2019. Бюл. № 34. 5 с. 
* Патенты США № 5489301, № 6106552. 
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Методом конечных элементов построены аналитические и численные модели имплантации кератопро- 
теза в роговицу. Исследовано напряженно-деформированное состояние роговицы в области контакта. Этим 
определяются цели настоящей работы. Для керапротеза первого типа: исследование его напряженно- 
деформированного состояния под действием внутреннего глазного давления в предположении, что внешний 
край опорной пластины закреплен. Для керапротеза второго типа: исследование напряженно- 
деформированного состояния роговицы в окрестности контакта с кератопротезом, внешний слой которого 
представляет собой конструкцию, позволяющую уменьшить травмирование мягких тканей. Этот промежуточ- 
ный (интерфейсный) слой моделируется как полый цилиндр с функционально градиентными механическими 
свойствами. 

В работе рассматриваются три модели. 

Г. — Аналитическая — для имплантированного протеза на основе изгиба составной круговой пластины. 

П. ’Конечноэлементная — для имплантированного протеза на основе составной круговой плиты. 

Ш. Конечноэлементная модель — для имплантированного протеза на основе составного сферического 
купола. 

Решаются задачи моделирования для: 

— кератопротеза первого типа (А), 

— кератопротеза второго типа (Б). 


Материалы и методы. Континуальная постановка задач. Общая математическая постановка исследу- 
емой проблемы (для задач А и Б) представляет собой статическую краевую задачу теории упругости для со- 
ставного изотропного тела. Для неизвестных и = (и, и, и.) — компоненты вектора перемещений система диф- 


ференциальных уравнений имеет вид [10]: 


о,, =0, 54 бы +щи). (1) 
где б,, Си, &и — компоненты тензоров напряжений, упругих постоянных и деформаций соответственно. 
Граничные условия ставятся для векторов смещений и напряжений на соответствующих поверхностях 
и: 
и, |5 =и, (х,0), ХЕ5,, (2) 
Г, |5 =0,И, 5, = 9(%,1), хе», , (3) 


где п, —_ Координаты единичного вектора внешней нормали. 


Кроме того, для задачи Б (рис. 2) подтела имеют различные свойства, а именно: 

— два из них (цилиндрический оптический протез и роговица) — однородные, с модулями упругости 
Еги Ек,; 

— функционально-градиентный модуль упругости интерфейсного слоя — Ех, = Ею (Г). 


Правая сторона закреплена, на левой задаются условия симметрии. На нижнюю границу действует рав- 
номерно распределенное давление, которое соответствует избыточному внутриглазному давлению по сравне- 
нию с атмосферным. На интерфейсных границах задаются условия непрерывности. 

Для исследования модели (Г) краевая задача (1)—(3) сводится к системе обыкновенных дифференциаль- 
ных уравнений (4), (5) и (6) для первого, второго и третьего участка соответственно относительно угла поворо- 
та нормали 6(’) при степенной зависимости от радиуса модуля упругости второго участка [11]: 


00) г? + (= в») г- 0(2) = -Ктз, (4) 
“0 (1) т? +а (= в(»)) т+а9>0(г) + (= в(»)) т — 9(`) = НВ (5) 
45007) |"? + (40065) " — 065) = +кза("2 — Фу, (6) 


где К =1/), р — цилиндрическая жесткость. 
Результаты исследования 
Аналитическое решение 
Общее решение системы (4)—(6) представлено ниже. 
Первый участок: 
0(7) =хС, + Кат? + — 
Третий участок: 
9() = - + тб. — --КЗа(4 ш(и)т2 — 72). (7) 
Второй участок: 


11 11 
—5а+=/а*-—4а92+4 —=а+=/а*-—4а92+4 
ыы 
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1К›ат"“ (67 — 72 )9) — т? + З7Р)а — в") 
ы (8+ (@, —3)а)а@.-1 


при 
ЕТ1 
ЕТ2 й. —_ шт (Ее) 


121-92) “= ШС.) 





Ет> —= Ет (2); 2> — 


Произвольные постоянные С, общего решения (7) определяются из граничных условий на правом кон- 
це и условий стыковки, причем для первого участка С, =0. 

Конечноэлементное решение 

Численное моделирование проведено в конечноэлементных пакетах АМЗУЗ и АСЕГАМ [12, 13]. 

Численные результаты для кератопротеза первого типа. На рис. 3 представлена САО-модель кера- 
топротеза с условиями закрепления и нагружения. 





$] 
0.000 9 000 (тити Хх 
Ти 


4,500 
Рис. 3. Геометрическая модель с условиями закрепления и нагружения 


Внешний край опорной пластины закреплен, на оптический элемент действует распределенное давле- 
ние. Нормальное внутриглазное давление составляет от 12 до 21 мм ртутного столба. Давление в 45 мм/рт. ст. 
является критичным. В этом случае операция еще может спасти глаз. При более высоком давлении глаз спасти 
практически невозможно. На данный кератопротез действует нагрузка в 46 мм/рт. ст. 

На рис. 4 представлено распределение осевого смещения (а) и напряжения по Мизесу (6) на конечно- 
элементной сетке модели. 





ты 
а 
Механика 


Рис. 4. Распределение на конечноэлементной модели характеристик напряженно-деформированного состояния: 
осевого смещения (а) и напряжений по Мизесу (6) 


Максимальные напряжения возникают в месте соединения оптического элемента с опорной пластиной, 
но их значения не превосходят предела прочности выбранного материала. Э55 
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Далее рассматриваются результаты расчетов для задачи Б в осесимметричной постановке. Так, на 
рис. 5—7 представлены компоненты напряженно-деформированного состояния (П). 





а) 0) 


Рис. 5. Расчет перемещений: конечноэлементная сетка на деформированной конструкции (а), 
распределение вектора смещений (6) 


О: 


О’. Н/мм" 30 ЕЕ О’. Н/мм" 


25 НИ 


20 





0,1 


т 


10 





А В 


од, Н/мм* 


в) г) 
Рис. 6. Расчет компонентов тензора напряжений. Зависимость от радиуса: радиальных напряжений на верхней границе (а), 
на интерфейсной границе с роговицей (6), угловых напряжений на верхней границе (в), на интерфейсной границе 
с роговицей (г) 
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Рис. 7. Распределение компонент напряжений внутри области: радиальных напряжений (а), угловых напряжений (6) 


Характеристики напряженно-деформированного состояния модели Ш представлены на рис. 8—11. 
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Рис. 8. Расчет перемещений купола: конечноэлементная сетка на деформированной конструкции (а), 


12 


Рис. 9. Расчет радиальных напряжений: распределение в области (а), на интерфейсной границе с роговицей (6) 
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Рис. 10. Расчет угловых напряжений: распределение в области (а), на интерфейсной границе с роговицей (6) 
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Рис. 11. Расчет касательных напряжений: распределение в области (а), на интерфейсной границе с роговицей (6) 


Обсуждение и заключения. Рассмотрены два типа кератопротезов и построены их механические и ма- 
тематические модели. Для первого типа кератопротеза с опорной пластиной исследовано его напряженно- 
деформированное состояние, создана конечноэлементная модель в АМЗУЪЗ. Показано, что при максимальном 
глазном давлении возникающие максимальные напряжения не превосходят пределов прочности выбранных 
материалов. Для второго типа кератопротеза построено аналитическое решение без учета кривизны и конечно- 
элементное решение в пакете АСЕГАМ с учетом кривизны роговицы. Рассчитаны характеристики напряженно- 
деформированного состояния, в том числе на интерфейсной границе с роговицей. Это позволяет оценить ее 
травмирование и выбрать геометрические параметры и градиентные свойства промежуточного слоя. 
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Частотная функция и функция демпфирования в оценке динамических т 
процессов в механических колебательных системах с симметрией | 


А. В. Елисеев у сЕ 


ФГБОУ ВО «Иркутский государственный университет путей сообщения» (г. Иркутск, Российская Федерация) 


Введение. Рассматривается новый подход в формировании методологического базиса системного анализа в 
приложении к задачам динамики механических колебательных структур. Целью исследования является 
разработка метода оценки свойств механических колебательных систем с учетом сил вязкого трения на основе 
частотных функций и функции демпфирования, зависящих от так называемого коэффициента форм связности, 
представляющего собой отношение характеристик обобщенных координат. 

Материалы и методы. Используемые графоаналитические методы оценки динамических свойств 
механических колебательных систем с двумя степенями свободы основаны на определении экстремальных 
значений частотных функций и функции демпфирования, определяемых из соотношений, связывающих 
кинетическую, потенциальную энергию и значения функции рассеяния энергии. Математические модели 
строятся на основе формализма Лагранжа, матричных методов, элементов теории функций комплексной 
переменной. 

Результаты исследования. Предложен метод построения частотных функций и функций демпфирования для 
класса механических колебательных систем с двумя степенями свободы на основе аналитических выражений, 
отражающих особенности соотношения потенциальной и кинетической энергии системы с учетом сил вязкого 
трения, представленных диссипативной функцией. Выведены общие аналитические выражения для частотной 
функции и функции демпфирования. Для механических колебательных систем с упруго-демпфирующими 
элементами, обладающими свойствами симметрии, проведен графоаналитический анализ экстремальных 
свойств соответствующих частотных функций и функций демпфирования. Представлены результаты 
численных экспериментов. Предложен критерий классификации частотных функций и функций 
демпфирования на основе топологических особенностей графиков соответствующих функций. 

Обсуждение и заключения. Разработанный метод построения частотных функций и функций демпфирования 
может быть использован для отображения динамических особенностей механических колебательных систем. 
Предложенный матричный метод построения частотно-демпирующей функции для системы двумя степенями 
свободы может быть распространен на механические колебательные системы, рассматриваемые в различных 
системах координат. 


Ключевые слова: механическая система, динамические связи, частотная функция, функция демпфирования, 
связность движений, экстремальные свойства, колебание, вязкое трение. 
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Введение. Методам использования механических колебательных систем в качестве расчетных схем в 
задачах оценки динамических свойств технических объектов, работающих в условиях интенсивных 
вибрационных нагрузок, уделяется существенное внимание [1-9]. К числу распространенных подходов к 
оценке динамических свойств механических колебательных систем можно отнести методы, основанные на 
экстремальных свойствах отношения потенциальной и кинетической энергии [10, 11]. Методы, построенные на 
основе энергетических отношений, получили развитие в использовании частотной функции, как функции 
коэффициента форм связности координат механической системы, для оценки особенностей динамики 
механических колебательных систем без учета сил трения [12—15]. 

Вместе с тем, методы оценки динамических свойств механических колебательных систем, с учетом сил 
вязкого трения на основе частотной функции, требуют детализации представлений в зависимости от величины 
вязкого трения. Это связано, в частности, с тем, что для систем с апериодическим движением понятие частоты 
колебания может терять смысл. 

Предлагаемая работа посвящена развитию метода оценки свойств механических движений на основе 
развития понятия частотной функции, когда вводится дополнительная функция демпфирования, отражающая 
особенности учета сил вязкого трения в зависимости от коэффициентов форм связности. 

Материалы и методы. Рассматриваются свободные движения механической упруго-диссипативной 
системы с сосредоточенными параметрами с двумя степенями свободы. Принципиальная схема системы 
представлена на рис. 1. 


У! У? 


[2 Бо Ь? 


ы оС оС А 


Рис. 1. Механическая колебательная система с учетом вязкого трения 
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Обобщенные координаты у, у» обозначают смещение массо-инерционных элементов т, т, 


относительно положения статического равновесия. Кинетическая энергия ТГ, потенциальная энергия П и 
функция рассеяния Ё имеют вид: 


Г. | В 
Т=оту +5тьу, , (1) 
Ре И 2 
В. т м. > (2) 
1 В 1 : . 2 1 ‚2 
г И +5%0,-Я) о | (3) 
Система уравнений Лагранжа второго рода имеет вид: 
Ч ОТ ОТ ОоЕ 
+ =0; 
и, д, 9 7 
Я ОТ ОШ ОЕ 
ати 
@ 0, Оу, 0, 
После подстановки выражений Т,П, ЕЁ система дифференциальных уравнений (4) принимает вид: 
ы +(5 +5), Б.у, + (№ +А) у, — Ко, = 0; (5) 
т, у, + (6 +5, )у, -Б У + (& + )у, —К у, =0. 


Формы свободных движений представленной системы (5) в общем случае определяются собственными 
числами характеристического уравнения с учетом их кратности. Рассматривается случай простых корней. 
Таким образом, пусть решение у = у (Ру =у,() системы (5) представимо в виде: 


р”, (6) 


а - |И . . 
где у - — вектор-решение, У = —щ числовой вектор, р=о+ 76 — комплексный параметр, Е — 
2 2 


переменная времени. Полагается, что начальные условия согласованы с видом искомого решения (6). 

Задача заключается в построении и оценке экстремальных свойств функций, отображающих 
характеристики собственных движений системы с учетом сил вязкого трения. 

Результаты исследований 


1. Построение частотной функции и функции диссипации на основе энергетического 
соотношения. Система (5) в обозначениях (6) имеет вид: 


бы -№Р-№ | (1) 
-6Р-№ т, р” + (Б, +Ь, р + К, || № 
Введем обозначения: 
= а —Б, с" —№ | (8 
О т, —/ ВБ -Ь, —«ю Ю№-К, 
С учетом (8) матричное соотношение (5) принимает вид: 
(р?А+рВ+С)У =0. (9) 


На основе матричного соотношения (9) могут быть получены различные скалярные уравнения, с 
учетом экстремальных свойств которых, в свою очередь, могут быть определены свойства решения р уравнения 


(9). В частности, скалярное умножение левой и правой части равенства (9) на вектор У приводит к скалярному 


выражению: 
р? (А7,У)+р(В7,У)+(СТ,У)=0. (10) 
С помощью подстановки р=о+ 0 приведем (10) к виду: 
(0*—@7 +2/06)(А7,7)+(в+ /)(ВУ,У)+(СУ,У)=0, (11) 
Пусть для координат вектора У выполнено соотношение: 
У, =а7, (12) 
где © — коэффициент формы связности. В этом случае вектор У может быть представлен в виде: 
та. (13) 
где @ = В — вектор, определяемый коэффициентом связности @. После подстановки (13) выражение (11) 


принимает вид: 
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(9—0 +2/о)А, + (0+ /)В, +С, =0, (14) 
где А = (Аа, 6) ‚ В = (Ва, 6) ‚В. = (Ва, 6) — скалярные функции коэффициента формы связности &@. 
Уравнение (14) может быть представлено в эквивалентном виде: 

<’ А, =о”А, +оВ, +С.; 
20%А, +®В, =0. 


Найдём решение системы (15) относительно ©,©, как функций ©. Функция (0) рассматривается как 


(15) 


частотная функция, о(9) — функция демпфирования. Особенности решения системы определяются знаком 
2 
выражения о’4А +оВБ, +С.. 
Учет диссипации. Уровень диссипации энергии может быть охарактеризован дискриминантом 


В’ -4А,С,. При условии: 














В” <4АС, (16) 
понимаемом как малость трения, решение (15) может быть представлено в виде: 
2 С. В. 
А, 24)’ 
{0 {0 (17) 
В. 
о=- 
24. 
Необходимо отметить, что при выполнении условий малости сил вязкого трения (16) выполнено равенство: 
С 
2 2 
ю +о =—“. 18 
7 (18) 
При условии «больших сил вязкого трения»: 
В’ >4А.С., (19) 
решение (15) может быть представлено в виде: 
«=0, (20) в (@) 5. (21) 
— : о = — — 
| 2А, 
Функция демпфирования имеет две компоненты с’ (0), ©,(@,) такие, что: 
В С 
и, > (22) 


Таким образом, в зависимости от уровня сил вязкого трения, возможны различные виды представления 


частотной функции и функции демпфирования. Если В ^-4А.С, <0, то определена частотная функция > (0) 
и одна компонента функции демпфирования ©(0). Если В^-4А.С. >0, то предполагается, что частотная 
функция ©’(9) принимает нулевые значения, а функция демпфирования имеет две различные отрицательные 
компоненты с’ (©), с, (0). Что касается условия В *-4А.С, =0, то оно требует отдельного анализа. Вместе с 


2 
тем, условие В ^-4ААС =0 может быть интерпретировано как граница двух различных режимов движения 
[02 [02 [02 


механической системы. 

Представленные аналитические выражения частотной функции, функции демпфирования и условий 
«малых» и «больших» сил вязкого трения могут быть детализированы при рассмотрении частных вариантов 
механических колебательных систем, полученных на основе системы с двумя степенями свободы. 

2. Частотная функция и функция демпфирования для механической системы с двумя степенями 
свободы. Рассматриваются варианты параметров для механической системы, представленной на рис. 1. 
Предполагается, что множество граничных параметров, разделяющих режимы движение при малых и больших 
силах вязкого трения, определяется уравнением: 


В.’ =4А.С., (23) 
где: 
А =т +т, о”, (24) 
В .=(Ь,+Ь, ог —-2а6, +В. +В, (25) 
С =(® +К,)а” =? ++. (26) 


Условия малости сил вязкого трения имеют форму: 
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В.’ <4АС.. (27) 


В этом случае частотная функция <’ (0) и функция демпфирования о(0): 


2 
(а) = (®-А,) о — 29% ++ [1 (+, — 296, +В, +В, 
т, и 2 т, +т, а” (28) 
в 
ое Вы 
2 т, + т. 


Условия больших сил вязкого трения имеют форму: 
В’ >4АС.. (29) 


При условиях (29) функции ©’ и о(а) имеют вид: 


_ 1 (6+6, )7 —-2а6 +В. +В _ 


2 т, +т, а” 


О (0) — 








— [16 -+Ь)а” —296, +6, +6, (+6)? -20% +6 +8. 0) 
2 т, +т, о т, +т, о 
я 
_2 
И ав 
2 т, + та 


1 (6 -Ь, а —296, +В, +В, 


й 
Е т, + т, о” | 


Представленные выражения отражают движения в форме экспоненциального убывания при отсутствии 





_ (+) — 20, + + 


2 
т, + т.о 





колебаний. 
3. Особенности частотных функций и функций демпфировании для симметричных механических 
колебательных систем. Рассматривается механическая колебательная система с упруго-демпфирующими 


элементами, на значения параметров которых наложены условия симметрии в виде ВБ =Б, =ВБ =Б, 


К =А, =А =К. Принципиальная схема представлена на рис. 2. 


У? У? 
Ь Ь Ь 
т тТ>? 
у: К 
ОО О: © 
\\\ \\\ 


Рис. 2. «Симметричная» механическая система 


Система дифференциальных уравнений (5) имеет вид: 


ы +2 — 9, +25 №, =0; ых 
ту, +2Бу, -Бу +2Ку, -Ку =0. 
Функции А., В., С, могут быть представлены выражениями: 

А =т +т, о”, (32) 

В .=2Ы о’ -@+1), (33) 

С =2К(а/ -@-+1). (34) 


На основе представленных компонент могут быть построены частотная функция, функция 
демпфирования и сформулированы условия малости сил вязкого трения. 
Учет сил вязкого трения. Условие малости сил трения могут быть представлены неравенством: 


Во С 
В | <—. (35) 





[0 
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После подстановки функций (32)—(34) условие малости сил трения (35) может быть представлено в 
виде: 
у, <М(а), (36) 
.. й [т + т, а.” .. .. 
где у, = — — обобщенный вязко-упругий параметр, М (а) =--————— — обобщенный массо-инерционный 
4К 2 а -@а+1 
коэффициент, зависящий от коэффициента формы @&. График функции М. для каждого фиксированного 7, 
определяет множество значений ©, при которых выполнено условие малости сил трения. 


В качестве примера на рис. 3 приведен график параметризующей функции М(а). Функция М(@) 


2 


т 
обладает глобальным минимумом М, и максимумом М,, при а —00 М(@)-— Е 


Ус 
(1) 








— 


— — = — — —— -_ — ——_—_—_— =—— 


-10 -5 0. 5 10 





{0 


Рис. 3. Параметризующая функция М (0): 1 — уровень глобального максимума М, , 3 — уровень глобального минимума 


т 
М, , 2 — уровень горизонтальной асимптоты — 


Изменение параметра ‘у, в интервалах (,М.) (М.Р), С,М,) (М9) определяет характерные 


интервалы коэффициента формы ©, в которых выполнены условия малости сил трения. 
При малых силах трения частотная функция и функция демпфирования имеют вид: 


— 2к(о” —9+1 _ Ба” —@-+1) и 


о (©) т, + та т +т.а? |’ 
1 2 1 2 
(37) 
(о? — а +1 
= 
т, + т, а, 


При больших силах трения, при которых реализуется апериодическое движение системы, частотная 
функция равна нулю, а функция демпфирования имеет две компоненты: 








Ф’ =0; 
(о? — а +1 (а? о -+1 | 2(о? —а-1 
ыы. (88) 
т, + т. а т, + т. т, + т. 
(а? — “+1 а, ЭЕы а 
(а-я, |5 2% ) 





й 2 
т, + т.“ т, + т, о, 


На основе аналитических представлений частотной функции и функции демпфирования могут быть 
определены характерные варианты и особенности экстремальных свойств соответствующих функций с учетом 
сил вязкого трения. 

Обсуждение и заключение. Интерес представляют характерные варианты частотных функций и 
функций демпфирования в зависимости от условий малости сил вязкого трения. Рассматриваются варианты 
значений величины 7,, определяющие характерные интервалы коэффициента формы &@, в которых 
выполняются условия малости сил вязкого трения. 
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1. Пусть 7, Е (0, М,). Рассматривается пример 7, = 0,1. В этом случае выполнены условия малости сил 
трения для любого коэффициента формы &е (-—0..оо). На рис. 4 и 5 представлены частотная функция 6’ (0) и 
функция демпфирования ©(0) для механической упруго-диссипативной системы с параметрами Б=1, т =3, 
т, =4,К=З. 





-0.5- в 
Рис. 4. Частотная функция ©’ (0): Рис. 5. Функция демпфирования ©(0,): 
(1) и (2) — экстремальные уровни (1) и (2) — экстремальные уровни 


Корни уравнения, представляющего собой равенство нулю соответствующего определителя 
Ар? +Вр+С|=0, (39) 
составляют р, =®,+10,, 1=1..4, где в = 0,91; ©, =1,56; ®, =-1,56; в, =-0,91; в =-0,14; в, =-0,44 ; 
с, = -—0,44; с, =-0,14. На рис. 4 частотная функция достигает экстремальных значений, равных квадратам 
частот ©,‘ 22,46 и в* = 0,82. На рис. 5 функция демпфирования достигает экстремальных значений, которые 
составляют ©; =-0,44 и оа,=-0,14. Свои экстремальные значения частотная функция и функция 


демпфирования достигают, когда коэффициенты формы составляют © * = -0,65 и ©,* =1,15. 


т 
2. Пусть у, Е(М,,—). Рассматривается механическая система с параметрами Б=4, 7, =1,33. На 

2 
рис.6 и 7 представлены соответствующие частотная функция и функция демпфирования. Множество 
коэффициентов форм, для которых выполняется условие малых сил трения, составляет: (—00, 0 ) (4,00), где ©, 
— корни уравнения М(о)=у,. Для параметров Р=4, т=3, т, =4, К=З3 корни характеристического 
уравнения (39) имеют вещественные части, представляющие собой коэффициенты диссипации, ©, = -—0,57; 
с, =—1,08; а, =-2,46; <, =-0,57 и мнимые части, представляющие собой частоты, ©, = 0,73; ©, =0; ®, =0; 


> —. 
©, = —0, 73 . На рис. 6 на интервале (—, © )\/(©,,с0) частотная функция имеет локальный минимум в,’ = 0,53. 


у. 
@ [9 


0.6 





(1) 
(2) 
-31 б 
9:1 
Рис. 6. Частотная функция : Рис. 7. Функция демпфирования : (3) -экстремальный 
(1) — экстремальный уровень в точке о. уровень в точке о ; (1), (2) — экстремальные уровни в точке 


* 
О на интервале двузначности 
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На интервале (©,о,) частотная функция равна нулю. В свою очередь, на рис. 7 функция 
демпфирования в интервале (0,0) двузначна и достигает одновременно двух экстремальных значений 
с, =—1,08 и о, =-2,46 в точке а,*. В области (—°,0 )\/(©,,00) функция демпфирования однозначна и 
имеет один локальный экстремум с, = -0,57 в точке @, *. 

3. Пусть уу Е М › ). Рассматривается механическая система с параметрами Б=6; т =3; т, =4; 

2 
К =3. Корни характеристического уравнения (39) имеют вещественные в, =-0,85; в, =-—0,56; в. я 4,74; 
с, =-0,85 и мнимые части ©, = 0,36; ®, =0; ®, =0; <, =-—0,36. Условия малости сил трения выполняются 
на интервале (0,0), где о = 0,63; @, = 2,37. На рис. 8 соответствующая частотная функция положительна 
только на интервале (9,0). Локальный экстремум частотной функции равен 6.” = 0,13. Вне интервала 
(0, %,) частотная функция равна нулю. 

Функция демпфирования на рис. 9 в интервале (—0,0/) двузначна и достигает одновременно двух 
экстремальных значений в точке а, *, которые составляют а, = 4,74 и в, =-0,56. В интервале (&,, оо) 
функция демпфирования также двузначна. В интервале (0,0) функция демпфирования однозначна и имеет 


один локальный экстремум о, = 0,85 в точке ©, *. 





(1) [9 
-4 -2 0 2 4 
[\, 
„1 —ы-{ 
(2) 
_-3 
| (3) 
] 
_-5 б 


в Рис. 9. Функция демпфирования: (1), (3) — экстремальные 


. * 
Рис. 8. Частотная функция: уровни в точке 0, на интервале двузначности, 


* 
вы к * 
(1) — экстремальный уровень в точке 0, (2) — экстремальный уровень в точке (©, 


4. Пусть у, Е (М, ). Для параметров Б = 6,52; т =3; т, =4; К=3 характеристическое уравнение 
(39) имеет только вещественные корни © =-0,50; в, =-0,88; в, =-0,96; в, =-5,26. Соответствующая 


частотная функция и функция демпфирования представлены на рис. 10 и 11. Интервал выполнения условий 
малости сил трения вырождается в пустое множество. 





(2) (1) (1 
> 
у 
в. 
(3) 
=Э 
(4) 

-4 
3 
-5 = 
ы 
(>) 
-6] св х 
Рис. 10. Частотная функция: Рис. 11. Функция демпфирования, образованная двумя ветвями: и 2 

вырожденный случай (1), (2) — экстремальные уровни «верхней» ветви в точках 0 и ©. ; 


* * 
(3), (4) — экстремальные уровни «нижней» ветви в точках 0 и © 
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На рис. 10 частотная функция равна нулю на всей числовой оси. На рис. 11 соответствующая функция 


демпфирования двузначна на всей числовой оси и имеет 4 локальных экстремума о, =-0,50; о, =-0,88; 


нЕ НЕЕ 


= 


* * 
с, 0,96; <, =-—5,26 в точках @ и а, 


(д) нь 








Рис. 12. Пиктограммы топологических особенностей графиков частотной функции и функции демпфирования: 
(а)-(г) — пиктограммы, сопоставленные графикам частотных функций, представленных на рис. 4, 6, 8, 10 соответственно; 
(9)-(3)— пиктограммы, сопоставленные функциям демпфирования на рис. 5, 7, 9, 1] соответственно 


Представленные на рис. 12 пиктограммы, сопоставленные графикам функций на рис. 4—11, отражают 
ряд топологических особенностей графиков частотной функции и функции демпфирования. К таким 
особенностям можно отнести форму графика в виде одной кривой, наличие на графике точек раздваивания 
одной кривой на две, наличие двух непересекающихся кривых или «кольца». 

Таким образом, достигаемые экстремальные значения частотной функции и функции демпфирования 
связаны с динамическими характеристиками механической колебательной системы с учетом сил трения. В 
частности, экстремальные значения построенных частотной функции и функции демпфирования связаны с 
собственными частотами и диссипативными коэффициентами затухающих колебаний. При этом вопрос 
существования экстремальных значений частотной функции и функции демпфирования, не совпадающих с 
квадратами частот собственных колебаний, требует дополнительного рассмотрения. Вместе с тем, можно 
предположить, что формы частотных функций и функций демпфирования, обуславливающие формы 
свободных движений механических колебательных систем с трением, представляют интерес для оценки более 
широкого круга динамических свойств. 

В плане практической реализации возможного управления колебательными режимами механических 
систем на основе коэффициента связности не существует принципиальных препятствий. К примеру, возможные 
системы управления динамическим состоянием могут включать датчики амплитуд колебания в контрольных 
точках вибрационной технологической машины или транспортного средства. Вместе с тем, вопросы 
построения подобных систем требуют детального рассмотрения широкого круга особенностей, связанных с 
техническим объектом. 

В заключение в качестве выводов представленных исследований можно отметить следующие моменты. 

1. Для механической колебательной системы с учетом сил вязкого трения разработан метод построения 
частотной функции и функции демпфирования, зависящих от коэффициента форм связности координат 
свободных движений. Показано, что совокупность экстремальных значений частотной функции и функции 
демпфирования отображает собственные характеристики упруго-диссипативной механической колебательной 
системы. 

2. Показано, что частотная функция и функции демпфирования для механической колебательной с двумя 
степенями системы с учетом вязкого трения может быть представлена в двух вариантах, определяемых условиями на 
величину сил вязкого трения для фиксированного коэффициента форм связности; для условий малых сил вязкого 
трения значения частотной функции принимают положительные значения, а функция демпфирования имеет одну 
отрицательную компоненту; если выполнены условия больших сил вязкого трения, то частотная функция принимает 
нулевые значения, а функция демпфирования имеет две отрицательные компоненты. 

4. Предложен метод построения возможных вариантов частотных функций и функций демпфирования для 
различных значений параметров системы на основе параметризующей функции, позволяющей определить области 
значений коэффициента форм связности, в которых выполняется условие малости сил вязкого трения. Предложен 
критерий классификации частотных функций и функций демпфирования в зависимости от топологических 
особенностей их графиков. 

5. Матричный метод построения частотно-демпирующей функции для системы двумя степенями 
свободы может быть распространен на механические колебательные системы, рассматриваемые в различных 
системах координат. 
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6. В качестве физической интерпретации коэффициента связности, используемого в частотной функции 
и функции демпфирования, может рассматриваться рычажная связь в форме передаточного отношения, 
выраженного через отношения амплитуд колебаний координат парциальных блоков. Рассматриваемое 
отношение, наравне со статическим состоянием, может быть определено для режимов установившихся и 
затухающих колебаний. Таким образом, развивается концепция, в рамках которой отправной точкой анализа 
механической системы является рычажная связь. 
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[#4 
Введение. Рассматривается метод решения задачи о бесконечной пластине, лежащей на упругом основании и 
периодически нагруженной силой, которая перемещается по произвольной замкнутой траектории и 
произвольному закону. 

Материалы и методы. Рассмотрен оригинальный метод решения задач теории упругости для пластин, 
нагруженных силой, движущейся произвольно по замкнутой траектории произвольной формы. Исследуется 
задача о бесконечной пластине, лежащей на упругом основании. Пластина нагружена нормальной силой, 
движущейся с переменной скоростью. Нагрузка раскладывается в ряд Фурье на временном отрезке, длина 
которого равна времени ее прохождения по траектории. Решение указанной задачи реализуется посредством 
суперпозиции решений задач, соответствующих нагрузке, задаваемой слагаемыми указанного ряда Фурье. 
Окончательное решение задачи представляется в виде отрезка ряда Фурье, каждое слагаемое которого 
соответствует решению задачи о действии на бесконечную пластину нагрузки, распределенной по замкнутой 
траектории движения силы. Для построения этих решений используется фундаментальное решение уравнения 
колебания бесконечной пластины, лежащей на упругом основании. 

Результаты исследования. Представлено решение задачи о бесконечной плоскости, по которой с переменной 
скоростью движется сосредоточенная сила. В качестве траектории рассматривалась гладкая замкнутая кривая, 
состоящая из дуг окружностей. Исследуется характер поведения перемещений и напряжений вблизи 
движущейся силы, а также изучается процесс распространения энергии упругих волн. Для этой цели 
рассматривается изменение вектора Умова-Пойтинга. 

Обсуждение и заключения. Полученные результаты могут быть использованы в расчетах при проектировании 
дорог. Исследование распространения энергии упругих волн от движущихся транспортных средств позволит 
оценить воздействие указанных волн на расположенные вблизи дороги строения. Анализ характера изменения 
перемещений и напряжений вблизи движущейся силы позволит оценить износ дорожного покрытия. 


Ключевые слова: бесконечная пластина, движущаяся нагрузка, произвольная замкнутая траектория, 
переменная скорость, энергия упругих волн. 
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Введение. Во многих областях науки и техники возникают задачи, связанные с распространением 
упругих волн. Целью данного исследования является изучение закономерностей распространения упругих 
волн, возникающих под действием подвижной нагрузки. Проблемы подобного рода ранее исследовались 
различными авторами в ряде работ, где рассматривались разнообразные постановки задач с подвижной 
нагрузкой и предлагались различные методы их решения. 

Часто при решении подобных задач для того, чтобы исключить время из числа независимых 
переменных, вводится подвижная система координат [1-5] или рассматривается квазистатическая постановка 
задачи [6-12]. В ряде работ используется метод конечных элементов [11—13]. Интересные результаты позволяет 
получить использование вариационных [14—16] или прямых методов [17-19]. Методы, основанные на 
применении фундаментальных решений соответствующих дифференциальных уравнений, использовались в 
работах [20-22] при решении задач теории упругости о силе, движущейся с постоянной скоростью. 

Постановка задачи. Рассмотрим дифференциальное уравнение, описывающее колебание бесконечной 
пластины, лежащей на упругом основании под действием вертикальной силы Р?: 


Р 
АИ +99 +ЕИ =—, (1) 
р 

ЕЧ” 
где (И — прогиб пластины; ДО = 12а—и2), Е — модуль Юнга; и — коэффициент Пуассона; 4 — толщина 

-Н 

. 2 _рР@. .._К. 

пластины; со = 1 ; р — плотность материала пластины; А = о ; К — коэффициент жесткости упругого 


основания. 

Будем рассматривать решение, для которого поток энергии направлен от источников возбуждения 
упругих волн в бесконечность и считаем, что сила Р перемещается по замкнутой траектории ‘у произвольным 
образом. Кроме того будем полагать, что Р = Р(5(Е)) ‚ где 5 — дуговая координата, отсчитываемая от некоторой 


фиксированной точки кривой у. Очевидно, что Р = Р(5(Р)) будет функцией периодической по Е с периодом Т, 
если 5(Г) также будет функцией периодической по Ес периодом Т. 


Материалы и методы. Рассмотрим фундаментальное решение уравнения (1). Его можно получить из 
выражения: 


АЗ еТОНИ + = би х,)6у-,)8, (9). (2) 


где 0. (1—т) = у 5(Е-т-иГ). 


П=—оо 
Решение уравнения (2) можно получить традиционными методами, применяя принцип предельного 
поглощения, и представить в виде ряда: 


Й’(х, ху, У›Ё-т)= у м», (Х, хо, У, Уз, ®,)е “7? р (3) 


К=Ь—о 


где м’, (х,х,, у, у, ®,) удовлетворяет уравнению: 
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1 | 
А? —с ум, + Ку, = о )5(у-у,)е“". 


Известно, что решение уравнения (1) можно представить в виде: 
Т/2 


И(х, у, в) = р | || И’(х, ху у ,Ё-т)Р(х,, у, тах, ау ах. 


-Т/2 В? 
Если подвижная сила представляет собой единичную сосредоточенную силу, которая описывается 
функцией: 
Р(5(1)) = 9(х — х,(5(1)))0( 7 — У, (5), 


то решение в этом случае будет иметь вид: 


й 
| 2 
Обь) => | бьет. (4) 
Определив ((х, у,г) ‚ можно вычислить перемещения и напряжения в любой точке пластины. 


2х 
Если выражение ©, = Е достаточно велико, то приходиться вычислять интегралы от быстро 


осциллирующих функций (3). Для этого использовалась квадратурная формула [23, 24]: 


и Ь 1 №Ь-—1 9х +1 ох } 
е —е 
| ем (хаха | ге" 5(х)ах= — —т У». В (мы -м 
Е й. › 
а а у е 
где й, — длины элементарных отрезков, на которые разбивается интервал [а; 6]; 5(х) — аппроксимация Дх) 


кубическим сплайном М, =5"(х,). 


Результаты исследования. Расчеты проводились для бесконечной пластины, нагруженной 
нормальной силой, которая двигалась по замкнутой кривой (рис. 1). 


У, м 





Рис. 1. Траектория движения сосредоточенной силы 


При проведении вычислений полагалось, что толщина пластины 4 = 0,25 м, параметр с = 221 м/с, 
модуль Юнга материала пластины ЁЕ= 232469 Н/м”, коэффициент Пуассона и= 0,36, коэффициент 
податливости упругого основания К=1,864 м”. Радиусы, определяющие форму траектории движения силы, 
принимались равными: А! =5м, К, =Зм, №: =1 м. В формуле (3) удерживалось 120 членов ряда Фурье, а при 
вычислении интеграла (4) интервал интегрирования разбивался на 120 равных подинтервалов. 

Закон движения силы по траектории описывался функцией: 


Р-н -7)) в. 


— › где «=>. РЕ 032 |. 


2511(91/) 2 


5(#) = 
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Г ь г ы 
Рассматривался момент времени Е = 5 ‚ В КОТОрыи движущаяся сила находилась в однои и тои же точке 


траектории при любом Т. При изменении параметра Т меняется скоростной режим передвижения 
сосредоточенной силы по траектории. 

Для анализа напряженно-деформированного состояния пластины вычислялись перемещения и 
напряжения в прямоугольной системе координат, связанной с движущейся сосредоточенной силой. При этом 
ось т данной системы направлялась по касательной к траектории т, а ось и совпадала по направлению с 
внешней нормалью к области, ограниченной траекторией п (рис. 1). 

В данной системе координат вектор перемещений и тензор напряжений могут быть представлены в 
виде: 

(0 =(Е-т+0п-п-+й’-К, 
$ = 51. 17+ 5-я +5. (+77), 


где К — нормаль к пластине. 
На рис. 2 представлено изменение по оси т компонент вектора перемещений ИИй,Ииё, тензора 


напряжений 5, эп, 5! (второй и третий графики) и изменение тех же величин вдоль оси и — Ии,Им,Иип и 


эм, опи, эти (первый и четвертый графики) при 2 = — ‚ скорости движения силы у = 49,3480 м/с, ускорении 


2 - 2 
движения по траектории у’, = 0 м/с и нормальном ускорении и’, = И. = 497,045 м/с’ (положение силы на 
1 


траектории отмечено красной точкой на рис. 1). 

Из представленных графиков видно, что значения перемещений Ии и И”значительно превосходят 
значения перемещений в плоскости пластины. То же можно сказать и о нормальных напряжениях, по 
сравнению с касательными. 








У\УМ, ИХ, Ип\, м УТ, Чат, ИТ, м 

0,5 0,5 
0 0 
-0,5 -0,5 
-1 -1 
-1,5 -1,5 
йо -2 
-2,5 -2,5 
ы. -3 

-10 -5 0 ь 10 -10 -5 0 5 10 

ЕТ ЗпТ,ЗиЕТ, н/м? С ЗМ, ие, ним? тм 

10 10 
0 0 
-10 -10 
-20 -20 
-30 -30 
-40 -40 
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-10 -10 -5 0 5 10 





Рис. 2. Изменение перемещений и напряжений (у = 49,3480 м/с) 
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При исследовании распространения энергии упругих волн вычислялся вектор плотности потока 
энергии: 


РЕ-(5 и +50, Я — (52 +5и,)7 , ой. 
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где 5, — компоненты тензора напряжений; ий, — производная по времени координат вектора перемещений. 


На рис.3 представлено распространение энергии упругих волн вблизи сосредоточенной силы, 
положение которой на траектории обозначено красной точкой. Длина изображенного вектора соответствует 
количеству энергии, проходящей через данную точку пространства в единицу времени, а направление вектора 
указывает направление переноса энергии. 


У, м 





-10 -5 0 5 Хм 


Рис. 3. Вектор плотности потока энергии (у = 49,3480 м/с) 


Расчеты показали, что при увеличении скорости движения силы характер изменения перемещений и 
напряжений, а также характер распространения энергии упругих волн несколько изменяется. На рис. 4, 5 
представлены результаты расчетов для скорости у = 246,7401 м/с, которая превосходит скорость распространения 
упругих волн в пластине — с=221 м/с. Влияние скорости движения сосредоточенной силы на характер 


распределения вертикальных перемещений И’ представлено на рис. 6, 7. 


УММ, ИХ, Чп\У, м УТ, Опт, (ТГ, м 
х10 х10? 





-10 -5 0 5 10 -10 -5 0 5 10 
УТ, Зи Т,ЗиеТ, н/м ЗМ, $ М,ЗиеТ, н/м" 





50 50 
0 0 
о -50 
-100 -100 
-150 -150 
-200 -200 
-250 -250 
-10 -5 0 5 10 -10 -5 0 5 10 


х, М 


Рис. 4. Изменение перемещений и напряжений (у = 246,740] м/с) 
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Рис. 5. Вектор плотности потока энергии (у = 246,7401 м/с) 


х107 





У, м 
х10 





У, м -10 


Рис. 7. Изменение вертикальных перемещений при скорости у = 493,480 м/с 
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На рис. 8, 9 представлены графики изменения максимальных значений перемещений и напряжений в 
зависимости от скорости движения сосредоточенной силы. Положение силы на траектории в рассматриваемый 


момент времени отмечено точкой на рис. Зи 5. 


УУМ, УТ, м 
х10° 
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Рис. 8. Изменение максимальных перемещений в зависимости от скорости движения сосредоточенной силы 
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Рис. 9. Изменение максимальных напряжений в зависимости от скорости движения сосредоточенной силы 


На рис. 10 представлено распространение энергии упругих волн около движущейся сосредоточенной 
силы. Касательное ускорение движущейся силы в этот момент равнялось у’, = 1,5503 м/с”. Расчеты проведены 
для момента времени #=7Т при том же, что и в предыдущем случае, законе движения силы вдоль траектории, 
следовательно при том же законе изменения скорости и ускорения движения силы. В этот момент времени сила 
находилась в точке траектории, указанной на рис. 10, а ее скорость равнялась нулю. На рис. 11 изображено 


Ь&р://уезбиК-доп$ва.га 


изменение по оси # компонент вектора перемещений Ил,Ий,ЦиЕ и тензора напряжений 51, 5ИЕ, 5Е. 


Увеличение ускорения движения силы также приводило к изменению перемещений и напряжений и 
влияло на характер распространения энергии упругих волн. На рис. 12, 13 представлены результаты расчетов 


376 для случая и = -155,0314 м/с*. 
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Рис. 10. Вектор плотности потока энергии (у», = 1,5503 м/с?) 
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Рис. 11. Изменение перемещений и напряжений (и, =1,5503 м/с”) 
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Рис. 12. Изменение перемещений и напряжений (и, = 155,0314 м/с”) 
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Рис. 13. Вектор плотности потока энергии (у›, = 155,0314 м/с”) 


378 


Галабурдин А. В. Бесконечная пластина, нагруженная нормальной силой, движущейся по сложной траектории 





На рис. 14, 15 представлены графики изменения максимальных значений перемещений и напряжений в 
зависимости от касательного ускорения движения сосредоточенной силы. 
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Рис. 14. Изменение максимальных значений перемещений в зависимости от касательного ускорения 
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Рис. 15. Изменение максимальных значений напряжений в зависимости от касательного ускорения 
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Обсуждение и заключения. Имеет место выраженная зависимость полей перемещений и напряжений 
от скорости и ускорения движения силы при рассмотренных выше пределах изменения этих параметров. 
Заметно зависит от скорости и ускорения также и характер распространения энергии упругих волн. 
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Были специально рассмотрены достаточно большие значения скорости и ускорения движения силы для 
того, чтобы проверить предложенный метод в подобных условиях. Полученные результаты позволяют сделать 
вывод о том, что метод достаточно устойчив в широких диапазонах варьируемых параметров. 

Применение предложенного метода вполне допустимо при решении и более сложных задач. Для этого 
необходимо, чтобы описывающие их дифференциальные уравнения допускали аналитическое построение 
фундаментального решения. Данный метод отличается экономичностью и простотой, так как использует для 
построения решения уже известные решения задач. 
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Обоснование гранулометрических характеристик рабочей среды 
при вибрационной обработке деталей с малыми пазами и отверстиями 
М.А. Тамаркин', Е. Н. Колганова', М. А. Ягмуров” 
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* ФГАОУ ВО «Северо-Кавказский федеральный университет» (г. Ставрополь, Российская Федерация) 





Введение. Представлена методика проектирования высокоэффективного технологического процесса вибрационной 
отделочной обработки деталей с малыми пазами и отверстиями. Решение основано на выборе гранулометрических 
характеристик обрабатывающих сред. Проанализирована форма поперечного сечения и геометрические размеры 
заусенцев на типовых деталях радиоэлектронной аппаратуры (РЭА). Разработана обобщенная модель заусенца. 
Определены методологические принципы выбора гранулометрических характеристик рабочих сред. 

Материалы и методы. Новая классификация и кодирование по конструктивно-технологическим принципам 
деталей РЭА позволят при проектировании технологического процесса их финишной обработки обоснованно 
выбрать оборудование, среды и режимы. 

Результаты исследования. Создана методика подбора гранулометрических характеристик рабочих сред с учетом 
основных технологических задач. Исходя из конструктивно-технологических особенностей деталей РЭА, предло- 
жены зависимости для определения размеров и формы обрабатывающих гранул. Определены критерии годности 
для оценки результатов вибрационной обработки. Отмечено, что одна из главных задач при виброобработке дета- 
лей с малыми пазами и отверстиями — обеспечение такой продолжительности процесса, при которой удаляются 
заусенцы, а шероховатость и другие параметры поверхности соответствуют техническим требованиям. В таком 
случае показателем качества следует считать точность линейных размеров обрабатываемых поверхностей. Коли- 
чественно данный критерий оценивается по специальному индексу, при расчете которого учитывается наибольший 
действительный размер до виброобработки, высота заусенца, наименьший допустимый размер после обработки и 
допуск, установленный техническими требованиями. Критерий эффективности процесса определяется как отно- 
шение индекса обрабатываемости к продолжительности обработки партии деталей или длительности цикла, при- 
веденной к одной детали. Предложенный критерий позволяет сравнивать процессы обработки при обосновании 
решения технологических задач. 

Обсуждение и заключения. Итоги исследования позволяют утверждать, что вибрационная обработка в рабочих 
средах органического происхождения способствует эффективному удалению заусенцев и скруглению кромок 
малоразмерных деталей радиоэлектронной аппаратуры, имеющих малые пазы и отверстия. 


Ключевые слова: вибрационная обработка, шероховатость, микронеровность поверхности, заусенцы, скругление 
кромок. 
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Введение. Заметно возросшая в последнее время потребность в прецизионной обработке деталей особенно 
актуальна в машино- и приборостроении, т. к. в этих отраслях годовые объемы изготавливаемых деталей превы- 
шают сотни тысяч штук. Следует особо подчеркнуть, что детали радиоэлектронной аппаратуры (РЭА) имеют до- 
статочно сложную конфигурацию наружного контура. Для большинства из них характерны такие нетехнологич- 
ные элементы, как пазы и отверстия малых размеров, глубокие отверстия, глухие отверстия с резьбой и т. п. Под- 
твердила свою эффективность финишная обработка подобных деталей в гранулированных средах [1, 2]. В настоя- 
щее время отсутствуют методики для проектирования отделочно-зачистных операций таких деталей, что ограни- 
чивает широкое внедрение и дальнейшее совершенствование вибрационной обработки. 

Материалы и методы. На основе конструкторско-технологического анализа деталей РЭА, изготавливае- 
мых на приборостроительном заводе, разработана классификация и выполнено их кодирование по конструктивно- 
технологическим принципам. Новый классификатор позволит на стадии проектирования технологического про- 
цесса финишной обработки деталей с учетом их конструктивно-технологических особенностей обоснованно вы- 
брать оборудование, среды и режимы операций. 

Особую роль в обеспечении качества деталей РЭА играет отделочно-зачистная обработка, выполнение ко- 
торой в условиях современного производства затруднено сложной конструктивной формой изделий. 

Одним из главных методов формообразования заготовок в машиностроении и приборостроении остается 
лезвийная обработка. Однако при любой механической обработке, основанной на резании, образуются такие де- 
фекты, как заусенцы, острые кромки и пр. Они возникают из-за действия законов механики сплошных сред, по- 
этому указанные недостатки невозможно исключить даже при использовании современных обрабатывающих цен- 
тров и оптимальных режимов обработки. В результате при изготовлении и эксплуатации деталей возможны функ- 
циональные, эстетические и эргономические проблемы [3, 4]. Это доказывает актуальность проблемы удаления 
заусенцев при отделочно-зачистных операциях деталей РЭА. 


Машиностроение и машиноведение 


383 


ЬИр://уезш-допза.га 


384 


А4уапсеа Епатеетиз Кезеагсй 2020. Т. 20, №4. С. 382-389. 155М 2687-1653 





Цель данной научной работы — совершенствование методики проектирования высокоэффективных тех- 
нологических процессов вибрационной обработки деталей с малыми пазами и отверстиями на основе подбора гра- 
нулометрических характеристик гибких рабочих сред. 

Для достижения поставленной цели необходимо решить следующие задачи: 

— проанализировать формы поперечных сечений и геометрических характеристик заусенцев на типовых 
деталях приборов, 

— разработать обобщенную модель заусенца, 

— обосновать методологические подходы при выборе гранулометрических характеристик обрабатываю- 
щих сред. 

Для изучения геометрии и параметров профиля сечения заусенцев микрошлифы рассматривали с помо- 
щью металлографического инвертированного микроскопа, оборудованного системой Тыхотей Рго. Как показало 
исследование, основная форма сечение заусенца — треугольник. При этом высота заусенца в среднем больше 
толщины его основания в 2-3 раза, а линейные размеры толщин не превышают 0,4 мм (рис. 1). 





в) 
Рис. 1. Сечение заусенцев у исследуемых образцов. Материал образцов: БрОЦ4-3 (а), ЛС-59-1 (6), АМгб (в) 


Выбор характеристик обрабатывающей среды в данном случае оказывает решающее влияние на качество 
поверхности деталей и производительность процесса [5-7]. На предприятиях применяют рабочие среды, гранулы 
которых можно классифицировать по геометрической форме и размерам, материалу и размеру режущих зерен, 
материалу связки зерен, структуре и способу производства. 

Для удобства проектирования технологических процессов отделочно-зачистной вибрационной обработки 
предложены системная классификация и кодирование гранул рабочих сред. 

Отделочно-зачистная обработка деталей РЭА предусматривает удаление заусенцев и подготовку поверх- 
ностей под покрытия. В приборостроении применяются различные покрытия: цинкование, кадмирование, никели- 
рование, хромирование, латунирование, палладирование, серебрение, золочение, пассивирование и т. д. На участки 
поверхностей электроконтактных деталей приборов наносят гальванические покрытия (серебро, золото), которые 
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характеризуются высокой электро-, теплопроводностью и химической стойкостью в условиях повышенной влаж- 
ности [8, 9]. Гальваника требует определенного качества заготовки. Не допускается наличие заусенцев и острых 
кромок. 

Как показал анализ конструкций, при отделочно-зачистной обработке деталей РЭА следует использовать 
среды, соответствующие перечисленным ниже требованиям: 

— высокая плотность (не менее 1,2 г/см^) при малом весе гранулы; 

— пригодность для обработки нежестких заготовок; 

— пригодность для обработки заготовок, имеющих сопряженные под углом поверхности и зоны с ограниченным 
доступом для гранулы обрабатывающей среды; 

— высокая износостойкость и способность сохранять форму в процессе обработки; 

— равномерная структура гранул. 

Учитывая эти требования и данные работ [10, 11], следует признать наиболее целесообразным использо- 
вание косточковых гранулированных сред. При дроблении косточек плодовых культур и скорлупы грецкого ореха 
образуются гранулы произвольных форм с краями в виде клина, благодаря которому гранула становится аналогом 
режущего инструмента. В этом случае для микрорезания доступны различные участки детали. 

Результаты исследования. Итак, среды из природных материалов представляются оптимальными для от- 
делочно-зачистной обработки деталей РЭА с малыми пазами и отверстиями. Широкому применению данного под- 
хода препятствует недостаточная изученность таких сред и недоступность методик проектирования отделочно- 
зачистных операций с их использованием. 

Как правило, при отделочно-зачистной вибрационной обработке деталей РЭА возникают три технологиче- 
ские задачи. 

Первая. Если на заготовках достаточно удалить заусенцы и скруглить кромки, нет необходимости обраба- 
тывать внутренние поверхности пазов и отверстий. Это наиболее характерная задача, решаемая в процессе обра- 


ботки деталей. В общем виде ее можно решить, применяя гранулы с размером А;,, превышающим размер самого 


гр> 
большого отверстия или паза [,, т. е. К’, > Ё. Это исключит заклинивание в них. 

Контакт обрабатывающей среды будем рассматривать как контакт единичной гранулы до полного удале- 
ния заусенца и получения требуемой величины скругления кромки (рис. 2). 














Рис. 2. Геометрическая схема контакта гранулы с кромками отверстия 


Учтем коэффициент износа гранулы К„, во многом определяющий расход обрабатывающей среды, а также 
производительность процесса и качество поверхностей. Оптимальный радиус гранулы рабочей среды определяется 
выражением: 

(+213)? +1 


2 


где й, — толщина основания сформированного заусенца, мм; йу — размер фаски, соответствующий техническим 
требованиям, мм; Г, — линейная или диаметральная величина обрабатываемого элемента, мм. 

Соблюдение данного условия гарантирует съем заусенцев по наружным поверхностям и обеспечение тре- 
буемого чертежом радиуса скругления кромок. 


Машиностроение и машиноведение 


385 


Ь&р://уезбиК-доп$ва.га 


386 


Айуапсей Епатеетиз Кезеагсй 2020. Т. 20, № 4. С. 382-389. 155М 2687-1653 





Вторая. Эта задача возникает при необходимости обработки и подготовки под покрытие внутренних по- 
верхностей пазов или отверстий. В этом случае критерием для выбора размера частиц должен быть размер 
наименьшего отверстия или паза Ги; детали. 

Для исключения заклинивания частиц рабочей среды в отверстиях и пазах их размер рекомендуется выби- 
рать из соотношения 

Пр = 0,6 + 0,7Ё т, 
где О; — характерный диаметр гранул обрабатывающей среды, мм. 

Однако при выборе размера частиц рабочей среды с использованием этого соотношения крайне важно 
проанализировать его применимость для обработки других отверстий и пазов детали, размер которых превышает 


Ё тт. Если их размер равен 20;, или З);,, то возможно заклинивание частиц рабочей среды, которое будет преры- 


гр› 
вать процесс обработки поверхности отверстия или паза. Для предотвращения этого явления размеры подлежащих 


обработке отверстий и пазов [ должны быть в диапазоне Риш < Ё < 1,7Б,,, 2,2, < Ё < 2,70, Ё > 3,2, 


(рис. 3). 


гр› 


гр> 











Рис. 3. Размеры отверстий, исключающих заклинивание тел рабочей среды: 
Ар — рекомендуемый размер, Ан — недопустимый размер 


Третья. Основная технологическая задача — обработка мест сопряжения поверхностей. Ее решение вызы- 
вает наибольшие сложности, т. к. гранулам сложно достичь поверхностей, находящихся под углом друг к другу. В 
качестве примера можно привести резьбу. Если в технических требованиях детали задан радиус сопряжения по- 
верхностей, то размер гранулы рабочей среды должен быть равен ему либо меньше: 

Ар 5 Г. 

Как показано в [12—14], из-за нарушения данного условия при вибрационной обработке сопряженных под 
углом поверхностей могут образоваться три условные зоны: 

— мертвая зона (в ней не происходит обработка), 

— открытая зона, 

— зона нестабильной шероховатости (рис. 4). 





Рис. 4. Не обрабатываемые зоны, возникающие при неправильном выборе размеров гранул 


Размеры мертвой зоны и зоны с нестабильной шероховатостью можно рассчитать по формулам: 


1 
гр 1800—а 
со5=——— 
2 
1800—а 
СВ = Вр ап“, 
2 


Для обработки труднодоступных поверхностей органическими частицами с клиновидными вершинами 


СО < В 1]; 


можно ввести понятие коэффициента проникаемости гранулы, который определяется из соотношения угла сопря- 
жения поверхностей к углу клиновидного ребра, образующего режущую кромку: 
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В. 
Учет коэффициента проникаемости при выборе формы гранул позволит решить сложную технологиче- 


Кир > 


скую задачу, связанную с отделочной обработкой мест сопряжений поверхностей. 

Для контроля результата отделочно-зачистной обработки необходимо сформулировать критерии годности 
деталей или их партий. Для деталей РЭА такими критериями будут параметры шероховатости поверхности и ста- 
бильность качественных характеристик обработанной поверхности (рис. 5). 
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Рис. 5. Определение индекса обрабатываемости: 1 — шероховатость поверхности после обработки; 2 — шероховатость поверхности 
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действительная; 3 — заусенец; ЕЙ’ — отклонение волнистости; ЕЁ — отклонение формы; Ам. — наибольший размер до 
обработки с учетом высоты заусенца; А„„ — наименьший размер после обработки; ТА — допуск на обрабатываемый размер 


Главное условие для решения данной технологической задачи: 
К, —6> [пт 
где А, — среднее статистическое показание величин значений принятых показателей качества, которые соответ- 
ствуют среднеарифметическим отклонениям профиля исходной поверхности; 6 — поле рассеивания значений ка- 
чественных показателей, которые соответствуют величинам съема металла при обработке; [В;„ш| — величина 
значения показателей качества в пределах допуска. 

Показатели качества в пределах партии деталей будут случайными величинами, распределенными в соот- 
ветствии с вероятностными законами. Поэтому критерии годности партии (т. е. оптимальности предложенной тех- 
нологии) можно задать предложенным выше неравенством. При этом нужно учитывать процент брака, не превы- 
шающий вероятности РТ [15]: 

В (Ви Е > 6) <=. 
где Р; — вероятность попадания 1-го показателя качества за пределы границ поля рассеяния; /); — расчетное зна- 
чение дисперсии К;. 

Одна из главных задач при виброобработке деталей с малыми пазами и отверстиями — обеспечение такой 
продолжительности процесса, при которой удаляются заусенцы, а шероховатость и другие параметры поверхности 
соответствуют техническим требованиям. Показателем качества будет точность линейных размеров обрабатывае- 
мых поверхностей. Количественно данный критерий оценим по индексу обрабатываемости (рис. 5): 


а — Атах Ани, 
ТА 
где А„„х — наибольший действительный размер до виброобработки, учитывающий высоту заусенца; Аш„ — 
наименьший допустимый размер после обработки; ТА — допуск на обрабатываемый размер, установленный тех- 
ническими требованиями (рис. 5). 
Критерий эффективности процесса определим как отношение индекса обрабатываемости к продолжитель- 


ности обработки партии деталей или длительности цикла, приведенной к одной детали: 


Со 
К. и Ри. 
Ц 


Предложенный критерий позволяет сравнивать процессы обработки при обосновании решения технологи- 
ческих задач. 

Обсуждение и заключения. В работе проанализированы форма поперечного сечения и геометрические 
размеры заусенцев типовых деталей приборов РЭА. Создана обобщенная модель заусенца, предложены методоло- 
гические принципы выбора характеристик гранул рабочих сред. 

Описаны способы подбора гранулометрических характеристик обрабатывающих сред в зависимости от 
главных технологических задач. Получены зависимости для определения размеров и формы обрабатывающих гра- 
нул, основанные на конструктивно-технологических особенностях деталей РЭА. Предложены критерии годности 
для оценки результатов вибрационной обработки. 
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Итоги исследования позволяют утверждать, что отделочно-зачистная вибрационная обработка в среде из 
дробленых косточек способствует эффективному съему ликвидов, в т. ч. заусенцев, а также скруглению кромок у 
деталей радиоэлектронной аппаратуры сложной конфигурации. 
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МАШИНОСТРОЕНИЕ И МАШИНОВЕДЕНИЕ 
МАСШИЛЕ ВОПО1УС АХО МАСНТУЕ ЭСТЕМХСЕ 








УДК 621 6 рз://4о1.0го/10.23947/2687-1653-2020-20-4-390-396 
Технологические особенности упрочнения коленчатых валов ОНО 
виброударным методом 


В. А. Лебедев', Ф. А. Пастухов', М. М. Чаава', Г. В. Серга [=] ый 
'ФГБОУ ВО «Донской государственный технический университет» (г. Ростов-на-Дону, Российская Федерация) 
*ФГБОУ ВО «Кубанский государственный аграрный университет» (г. Краснодар, Российская Федерация) 





Введение. Рассматриваются технологические особенности обработки коленчатых валов виброударным методом 
поверхностного пластического деформирования (ШПД), получившего широкое применение в технологии 
изготовления деталей машин. Целью исследований является обоснование эффективности влияния виброударной 
упрочняющей обработки на повышение качества и эксплуатационных характеристик коленчатых валов (КВ). 
Материалы и методы. Методически исследования включали в себя обоснование технологической схемы 
виброударной обработки и разработку расчётной модели оценки влияния обработки на изменение макрогеометрии 
(коробления) КВ. 

Результаты исследования. Разработаны технологические схемы объемной виброударной отделочно-упрочняющей 
обработки КВ с применением вибростанка с Ц-образной рабочей камерой, габаритные размеры которой соизмеримы 
или превышают габаритные размеры подвергаемого обработке КВ и обеспечивают расположение вала таким 
образом, чтобы его основная ось, совпадающая с осью коренных шеек, находилась в зоне расположения условной 
оси вращения массы рабочей среды. Исследованы параметры качества поверхности при обработке их на 
виброустановке УВГ 4Х10 по апробированным методикам с применением специальной оснастки. Установлено, что 
виброударная упрочняющая обработка (ВиУО) позволяет за счет пластической деформации микронеровностей 
получить качественно новый микрорельеф поверхности и снизить ее исходную шероховатость, существенно 
повышает микротвердость поверхности коренных и шатунных шеек КВ и изменяет напряженное состояние их 
поверхностного слоя. Предложена расчётная зависимость для оценки суммарного коробления КВ, упрочнённого в 
процессе ВиУО и подтверждена еб адекватность. Показано, что коробление вала после ВиУО обусловлено 
разнонапряжённостью шатунных и коренных шеек КВ на уровне К; = 0,6. 

Обсуждение и заключение. Виброударная обработка КВ обеспечивает улучшение геометрических и физико- 
механических параметров поверхностей шатунных и коренных шеек. В результате обработки всех поверхностей 
вала коробление не превышает допустимых значений, устанавливаемые техническими требованиями. Это позволяет 
сделать вывод 06 эффективности рассматриваемого способа упрочнения КВ с целю повышения их 
эксплуатационных свойств. 


Ключевые слова: коленчатый вал, поверхностное пластическое деформирование, виброударный метод, упрочнение, 
качество поверхности, коробление. 


Для цитирования: Технологические особенности упрочнения коленчатых валов виброударным методом / 
В. А. Лебедев, Ф. А. Пастухов, М. М. Чаава, Г. В. Серга // Адуапсеа Епошеегие Везеагсв. — 2020. — Т. 20, № 4. 


— С. 390-396. № рз://401.0ге/10.23947/2687-1653-2020-20-4-390-396 
© Лебедев В. А., Пастухов Ф. А., Чаава М. М., Серга Г. В., 2020 





Тесппо10о1са| Геабиге$ оГ сгапК5Вай Вагдепто Бу уфгайоп зВоск тео4 


У. А. Г.еЪедеу', Е. А. РазбикКПоу’, М. М. СПаауа', С. У. бегоа” 


'Роп З{ае Тесшиса! Ошуегзйу (Возюу-оп-Ооп, Визап Еедеганоп) 
`Кибап Зе Аэтамап Ошуегзйу (Ктазподаг, Вазз1ап Еедеганоп) 


[итодисноп. Те ф1есвио]оз1са| Геаагез оЁ Ше ргосеззше стапкзВаЙ$ Бу Фе угаНоп зВоск те фо4 оф за! асе р1азИс 
деГогтайоп (ЗРО), \ВлеЬ 1$ ул4ейу чзе4 ш Фе 1есвпо]осу оГ тапасвигше тасбше ра, аге сопз14еге4. ТВе 


Лебедев В. А. и др. Технологические особенности упрочнения коленчатых валов виброударным методом 





гезеагсь об]есйуе 1$ о ра5ЯБу Фе еЁсепсу оЁ Фе шйЙчепсе о{ уБгаНоп воск Ваг4епие фгеайтеп оп пиргоуте Фе 
ацау ап4 ре{огтапсе о? сгапКзВай$ (СЗ). 

Маепаб апа Мефос4$. Тбе тефо4о]огл1са| за ез шс[аде4 е уаП4айоп о? Ве убгаНоп Воск ргосеззте Но\усваи 
апа Фе 4еуе|ортеп ог ап апайуйс то4е| Гог аззеззте фе еНес® оЁ ргосеззте оп фе свапое ш фе тасгосеотену 
(\уаграее) ое СЪ. 

Кезий5. Ио\усВаг$ Вауе Бееп 4еуе]оре4 юг уо[атейлс уганоп звоск Ни1$Вте ап4 Вагдепие фгеайтеп ог С$ ие 
а угаюгу тасбше ми а О-5ВареЯ у’огК свашфег. [5 оуегаЙ аппеп$10и$ теазиге аПКе ог ехсее4 Фе оуегаП 
Аптепз1оп$ оЁ Фе СЗ Бешс ргосеззеа, ап епзиге 1е 1осайоп о{ фе зваЁ ш зисВ а \уау ай 1$ тат ах1$, сошс1те 
УМЕ Ще ах15 оР Фе таш Беагис ]огпа[$, 15 ш Ще топе оЁ Фе сопд опа! ах1$ о уогкше та$$ гобайоп. ТБе зигЁРасе 
ацаШу рагатаегз \еге шуезисае дигие Фет ргосеззш® оп Фе ЦУС 4Х10 угаюгу ип ассогашо фо Фе ргоуеп 
{есби1ачез изше Фе 4ефсае4 1юоПпе. [ Ваз Бееп езабИзВе Фа{ Фе уганоп зВосК Вагд4епте ‘геайтете (У1НТ) 
епа ]е$, Чие © р!азис деогтаНоп о? пл1сгогоце пез, 0 оМаш а дааШайнуе[у пе\и зигЁасе писгогепе{ ап © гедисе 
1$ па| гоч2рпез$, ю ЧгазисаПу шстеазе бе зигРасе писговаг4пез$ ог Фе С$ таш ап4 то4 ]оигпа[$; аё а и 
свапоез Ше згеззе зе оГ пет зигЁасе Тауег. А са1си]ме4 дерепдепсе 1$ ргорозе4 ю аз$ез$ пе ю{а| \уагразе оЁ е 
СЗ збепоепе ипаег Фе УНТ, ап4 15 адеаиасу 15 сопйгтеа. [ 15 зВо\’й фа Фе угаграсе оЁ Фе зай айег Фе 
УНТ 15 дие ю Фе аШегеп (епзлоп ое го4 ап4 тата ]огпа[$ о Фе СЗ а Фе 1еуе] оЁ К‚ = 0.6. 

О5си5яюп апа Сопсияоп5. Утгайоп Воск фгеайтепт ог СЗ ргоу!14ез ап пиргоуетептё ш Фе сеотейлс апа 
рвуз1сотесватса| рагатеег$ оЁ загЁасез оЁ Ше го4 ап та ]оигпа15. Аз а геза ог ргосеззше аП зигРасез оЁ Фе 
вай, ууаграсе 4оез по{ ехсее4 Ше регил1$51Ше уаез езба Вед Бу Ше 1есби1са1 гедитетеп 5. 50, \е сап сопсшае 
оп Фе еЁйслепсу оЁ Фе сопз14егеа тео4 о{ Вагдетите СЗ уу Фе апп оЁ шсгеазше ег орегайопа] ргорегаез. 


Кеунога$: стапкзВай, за асе р]азис деогтайоп, угайоп зВоск тево4д, Ваг4епто, зиГРасе диаЩу, \агразе. 


Еог сйапоп: У. А. ГеБедеу, Е. А. РазбаКВоу, М. М. СБаауа, С. У. Зегоа. Тесвпо]ое1са| Теаагез оЁ сгапкзВай 
Багдепие Бу угайоп Воск штефо4. Адуапседа Епошеегие КВезеагсв, 2020, уо1|. 20, по.4, рр. 390-396. 


№рз://401.0го/10.23947/2687-1653-2020-20-4-390-396 


Введение. Комплекс мер по сохранению точности КВ, достигнутой формообразующей обработкой, 
включает совершенствование способов получения заготовок, механической обработки, а также введение в 
общий технологический процесс операций упрочняющей обработки, предназначенных для повышения 
эксплуатационных свойств этих изделий. Одними из распространённых способов ШПД, упрочняющих галтели 
и повышающих сопротивление усталости КВ, являются обкатка и чеканка [1-7]. Однако упрочнение КВ 
методами обкатки и чеканки галтелей сопровождается их короблением, которое приводит к увеличению биения 
коренных шеек, нарушению исходной геометрической формы, что требует проведения дополнительной 
операции правки. В этой связи разработка новых технологических способов отделочно-упрочняющей 
обработки методами ППД, повышающих надёжность и долговечность КВ в сложных условиях эксплуатации, 
является весьма актуальной. Особый интерес в этом направлении представляет вибрационная ударная 
обработка, получившая широкое практическое применение в технологии изготовления деталей машин [8—15]. В 
этой связи целью исследований является обоснование эффективности влияния виброударной упрочняющей 
обработки на качество и эксплуатационных характеристик КВ. 

Материалы и методы. Для достижения поставленной цели поставлены следующие задачи: 

— разработать технологические схемы виброударного упрочнения всех элементов КВ; 

—Щ исследовать влияние виброударной обработки на геометрические и физико-механические 
параметры поверхности шатунных и коренных шеек; 

— предложить расчётную модель оценки изменения макрогеометрии (коробления) КВ в процессе 
виброударной обработки и подтвердить еб адекватность. 

Исследования проводились на виброустановке УВГ4Х10 по апробированным методикам. Для 
определения параметров качества упрочнённой поверхности применялось специальное приспособление, 
имитирующее коленчатый вал. Результаты измерений обрабатывались методами математической статистики. 
Коробление исследуемых валов определялось путем замеров радиального биения шеек до и после упрочнения 
индикатором с ценой деления 0,01 мм. При этом КВ устанавливался коренными шейками на призматические 
опоры и поворачивался от руки. 

Результаты исследования. На рис. | представлены технологические схемы объемной виброударной 
отделочно-упрочняющей обработки КВ на базе вибростанка с Ч-образной рабочей камерой, габаритные 
размеры которой соизмеримы или превышают габаритные размеры подвергаемого обработке изделия. Вал в 
процессе обработки устанавливается таким образом, чтобы его основная ось, совпадающая с осью коренных 
шеек, находилась в зоне расположения условной оси вращения рабочей среды. В результате такого 
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расположения вала исключается центральная зона малой активности, а смещённые относительно основной оси 
вала цилиндрические поверхности коренных и шатунных шеек подвергаются обработке в зонах среднего и 
максимального давления. КВ устанавливается в ложемент по направляющим пазам и удерживается на опорах, 
прикреплённых к нему. Ложемент с валом погружаются в камеру с рабочей средой и фиксируется на стойках, 
прикреплённых к станине. Опоры позволяют валу прокручиваться вокруг своей оси и равномерно упрочняться 
под динамическим воздействием рабочей среды (рис. 1 а-6) или путём сообщения ему встречного или 
попутного вращательного движения с помощью автономного привода (рис. | в) 
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Рис. 1. Технологические схемы обработки КВ с вращением под действием рабочей среды (а, 6) и с помощью 
дополнительного привода (в): 1 — станина; 2 — пружины; 3 — рабочая камера; 4 — ложемент; 5 — коленчатый вал; 
6 — рабочая среда; 7 — вибратор, 8 — дополнительный привод 


Для определения предпочтительных режимов обработки КВ выполнены исследования геометрических 
и физико-механических параметров поверхности шатунных и коренных шеек с помощью цилиндрических и 
кольцевых образцов из сталей 45 и 40Х с исходной шероховатостью К. = 0,16-0,42 мкм. Обработка 
осуществлялась рабочей средой, состоящей из стальных шариков диаметром 3—6 мм, по базовой 
технологической схеме, приведённой на рис. 1 а-6 с применением специальной оправки (рис. 2) при различных 
амплитудно-частотных характеристиках и продолжительности процесса. 





а) 6) 
Рис. 2. Экспериментальная оправка, имитирующая коленчатый вал: 
а — общий вид; 6 — монтаж в рабочей камере вибрационной установки 


Установлено, что в процессе ВиУО на поверхности образцов формируется качественно новый 
микрорельеф, среднеарифметическое отклонение профиля которого меньше, чем на исходной поверхности. На 
изменение среднеарифметического отклонения профиля наиболее существенное влияние оказывают амплитуда 
колебаний рабочей камеры и продолжительность обработки. Так за один и тот же промежуток времени 
обработки 1=20 мин. при частоте колебаний !=25Гц увеличение амплитуды от 2 до Змм 
среднеарифметическое отклонение профиля поверхности понизилось в 1,7 раза, а изменение частоты колебаний 
от 20 до 30 Гц при амплитуде А = 3 мм понизило среднеарифметическое отклонение профиля поверхности в 1,1 
раза. 

С учетом технологических возможностей вибростанков установленные закономерности позволили 
обосновать амплитудно-частотные параметры, обеспечивающие наиболее эффективное энергосиловое 
воздействие на деталь в процессе ВиУО на уровне А = 3 ммиЕ= 25 Гц. Обработка образцов на этих режимах в 
течение 20 мин, как следует из рис. 3, обеспечивает снижение среднеарифметического отклонения профиля 
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исходной поверхности в 2,6 раза или на 60%. Увеличение продолжительности обработки до 40 мин к 
существенному изменению рассматриваемого параметра не приводит. 
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Рис. 3. Зависимость среднеарифметического отклонения профиля поверхности 
от продолжительности виброударной обработки 


Подтверждением эффективности выбранных режимов являются результаты исследования физико- 
механических характеристик поверхностного слоя после обработки стальными шариками в режиме А = 3 мм, 
= 25 Гц (рис. 4, 5). Оценка микротвёрдости проводилась на твердомере ПМТЗ, напряжённое состояние 
оценивалось по величине остаточных сжимающих напряжений, определяемых по методике Давиденкова путём 
разрезки кольцевых образцов. 
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Рис. 4. Влияние продолжительности обработки на микротвёрдость поверхностного слоя для материала шариков: 
| — сталь 40Х; 2 — сталь 45 
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Рис. 5. Зависимость остаточных сжимающих напряжений в. от продолжительности ВиУО образцов, 
установленных на втулках (1) и на валу (2) 
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Анализ напряжённого состояния образцов, упрочнённых в процессе ВиУО (рис. 5), показал, что 
остаточные сжимающие напряжения образцов, закреплённых на втулках и имитирующих шатунные шейки КВ, 
на 35 % превышает таковые образцов, закреплённых на валу и имитирующих коренные шейки. Это 
обусловлено различием интенсивности воздействия обрабатывающей среды в различных зонах рабочей 
камеры. 

Следующий этап исследований предусматривал разработку расчётной модели оценки изменения 
макрогеометрии (коробления) КВ в процессе виброударной обработки и экспериментальное подтверждение её 
адекватности. Емельяновым В. Н. [1] предложен графоаналитический метод анализа макродефомации КВ после 
упрочняющей обработки методом ППД (рис. 6). 
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Рис. 6. Схема коробления КВ после ППД галтелей коренных и шатунных шеек 


Суть метода заключается в определении смещения конца вала дд в сечении АА в зависимости от его 
геометрических размеров: 7 — расстояние от оси коренной шейки (КШ) до оси шатунной шейки (ППП); й — 
толщина щеки кривошипа; 5 — расстояние от оси КШ до конца КВ; [, — расстояние между осями КШ одного 
кривошипа; С — расстояние от щеки кривошипа со стороны ШШ до оси КШ. Смещение дд соответствует углу 
отклонения оси вала от горизонтали, величина которого зависит от углового положения щеки кривошипа 0, 
которая жёстко связана с коренной и шатунной шейками. Взаимосвязанные параметры бд иос@ зависят от 
средней величина остаточных напряжений сжатия о. в поверхностном слое после ПД и глубины их залегания 
бс. Кроме того, в работе [1] показано, что величину коробления многокривошипного КВ в основном определяет 
степень упрочнения галтелей КШ и ШШ, прилагающих к крайней щеке крайнего кривошипа. 

С использованием рассмотренного метода получено выражение для расчёта суммарного коробления 


__ бы т ЬБ 


КВ, упрочнённого ВиУО: 
ОА Е -7'т' Ку" Одт " Ккоз [-Кн(Ё + — С) т (Е - С)], (1) 


где Е — модуль упругости первого рода; у — коэффициент Пуассона; от — динамический предел текучести; 
Кко2 — коэффициент корректировки величины остаточных сжимающих напряжений ШШ в зависимости от их 
диаметральных размеров и деформационных параметров энергосилового воздействия части рабочей среды; 
К, = 1,1-1,5 — коэффициент корректировки глубины залегания остаточных сжимающих напряжений; К, — 


коэффициент, учитывающий отличие напряженного состояния КШ от напряженного состояния ШШ. В 
соответствии с результатами настоящих исследований этот коэффициент зависит от удаления от стенок 
рабочей камеры и составляет К’, = 0,6. 

Для подтверждения адекватности модели (1) проведены экспериментальные исследования на 5 натурных 
КВ, прошедших полную механическую обработку. Материал коленчатого вала — сталь 45, твёрдость после 
отжига 180-228 НВ, шатунные шейки диаметром 25 мм подвергнуты закалке ТВЧ на глубину 2—4 мм до 
твердости 52 — 65 НВС, при этом галтели радиусом 2 - 3,2 мм с шероховатостью А, = 1,6 мкм остаются без 


термообработки. Механические характеристики материала КВ: у = 0,25; Е = 2.10”МПа; о. = 360 МПа. Размеры 
КВ: и = 37,5 мм, й = 27 мм, Г = 254 мм, С = 60 мм. 

Упрочнение КВ проводилось по базовой технологической схеме (рис. 1 а-6) на виброустановке 
УВГ 2х50 с объемом рабочей камеры 50 дм”. Рабочая среда — смесь закаленных полированных шариков 
диаметром 3—6 мм из стали ШХ15 твердостью 60 — 62 НКС. Параметры обработки: амплитуда вибраций — 
3 мм; частота вибрации — 25 Гц; время обработки — 20 мин. 

Коробление валов определялось путем замеров радиального биения коренных шеек до и после 
упрочнения индикатором с ценой деления 0,01 мм. КВ при этом устанавливался коренными шейками на 
призматические опоры и поворачивался от руки. На рис. 7 приведены экспериментальные и расчётные 
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значения радиального биения КВ после ВиУО в сечениях конца вала, удалённых от крайней щеки крайнего 
кривошипа на расстояния: А—А = 123 мм, Б-Б = 77 мм, В-В = 17 мм. Из этих данных следует, что отклонение 
фактических величин коробления КВ от расчётных составляет в среднем не более 15—20 %, что позволяет 
рекомендовать зависимость (1) для расчёта коробления КВ при виброударной обработке. 
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Рис. 7. Зависимость радиального биения КВ от местоположения контролируемого сечения: сплошная 
линия — экспериментальные значения, пунктирная — расчётные 


Обсуждение и заключения. Виброударная обработка обеспечивает улучшение геометрических и 
повышение физико-механических параметров поверхностей шатунных и коренных шеек КВ. В результате 
обработки всех поверхностей коленчатого вала этим методом величина коробления не превышает допустимых 
значений, установленных техническими требованиями. Это позволяет сделать вывод об эффективности 
рассмотренного способа упрочнения КВ с целю повышения их эксплуатационных свойств. 
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Введение. Представлены результаты исследований процесса единичного взаимодействия при центробежно- 
ротационной обработке деталей в абразивной дискретной среде. Одновременно с численным моделированием 
проведены эксперименты на центробежно-ротационной установке и исследованы максимальная глубина 
внедрения в поверхность детали, размеров единичных следов, съёма металла за один удар абразивной гранулы. 
Исследован съём металла с обрабатываемых деталей в зависимости от режимов обработки, характеристик 
абразивной частицы и обрабатываемого материала. 

Материалы и методы. Учтены зависимости для определения съёма металла с обрабатываемых деталей 
(сталей 45, меди МОБ и алюминиевого сплава Д16Т) в зависимости от зернистости (№) абразивных частиц, 
рассматривается процесс единичного взаимодействия абразивной частицы с поверхностью заготовки в рамках 
динамической контактной задачи теории упругости. Авторами проведено конечно-элементное моделирование 
рассматриваемых конструкций в САЕ пакете АМЗУЪ. 

Результаты исследования. Представлены результаты теоретических и экспериментальных исследований 
процесса съёма металла с обрабатываемых деталей в зависимости от зернистости абразивных частиц. Описана 
методика их проведения, применяемый инструмент и оснастка. Сопоставлены результаты теоретических и 
экспериментальных исследований. Установлена их высокая сходимость. Подобраны абразивный инструмент и 
режимы обработки. 

Обсуждение и заключения. Построенные в работе зависимости позволяют определить рациональные значения 
технологических параметров процесса ЦРО, и могут быть использованы при проектировании технологических 
процессов ЦРО. Следовательно, для достижения желаемого качества поверхности можно сэкономить 
временные и финансовые ресурсы. 


Ключевые слова: съём металла с поверхности детали, глубина внедрения, единичное взаимодействие, расчёт 
удаления металла, центробежно-ротационная обработка. 
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Введение. В машиностроении технология шлифования всегда обеспечивает высокую точность 
поверхности и является последним шагом обработки поверхности деталей. Материалы, обладающие 
повышенной термостойкостью твердостью и высокой прочностью, могут обрабатываться с использованием 
технологии шлифования [1]. Для повышения эффективности процесса центробежно-ротационной обработки 
(ЦРО) необходимо оптимизировать модель фрикционного взаимодействия между абразивными частицами и 
поверхностью заготовки. 

Обработка в среде свободного абразива позволяет обрабатывать детали различной формы, размеров и 
материалов с использованием простого и надежного оборудования. Температура в зоне обработки значительно 
ниже, по сравнению с процессом шлифования, а обработка сопровождается подачей технологической 
жидкости. В [2,3] анализ термической текстуры и прогноз износа заготовки были выполнены методом 
конечных элементов и при помощи гауссовского процесса. Исследования, представленные в работах [4, 5], 
показали, что при скольжении частицы оксида алюминия происходит разрушение поверхности заготовки и 
съём металла с детали. В данной статье рассматривается процесс единичного взаимодействия заготовки с 
абразивной средой при ЦРО. Одновременно с численным моделированием проведены эксперименты 
исследования глубины внедрения абразивной гранулы в поверхность детали. Исследовано влияние зернистости 
абразивной частицы на съем металла с обрабатываемых деталей. 

Исследования по определению влияния технологических параметров обработки ЦРО на качество и 
производительность изучено недостаточно глубоко. Это сдерживает широкое внедрение в производство 
технологии ЦРО. Для решения этой задачи необходимо получение теоретической модели процесса ЦРО, 
позволяющей прогнозировать результаты обработки еще на стадии проектирования. 

Описание установки. Принцип работы ЦРО (рис. 1) состоит в следующем: заготовка 4 и абразивные 
частицы загружаются в рабочую камеру 1; абразивные частицы и заготовка движутся по спиральной орбите; 
вращающееся днище (ротор) 2 соединено с двигателем; внутренняя поверхность днища покрыта износостойким 
материалом для уменьшения трения. 





технопогическая № 
жидкость 





т зы раны кодек сочые и ==, 


Рис. 1. Схема центробежно-ротационной обработки: 
| — цилиндрическая обечайка; 2 — ротор; 3 — абразивные частицы; 4 — заготовки 
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На рис. 2 изображена трехмерная модель взаимодействия абразивной частицы с заготовкой. 

В работе [6] изучен процесс единичного взаимодействия абразивной гранулы с деталью. Особенности 
этого процесса представлены в работах [7-10]. Из работы [11] известно, что температура в зоне обработки 
невысока и не приводит к изменению структуры поверхностного слоя заготовки. В данной работе исследуется 
контактное взаимодействие абразивной частицы с поверхностью заготовки, абразив движется со скоростью у, 
угол контакта с поверхностью обрабатываемой детали а = 15°—25° (рис. 2). 

Исследование процесса единичного взаимодействия. Решение задачи теоретического 
моделирования единичного взаимодействия позволит исследовать влияние технологических параметров на 
ЦРО. Для создания математической модели съема металла с поверхности детали необходимо описать влияние 
факторов на форму и размеры следов взаимодействия гранулы с поверхностью. 





ЪдЪЫ 


Рис. 2. Модель абразивной частицы и заготовки: 
| — абразивная частица; 2 — заготовка 


Множество работ было посвящено процессу единичного взаимодействия [12, 14, 15]. Анализ работ [12— 
14] показал, что взаимодействие абразивной частицы с поверхностью заготовки происходит следующим 
образом: при столкновении движущейся частицы с поверхностью заготовки на некоторый угол а на 
абразивную частицу действует сила сопротивления Р, которая состоит из касательной Р’‚ и нормали Р». Задача 
сводится к движению кругового конуса с постоянной скоростью у > 0. В работах [5,6] определена 
максимальная глубина внедрения при единичном взаимодействии. 


их = 27 р К Ш 9 (1) 





где р„ — плотность материала; К, — коэффициент, учитывающий влияние шероховатости поверхности 
заготовки на площадь фактического контакта; с — коэффициент несущей способности контактной 
поверхности; о, — предел текучести материала детали; А — радиус абразивной частицы; © — угол встречи 
абразивной частицы с поверхностью заготовки; Г», — эффективная скорость движущейся абразивной частицы, 
определяемая по формуле: 


и 


о (2) 


где К, — эффективный радиус ротора; « — угловая скорость вращения ротора; Ку — обобщенный 


коэффициент эффективной скорости. 
Исходя из работы [4], при столкновении абразивной частицы с заготовкой с образованием 
эллиптического следа обработки значения полуосей а, Б находятся по формуле: 


= \/Е° -(В-й). 
п 
а= 5 (2189 — ЛА, +6, (3) 


где / — коэффициент трения частицы по поверхности заготовки. 
Разрушение поверхностного слоя происходит за счет микрорезания. Поэтому при расчете съема 
металла достаточно учитывать количество взаимодействий Л№>, приводящих к микрорезанию [16]: 
О=М,4, 
где О — съём металла с поверхности заготовки. 
Подставив значение /,, из зависимостей выше, получим: 





2 ет 


5 
О=РВощ р при 5, > 4К” (4) 


О’=Р'Рош при 5, <4В” (5) 


ет 
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Конечно-элементные модели абразивных частиц и деталей были построены в программе АМЗУЗ для 
того, чтобы исследовать максимальную глубину внедрения при единичном взаимодействии на основе 


параметров технологии ЦРО. 
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Рис. 3. Максимальная глубина внедрения при единичном взаимодействии 


На рис. 3 представлено распределение вертикального смещения при ударе и скольжении фрагмента 
абразивной частицы по поверхности заготовки, рассчитанное в АМЗУЪ5, с учетом пластичности материала 
детали (билинейная модель). 

Методика экспериментальных исследований. Исследование процесса съема металла с поверхности 
заготовки проходили на цилиндрических образцах из различного материала: сталь 45 алюминиевый сплав 
Д16Т, медь МОБ, 12 штук на образец (рис. 4). Твердость НВ измеряли на твердомере по Бринеллю. Твердость 
НВ и предел текучести материалов образцов показаны в таблице 1. 

В качестве абразивных сред были выбраны: призмы белые ПТ 15х15 условно приравнивались к 
зернистости 25 (рис. 5 а); фарфоровые шары диаметром 10 мм условно приравнивались к зернистости М60 
(рис. 5 6); ПТ 25х25 условно приравнивались к зернистости 12 (рис. 5 г); конус абразивный бело-зеленый 
зернистость 8, 4=30 мм, В=30 мм (рис. 509). 

Абразивные гранулы и заготовки загружаются в рабочую камеру ЦРО. Обработка выполняется с 
частотой вращения ©=12 об/с. Каждые 30 минут обработку останавливают, образцы вынимают из камеры, 
тщательно промывают и сушат. Для определения съема использовались весы АО 200. Во избежание коррозии 
использовался 0,2 % раствор кальцинированной соды. 





Рис. 4. Образцы для определения съёма металла с поверхности заготовки. 


Материалы: а — сталь ХВГ, 6 — алюминиевый сплав Д16Т, г — медь МОБ 
Таблица 1 


Размеры и механические свойства материалов образцов 


ь 6 П т? 
Материалы образцов о Твердость, НВ ВЕРЕС 
а МПа 


Алюминиевый сплав Д16Т < 30х50 80—82 240 
Медь МОБ © 15х20 56—59 180 
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Рис. 5. Абразивные частицы использовались в ЦРО: а — призмы белые ПТ, 
6 — фарфоровые шары, г — ПТ 25х25, д — конус абразивный бело-зеленый 


Результаты экспериментальных исследований. Проведены исследования зависимости процесса 
съёма металла от зернистости абразива (№). Для комплексной проверки теоретической модели сравнение 
производилось по результатам обработки 12 образцов каждой марки. Результаты теоретических расчетов с 
использованием зависимостей (4), (5) сравнивались с результатами экспериментальных исследований. По 
результатам экспериментальных и теоретических исследований построены графики. 

На рис. 6-8 представлены зависимости съема металла с поверхности детали от зернистости абразивных 
частиц М60; 8; 12; 25, материал детали — сталь 45; алюминиевый сплав Д16Т, медь МОБ, режим обработки — 
число оборотов ротора 2=1206/с. 
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Рис. 6. Зависимость съёма металла от №, материал детали — сталь 45 
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Рис. 7. Зависимость съёма металла от №з, материал детали — Д16Т 
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Рис. 8. Зависимость съёма металла от №з, материал детали — МОБ 
2:1 
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Рис. 9. Сравнение экспериментальных и теоретических зависимостей съёма металла от №: | — медь МОБ; 
2 — алюминиевый сплав Д16Т, 3 — сталь 45 
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Обсуждение и заключения. По результатам теоретических и экспериментальных исследований 
найдены зависимости, представленные на рис. 6-9. Проанализировав результаты проведенных исследований, 
можно сделать следующие выводы: 

1. Съём металла О пропорционален увеличению зернистости абразивных частиц. При обработке с 
©=12 об/с для стали 45 при №=М60 О=0,169 гр, при №=25 О=0,37 гр; для алюминиевого сплава Д16Т при 
№ з=М60 О=0,267 гр, при №з=25 О=0,797 гр; для меди МОБ при №з=М60 О=1,065 гр, при №з=25 О=1,789 гр. 

2. Влияние зернистости абразивной среды, механических свойств материала и режима обработки 
правильно отражает теоретическую модель единственного взаимодействия. 

3. При сравнении результатов экспериментальных исследований с теоретическими данными разница 
составляет менее 20 % (с учетом съёма металла от зернистости абразивных частиц М№з). 

4. Внедрение результатов исследований в производство позволило повысить производительность 
обработки в 1,5—2 раза при требуемом качестве обработанной поверхности и снижении износа рабочих сред. 

Полученные результаты можно использовать для повышения эффективности технологического 
проектирования ЦРО и расширения технологических возможностей. 
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Подход к прогнозированию ущерба из-за неблагоприятного 
стечения обстоятельств, связанного с неразличимостью 
исходных данных 
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›? АО «ВНИИ “Градиент”» (г. Ростов-на-Дону, Российская Федерация) 
°ФГБОУ ВО «Донской государственный технический университет» (г. Ростов-на-Дону, Российская Федерация) 


Введение. При администрировании сложных многопараметрических систем управленческие решения часто 
принимаются в условиях неопределенности. Остро стоит проблема снижения вероятности нежелательных со- 
бытий и уменьшения возможного ущерба. Эффективность прогнозирования ущерба для сложных систем 
напрямую зависит от качества методов обработки, систематизации и количества входных данных. Необходимо 
совершенствовать методы оценки и прогнозирования ущербов и разрабатывать новые подходы и критерии ста- 
тистического прогнозирования ущерба и оценки надежности системы. Решение таких задач осложняется боль- 
шим числом показателей, неопределенностью данных, короткими рядами наблюдений, неполнотой исходной 
информации, недостаточно развитым научно-методическим аппаратом. Существующие методы прогнозирова- 
ния ущерба в системах потенциально опасных объектов не учитывают причины происшествий, случившихся 
из-за неблагоприятных стечений обстоятельств. Как следствие, управленческие решения принимаются на осно- 
вании недостоверных результатов прогнозирования. В связи с этим актуальной научной задачей представляется 
разработка методов и методик для формирования целесообразных управленческих решений, свободных от ука- 
занного недостатка. 

Основная цель работы — рассмотрение частной задачи для прогнозирования ущерба из-за неблагоприятного 
стечения обстоятельств, связанного с неразличимостью исходных данных. Задачи: рассмотреть такой вид не- 
определенности, который включает в себя неразличимости истинного состояния системы и реального значения 
ее количественной характеристики; сформулировать задачу комбинаторики для случая, когда составной весьма 
опасный признак определяется совместным проявлением двух и более простых признаков. 

Материалы и методы. В условиях множественной неразличимости в качестве исходных данных использованы: 
множество неразличимых исходов с достоверной информацией о факте реализации события и неопределенно- 
стью отнесения события к некоторому типу; семейство множеств, имеющих одинаковое число элементов. 
Учтены Декартово произведение семейств соответствующих множеств и фактическое значение группы состав- 
ного потенциально опасного фактора с составным весьма опасным признаком. Представлена результирующая 
моноэлементная группа неразличимости, которая является также возможным событием, полученным в резуль- 
тате пересечения двух необходимых событий. 

Результаты исследования. Установлено, что задача прогнозирования ущерба из-за неблагоприятного стечения 
обстоятельств соответствует задаче комбинаторного типа, состоящей в перечислении всех множеств- 
аргументов. Полученный диапазон, представляющий собой элементную группу неразличимости, характеризует 
меньшее и большее возможное значение численности группы потенциально опасного фактора с составным 
весьма опасным признаком. Показано, что сформулированные комбинаторные задачи без существенных изме- 
нений применимы к задачам в обобщенном виде, когда составные весьма опасные признаки определяются с 
применением не только операции пересечения, но также объединения и разности, а исходными не обязательно 
будут группы объектов с простыми признаками. 

Обсуждение и заключения. Полученные результаты могут быть использованы при построении сложных систем 
мониторинга, моделирования и прогнозирования многопараметрических объектов и динамических систем. 


Ключевые слова: неразличимость, вероятность, математическая модель, риск, случайное событие, вероятность 
происшествия. 


Для цитирования: Золотухин, В. Ф. Подход к прогнозированию ущерба из-за неблагоприятного стечения 
обстоятельств, связанного с неразличимостью исходных данных / В. Ф. Золотухин, А. В. Матершев, 


Информатика, вычислительная техника и управление 


> 
=> 
л 


ЬИр://уезш-допза.га 


> 


А4уапсеа Епятеетиз Кезеагсй 2020. Т. 20, № 4. С. 405—413. 155М 2687-1653 





Л. А. Подколзина // Адуапсеа Епошеегие Везеагсв. — 2020. — Т. 20, № 4. — С. 405—413. В рз://Ао1.0ге/10.23947/ 
2687-1653-2020-20-4-405-413 





Финансирование: работа выполнена при поддержке гранта РФФИ № 19-01-00357. 


© Золотухин В. Ф., Матершев А. В., Подколзина Л. А., 2020 


Ап арргоасй Фо Гогеса$Ипо датасе дие 0 ип?ауогаШе сгситубапсе$ а5$остжей уиП 
п915Япои$ пару о? 5оигсе даа 


у. Е. Хо1ошКыш', А. У. Маегзвеу”, Г. А. РодКои та” 
5? $С “УМП ‘Отафепе” (Возюу-оп-Ооп, Виззап Еедеганчоп) 
оп Зе Тесфиса! ОшхуегзИу (Воз®юу-оп-Ооп, Визап Еедеганоп) 


[итодиспоп. У\УПеп ад “ети сошрех ши!@-рагатаег зуз{еп15, тапасетеп( 4ес1510п$ аге ойеп таде ипдег ип- 
сецатифу. ТБеге 1$ ап асще рго ет оГ тедасйоп оЁ Фе ПКеПвоо4 о! ипуаще4 еуеп5 ап4 пуйсайоп оЁ ро551е дат- 
асе. ТВе еЁРслепсу оЁ ргед1сипе датасе ю сошр!ех зубетз 4ереп4а$ диесЧу оп Фе дааШу оф ргосеззше тео4$, 
зузетайнтаНноп, апа 1е атоций оЁ шриЁ даа. П 1$ гедае4 © пиргоуе те о4$ Юг аззеззше ап4 ргед1сЯп> Чатасе 
ап ю деуеюр пе\и арргоасВез ап4 сгцема Юг зайз$Иса| Гогесазипе ог датазе ап4 еуамайие Фе зуфет тепа иу. 
Тре зоной ю зисВ ргоеилз 15 сошрПсае4 Бу а 1агое патЬег ог 1 1саюгз, даа ипсецалбу, зВой земез оЁ обзегуа- 
Ноп$, шсотр|ее ша тпогтайоп, шзаИсепт Ну деуеюоре4 зслепиЙс тефо4о|о21са] аррагайа$. Ех15Япе тео4$ Юг 
ргефсипе датасе ш Фе зу$етз оЁГ р\епнаПу дапзегом$ об]ес{5 4о по! {аКе шю ассоцп{ Ще саизез оЁ асс14еп{$ фа 
Варрепе4 4дае ю итауогаШе сиситзапсез. Аз а сопзеадчепсе, тапасетепе 4ес15$101$ аге таде ироп ипгейаЫе Гоге- 
сазние гези$. ш {1$ гегаг4, ап игоеп зслепийс 1а$К 1$ е 4еуеортепЕ оЁГ тео4$ апа {есби1ачез Юг Фе Ююгтайоп 
оЁ у1а е тапагетеп{ 4ес1510п$, Нее Нот 1$ зВойсоптше. Те пла]ог зу оБ]есйуе 1$ №0 соп$14ег а рагасшаг ргоБ- 
[ет Юг ргед1сип® датасе дие ю ип ауогае стсит${апсез аззос1ае4 уп Фе ша1$Ипеи1$Ва у оЁ фе ша] даа. 
ТБе {азК$ аге © сопз14ег 11$ Кш4 оЁ ипсецаииу ус шеадез ш1зниои5Ваб у оЕ Фе {чае зубет сопд!1оп ап 
Фе геа| уае о? {5 ацапшайуе свагасет1$8с$; фо Юогти|аже а сот шаюпа| рго ет Юг Ще сазе уПеп а гаег дапсет- 
ои5 сотрозЦе Геахге 15 Чаегилшеа Бу Фе ]ошё тапПебаноп оЁ мо ог тоге зпаре Геагез. 

Маепа! апа Метод. Опдег Ще сопд11оп$ оР шаре ш1зНиоизВаб у, Фе оПомлие \маз изе4 аз Ще зоигсе Аа- 
{а: а зе! ог ша15Ниеил$Ва е очсотез у тепа ]е шогтаноп оп Фе еуе{ шуапсе ап ®е ипсецалибу о? а5$12п112 
Фе еуеп: 10 а сецат фуре; а шПу оР 5е{$ Вауше Фе зате питфег о? @аетеп5. Те Сацезлап ргодисй ог Фе атшез 
ог Ше соггезроп4тсе 565 ап Ше асаа| уа[ае о Фе огоир оЁГ а сотроипа роепнаПу дапоегоч$ Гасюг \и а сот- 
рочпа гаег дапоегоч$ Геабиге аге фаКеп шю ассоии. ТВе гез пе топо-@етепе Ёа77у отоцр 1$ ргезеще4, \мсь 1$ 
а[50 а роз51е еуепё гези ие Нот фе пиегзесноп оЁ мо песеззагу еуеп5. 

Кези!15. [ 15 еба Вед а Ве рго ет оЁ рге сих датазе ие ю итауогае стсатзапсез соггезроп4$ №0 а сот- 
Ыпаопа1-уре рго ет, с соп3155 ш епатеганпе аП зе оЁ агоитеп!. ТВе геза о гапое, \Н1сьЬ 15 ап е|етета1 
отопр оЁ ша15ипеи1зва Шу, свагасеттез фе зтаег ап ]агоег роз51Ые уа[аез оЁ Ве 517е оЁ Ше этоир оГа рофепнаПу 
Чапоегоиз$ Тасог УИ а сотроз(е гафег дапоегои$ Геабиге. Ц 15 зво\’а фа Фе Ююгии|=же4 сотбшаюпа| ргоет$ 
оц $1от1Нсап сВапоез$ аге аррИса е №0 ргоепл$ ш а эепега1те4 юг, увеп сотроз(е гаег дапхегом$ Ееаагез 
аге Чаегилше4 изшх по{ оу Ше орегайоп оЁ пуегзесйоп, БиЁ а[50 иш@пс апа ЧШегепсе; ШегеБу, Фе шйа| огопрз 
аге по{ песезза у ве об]ес{$ у зпаре Геагез. 

Обси5яоп апа Сопсибтопб. ТВе гези 5 оМатед аге Ююсизе4 оп Фе сопзгасноп оЁ апа[уйса| а1еоги таз Гог езбаБ$|- 
шо ша1знпоизраб у ипаег фе топцогше, то4дейпе, югесазипо $ае-геаже4 ргосез5ез ап сошр!ех упатс пи - 
рагате{ег об]ес5. 


Кеунога$: тша15ипои1ва Шу, ргобаб ку, та фетайса| под], г13К, гап4от еуепф, асс14еп{ роепна1. 


Еог спаноп: У. Е. ГоофКыю, А. У. МмегзВеу, Г. А. Ро4Ко]7ща. Ап арргоасВ фо югесазИпе датасе 4ие © 
ипауогае стситзапсез аззослае у шазипешизва у оР зоигсе ааа. Адуапсеа Епотеегте КезеагсВ, 2020, 


уо1. 20, по. 4, р. 405—413. В рз://401.0г2/10.23947/2687-1653-2020-20-4-405-413 
Еипат> трюгтапоп: Ще гезеагсй 1$ зирроце4 Бу КЕЕ! огап{ по. 19-01-00357. 
Введение. Деятельность в техногенной сфере предполагает риски и критические ситуации. Это могут 


быть разрушения систем (техника), потеря управления (военное дело), банкротство (экономика). Проблема про- 
гнозирования актуальна для многих отраслей и напрямую связана с необходимостью совершенствования, раз- 
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вития, разработки и применения математического аппарата средств управления сложными многопараметриче- 
скими системами. 

Остаются актуальными вопросы одновременной обработки динамических массивов разной степени 
структурированности. Прообразами математических моделей, содержащих подобные структуры, являются 
также оптимизационные задачи практического распределения ресурсов в условиях возможных трудно форма- 
лизуемых воздействий [1]. Высокая степень неопределенности процессов снижает целесообразность использо- 
вания ресурсоемких алгоритмов распределения. В то же время необходимо получать множества альтернатив- 
ных решений. Это особенно важно в ситуациях нечеткости и противоречивости информации об имеющихся 
угрозах [2]. Если объект или процесс непредвиденно изменяется, то в момент принятия решения может отсут- 
ствовать адекватная математическая модель. При этом отказ сложных производственных объектов нередко ста- 
новится причиной техногенных чрезвычайных ситуаций с серьезными экономическими, экологическими и со- 
циальными последствиями, что обусловливает необходимость совершенствования математических основ ана- 
лиза рисков [3, 4]. 

Стратегические решения по управлению сложными многопараметрическими системами принимаются в 
условиях неопределенности. Цель управления современными рисками — недопущение критической ситуации. 
В случае же реализации негативного сценария необходимо минимизировать потери. Для прогнозирования 
ущерба и принятия адекватных решений ведется постоянный мониторинг, выявляющий основные факторы реа- 
лизации критической ситуации |5]. 

Необходимо ранжировать и систематизировать риски по степени влияния на деятельность защищаемо- 
го объекта. Любая система характеризуется множеством параметров, которые зачастую носят случайный харак- 
тер и подвержены изменениям. Эффективность прогнозирования ущерба для сложных систем напрямую зави- 
сит от качества методов обработки, систематизации и количества входных данных. Поэтому для оценки време- 
ни и условия наступления критической ситуации применяются комбинированные методы прогнозирования, 
содержащие экспертную, аналитическую, имитационную части и использующие аппарат теории вероятно- 
сти [6, 7]. Критические значения можно определить аналитически по результатам обработки эксперименталь- 
ных данных. Такая работа зачастую выполняется в условиях недостатка информации, поэтому возможны 
ошибки. 

Количественной характеристикой степени опасности (или безопасности) систем и ситуаций является 
риск. Риском здесь и далее будем считать ущерб (последствия) от реализации возможного события в некоторой 
сложной системе [8]. Ввиду сложности и разнородности факторов, которые влияют на систему и должны быть 
учтены в процессе принятия решений, необходимо использовать совокупность методов анализа и обработки 
информации для оценки ущерба, полученного вследствие нарушений работоспособности систем. 

Зачастую из-за сложности и дороговизны контроля часть параметров, связанных с техническим состоя- 
нием системы (далее — ТСС), является неразличимой. В общих случаях неразличимость рассматривается как 
неопределенность состояния, которая может проявляться не только в будущем, но и в настоящем (или ближай- 
шем прошедшем) времени [9-10]. Неразличимость описывается в терминах теории возможностей, в которой 
разработаны методы получения простых числовых характеристик — оценок различных ТСС. Также это соот- 
ветствует принципу принятия решений, основанных на фактах. Согласно этому принципу, неразличимость со- 
стояний системы вызвана неопределенностью знаний о ней и представляет конкретизацию неопределенности 
знаний о функционировании, физическом поведении. Неопределенности могут иметь экзогенный и эндогенный 
характер, обусловленный соответственно внешними и внутренними воздействиями [11]. Каждый из видов не- 
определенности может существенно ухудшать точность принимаемого решения. 

Современное состояние проблемы. Значительное количество работ посвящено отдельным аспектам 
проблемы прогнозирования и оценки рисков при неразличимости данных. Однако эти изыскания лишены об- 
щего методологического базиса. В исследованиях [12—27] предлагаются частные решения, но общий математи- 
ческий аппарат все еще не разработан. Так, в [14—15] принимается, что оценка каждого из неразличимых между 
собой критериев равна среднему арифметическому их номеров. За ранг каждого из неразличимых критериев 
берется номер всей группы как целого объекта в упорядочении. 

В статье [16] предложено унифицировать подходы к управлению комплексной безопасностью различ- 
ных систем, а также показана ситуация, когда эксперт не различает некоторые критерии. 

В [17] рассматривается задача построения интервальных оценок для неизвестной вероятности при 
наличии в результатах опытов множественных неразличимых исходов. Предлагаются два способа решения: 
учитывать все неразличимые исходы или же отбрасывать их. В обоих случаях возможны огрубление результата 
и ошибки. 
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В работе [18] уточняются неопределенность и неразличимость, возникающие при диагностике состоя- 
ний энергоустановок. Неразличимость понимается как неопределенность состояния управляемого объекта для 
наблюдателя. При этом принимается, что неопределенность снижается по мере движения вниз по иерархиче- 
ским рангам управления энергетическими системами. Кроме того, предложено ввести некоторый порог и счи- 
тать решения неразличимыми в случае, если квадрат разности искомого значения не превышает установленно- 
го порога. 

Аспект неразличимости упоминается в свете развития теории приближенных множеств в [19]. Раскры- 
ваются понятия нижней, верхней аппроксимации и пограничного региона. Это позволяет создавать решающие 
правила «если ..., то» и оперировать только фактами, без допущений. Ссылаясь на [20], автор уточняет, что на 
основе теории приближенных множеств разработан один из главных методов принятия решений в области мно- 
гокритериальной оптимизации. 

В статье [21], посвященной безопасной эксплуатации грузоподъемных кранов, предложен метод экс- 
пертной оценки частоты неблагоприятного события, позволяющий разработать рекомендации по уменьшению 
риска. 

В [22] рассматриваются проблемы военно-оперативных исследований в следующем контексте: военные 
корабли стремятся сдерживать и пресекать морской разбой, а модель движения пиратов основывается на про- 
гнозе вероятности пиратства и на допущении Маркова. Автор использует решение задачи потока с минималь- 
ными затратами. Количество поисковиков не имеет значения. При этом предполагается, что они идентичны и 
неразличимы. 

В исследовании [23] изучены операторы Т-неразличимости с геометрической точки зрения как частный 
случай обобщенных метрических пространств для дальнейшего применения при изучении нечетких подгрупп. 

В [24] показано, что существующая неразличимость данных из множества наблюдений не позволяет 
давать точные оценки состояния системы. Значит, точный прогноз должен основываться на плотности вероят- 
ности неразличимых состояний. Эту плотность можно рассчитать, сначала вычислив оценку максимального 
правдоподобия состояния, а затем ансамблевую оценку плотности состояний, которые неотличимы от состоя- 
ния максимального правдоподобия. 

В [25] представлена характеристика функций, которые позволяют объединить операторы частичной Т- 
неразличимости (отношений) в новый набор. В [26] рассматривается агрегация операторов частичной Т- 
неразличимости и частичных псевдометрик. Анализируется агрегация набора операторов частичной Т- 
различимости и демонстрируется связь между функциями, которые: 

— объединяют операторы частичной Т-неразличимости, 
— сохраняют частичные Т-псевдометрики в процессе агрегации. 

В работе [27] обосновывается связь неразличимости с нечеткими подмножествами. Доказывается, что 
базовое для них отношение является решеточным изоморфизмом. 

Итак, известные подходы к проблеме неразличимости воздействующих факторов функционально огра- 
ничены. Они не обеспечивают необходимый уровень достоверности при принятии управленческих решений. 
Очевидны слабые места их научно-методической базы. Это в целом сдерживает развитие систем прогнозирова- 
ния с учетом фактора неразличимости и доказывает востребованность разработки и совершенствования мате- 
матического аппарата. 

Таким образом, важная задача управления рисками — прогнозирование ущерба, вызванного неблаго- 
приятным стечением обстоятельств ввиду неразличимости исходных данных. Ее решение представляет значи- 
тельный теоретический и практический интерес для многих сложнопараметрических динамических систем. 

Существующие методы прогнозирования ущерба в системах потенциально опасных объектов не учи- 
тывают причины происшествий, случившихся из-за неблагоприятных стечений обстоятельств. Как следствие, 
управленческие решения принимаются на основании недостоверных результатов прогнозирования. В связи с 
этим представляется актуальной научной задачей разработка методов и методик для формирования целесооб- 
разных управленческих решений, лишенных указанного недостатка. 

Постановка задачи исследования. Ущерб, вызванный влиянием опасных факторов и реализацией 
угроз, не всегда удается заранее спрогнозировать из-за сложности различения обычных состояний многопара- 
метрических объектов и систем, поведение которых не детерминировано. В результате возникает необходи- 
мость исследования и оценивания двух независимых видов неопределенности: неразличимости и недетермини- 
рованности. Рассмотрим модель прогнозирования из-за неблагоприятного стечения обстоятельств, связанную с 
неразличимостью исходных данных. 

Покажем, что в терминах теоретико-множественного моделирования математическая модель неблаго- 
приятного стечения обстоятельств представляет задачу комбинаторики с общими исходными данными. Сфор- 


Золотухин В. Ф. и др. Подход к прогнозированию ущерба из-за неблагоприятного стечения обстоятельств 





мулируем эту задачу для случая, когда составной весьма опасный признак определяется совместным проявле- 
нием двух и более простых признаков. 

Материалы и методы. Исходные данные: 

— множество неразличимых исходов И с числом элементов и, содержащее достоверную информацию о 
факте реализации события и неопределенность отнесения события к некоторому типу; 


— семейство А, множеств типа А;, имеющих одинаковое число и, элементов. 
В условиях множественной неразличимости мы вынуждены для каждого 1 рассматривать все множе- 


ства типа А;, потому что любое из них может представлять группу объектов, которые приобрели 1-й простой 
признак. 


Рассмотрим Декартово произведение семейств Ах, х...х А,. Его элементами являются все наборы ( 


А,4, ..., А.) множеств типа 4,4, ..., А. 


Рассмотрим ситуацию, когда множество И = 1, 2,3,4}, К =2 множества типа А содержит по и =3 


элементов, а типа А, — по и, = 2 элементов. Допустим, 4 г = $1, р. А, 5 = $1, 4}. Тогда 


— фактическое значение группы составного потенциально опасного фактора с составным весьма опас- 


ным признаком. В этом случае семейства А И А, имеют вид: 
А = {{1,2,3}, {12,4}, {13,4}, {2,3,4}}, (1) 


А, ={11,2},{1,3}, {1,4}. {2,3}, {2,4), {3,4}. (2) 


Состав групп потенциально опасного фактора с простыми признаками имеет вид (3, 4), представляю- 
щий мономножественные группы неразличимости: 


Ав: 12,3, 4}, (3) 


А, 


в:{2,3,4}- (4) 


Как видно, полученные семейства (1) и (2) содержат по С; =4,С; = бэлементов, представляют поли- 


элементные группы неразличимости и являются необходимыми событиями. 
Декартово произведение этих семейств будет представлять совокупность, содержащую следующие 


3 3 Я Я 
п. =С, -С. = 24 различные пары элементов А и 4, : 


(1,2,3} {1,2} „({1,2,3} {53} ,..-((1,2,3} {3,4}, 
({1,2,41 {1,2} ),...({1,2,4} {1,4}),...,(1,3,4} {3,4}, (5) 
Е В и ‚({2,3,41 {3,4}. 

Можно увидеть, что А х А, является семейством различных пар множеств типа А, А, ‚ имеющих оди- 
наковое число элементов — соответственно 3 и 2. При А>2 Декартово произведение А х А, а А, будет 
семейством наборов, каждый из которых содержит по одному множеству типа 4, 4,,..., А,. 

Будем называть А-наборами упомянутые наборы — элементы Ах, х...х А, . Каждому А-набору поста- 
вим в соответствие пересечение входящих в него множеств, то есть выделим в таких множествах общие эле- 
менты. Например, пересечениями для шестой и девятой пар в (5), имеющих вид (1, 2,3} {3,4} и (1, 2. 4} {1,4} - 
будут соответственно множества 11 м 3} ®. {3,4} = {3} и $1, 2,4} (Г 11,4] = {1,4} | 


Общий результат (5) будет иметь вид: 





1,2} 0,3} № 23102} 2 

1,2} {1} 154} {2} {2,4} 4 (6) 
034 8 04 849 | 

{2} 13} 4} 12,3} {2,4} 13,4} 
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Выражение (6) является результирующей моноэлементной группой неразличимости, а также возмож- 
ным событием, полученным в результате пересечения двух необходимых событий (1) и (2). 
Теперь рассмотрим функцию Г с аргументами, являющимися упомянутыми пересечениями множеств 


типа А ПА, ^...С А, ‚ и со значениями, равными числу элементов этих пересечений: 
(АА о... СА, ) =АбА п... СА ЕР, (7) 


где К — число элементов соответствующего множества. 


Для рассмотренных выше шестой и девятой пар: 
1({1,2,3} 93,4} = /({3})=1 
#(@,2.4}5,43) = 7 (6.4) 
Теперь можно увидеть, что речь идет о задаче комбинаторного типа, состоящей в перечислении всех 


множеств-аргументов для получения по ним различных значений. Такое перечисление можно найти после 


определения меньшего и большего значений та ‚ которые совпадают с меньшим и большим возможными значе- 


ниями численности группы объектов с составным весьма опасным признаком. 

Результаты исследования. Рассмотрим следующее определение: С — возможное при осуществлении 
Е событие, если Е ПС 0. 

Будем исходить из того, что понятие возможности события обычно связано с предположениями в не- 
определенных обстоятельствах. Поэтому применение этого понятия неуместно, если известны Ё, С и осуще- 
ствившийся исход х ЕЁ. С появлением Е ясно, что С либо имеет место, если хЕ С, либо нет — в противном 
случае. При этом неопределенность отсутствует. 

Самым распространенным, продуктивным и достаточным для наших целей является понятие возмож- 
ного события, обусловленное неопределенностью в виде неразличимости появившегося исхода среди других 
исходов ЕР. В связи с этим о С как о возможном событии говорят при перечисленных ниже дополнительных 
условиях. 

Условие 1. Известны события Ё и С и, следовательно, их пересечение Е ПС + © (по определению, это 
не пустое множество). 

Условие 2. Не установлен исход х;, из-за которого происходит Е, если Е = {х1, хо, ..., хи} принадлежит 
более одного элемента. Известно лишь, что х; — один из элементов А, но не установлено, какой именно, т. е. 
Н!х, ЕЁ. 

С учетом условий 1 и 2 конкретизируется общепринятый смысл утверждения, что С — важное собы- 
тие. Итак, С либо осуществлялось, если неразличимый фактический исход х; Е ЕП С (а значит, х; Е С), либо не 
осуществилось — в противном случае (если х; & ЕП С). Рассмотрим пример, иллюстрирующий возможное со- 
бытие такого рода. 

Пусть Ё и С состоят соответственно в появлении четного числа и больше трех очков при бросании иг- 
ральной кости, то есть Е = {2, 4, 6}, С= 44, 5,6}. Здесь Е ПС = {4,6} = 0. Известно, что Е произошло, но не 
установлено, какой из его исходов имел место. В результате соблюдены выше приведенные условия, согласно 
которым С — возможное событие: оно осуществилось, если имел место либо исход х.› = 4, либо исход хз = 6, 
и не осуществилось, если реализовался исход Хх! = 2. 

Важным частным случаем возможного является необходимое событие С. Оно имеет место, если 
ЕПС = © (условие 1), соблюдается условие 2 и ЕС С, то есть в случае, когда все исходы Е принадлежат С. Ка- 
кой бы из них ни появился, он приводит к осуществлению С. 


В данном примере по п = 4, п; = 3, и> =2 найдем, что { принимает два значения: 1 и 2. Это свидетель- 


ствует о том, что имеются лишь одно- и л двухэлементные пересечения А г. множеств типа 


А, со, И = 41, р. 4} ‚ т.е. в группе потенциально опасных объектов с составным весьма опасным призна- 


ком может быть либо 1, либо 2 объекта, и группа неразличимости (8) является числовой моноэлементной груп- 
пой неразличимости: 


{2}. (8) 
Следующей вытекающей из рассматриваемой модели неблагоприятного стечения обстоятельств явля- 
ется частная задача, которая заключается в том, что каждое полученное выше значение / нужно разделить на 


п, т. е. найти ее нормированные значения — набор альтернативных частот приобретения произвольным объек- 
том составного весьма опасного признака в условиях множественной неразличимости и недетерминированно- 
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сти. Последняя состоит в том, что частные иди представляют частоты появления на произвольном объекте 1-го 
простого признака для всех 1 = 1, ..., К. 
Сформулируем цель решения следующей комбинаторной задачи, соответствующей частотной задаче: 


определить число 5 элементов Ах А, х...х А, , найти 5(т)/5 для всех и. 


Частные 5(7)/5 для всех г представляют собой вероятности того, что число объектов с составным весьма 
опасным признаком будет равно определенному значению г. Для вычисления таких вероятностей в этой связи 
применим зависимость: 


и т. 
5=С" С»... С". (9) 

Аналогично, расширив рамки частной задачи, можно определить представляющий особый интерес 
набор альтернативных частот (отсутствия на произвольном объекте составного весьма опасного признака). Его 


можно найти, заменив в исходных данных множества типа А, А, на их дополнения типа (\А\, (ЛА), с числом 


элементов и — и! = [ии- и> = 2 соответственно. При этом на предварительном этапе аналогично определяются 
семейства, Декартово произведение семейств, элементам которых ставятся в соответствие пересечения мно- 
жеств типа (ЛАП 0\4>), и функция, значениями которой являются числа элементов таких пересечений. Опре- 
делив эти числа, решим комбинаторную задачу, соответствующую одной частной задаче, а после нормализации 

— другой. В терминах моделирования первой задачи применительно к новым исходным данным вторая 
представит: 

— множество И с числом элементов и, 

— различные множества (ЛА! с одним и тем же числом элементов и — ит, 

— различные множества (ЛА› с количеством элементов и — и> каждое. 

Важное значение для практики имеет соответствующая проблеме комбинаторная задача с новыми ис- 
ходными данными. Для ее решения целесообразно использовать результат следующих рассуждений. Пусть 
наряду с набором упорядоченных по возрастанию значений / получен набор аналогично упорядоченных значе- 
ний — такой же, как } функции ©, аргументы которой соответствуют новым данным. Можно показать, что пер- 
вый и второй наборы имеют одинаковое число элементов. Пусть /» и ®, — члены этих наборов с номером й. 
Можно убедиться, что члену ©, второго набора соответствует такое же число аргументов функции ©, как и чле- 
ну /, первого набора, и число всех аргументов (вариантов) Ги © одинаково, то есть равно 5. 

Следовательно, произвольному значению 2, функции © соответствуют в терминах модели частные типа 
5(7)/5, где г = [». Эти частные представляют соответствующие вероятности наличия групп объектов без состав- 
ного весьма опасного признака различной численности, если вероятности появления различных групп объектов 
одинаковой численности с простым признаком 1 для всех 1 = 1...К одни и те же. 

Можно показать, что сформулированные комбинаторные задачи без существенных изменений приме- 
нимы к задачам в обобщенном виде, когда составные весьма опасные признаки определяются с применением 
не только операции пересечения, но также объединения А ‹/А, ‹..../А,, разности А\/А,, А, \А ит. п., 


а исходными не обязательно являются только группы объектов с простыми признаками. 

Обсуждение и заключения. Итак, установлено, что задача прогнозирования ущерба из-за неблагопри- 
ятного стечения обстоятельств соответствует задаче комбинаторного типа. Она состоит в: 

— перечислении всех множеств-аргументов, 


— выводе различных значений }, 
— определении меньшего и большего значений {, 


— получении диапазона возможных значений функции /, в котором находится его фактическое значе- 
ние. 

Этот диапазон представляет собой элементную группу неразличимости, характеризует меньшее и 
большее возможное значение численности группы потенциально опасного фактора с составным весьма опас- 
ным признаком и включает в себя несколько операций. 

Полученные результаты ориентированы на построение аналитических алгоритмов установления нераз- 
личимости в процессе мониторинга, моделирования и прогнозирования процессов, связанных с состоянием, и 
сложных динамических многопараметрических объектов. 
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Сткема оГ еуаша_по апотеще4 геаШу аррНсайоп$ 


М. В. АЫуаеу, А. М. АБПаКитота, 1. 5. зеаДатебоуа 
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[итодисноп. Тфе Не оЁ{ ааотеще4 геаШу (АК) 1$ ого\лис гар1 Чу ап4 Баз огеа{ ауапсез ш ицегасноп, паутеайоп 
ап гасК1п®. Мо\уадауз Шеге аге а [01 оЁ Неп4$ юг АК аррИсанопз ш ЧШегепЕ агеаз (едисаноп, ещецаштитепь, Би$1- 
пез$, тесте, ес.). Но\еуег, Феге 15$ а 1асК оф тезеагсВ фо ргоу14е Фе еуашанпе АВ аррз Нате\могК © зиррой 4е- 
уеорегз \Пеп сгеайпо зийа Ме АВ аррйсайопз$ Юг зрес1Йс пее4$. \е ргоу14е а ргасйса| арргоасВ © ацапиу зоте 
ог Ше АК арр|саноп$ Ееаге$. \Ме Юси$ оп Ше деуеортеп: оЁ сгцепа Гог еуашайпе ааотешще4 геаШу аррИсаНоп$. 
\е 415си5$ фе сгцепа оЁ споозше аппепз1о1$ Гог ай расе зисВ аз $апдага$ Гог АК, 1001$ Юг АК 4еуе|ортетф, пау- 
1сайоп ап4 маске, сомеп тапасетет, изаШиу. \Ме ргоу14е апа[у$15 ап4 еуаайоп ог АК аррз @гоиэВ еасВ сВаг- 
ацегзис изе ои14ейпез$ у мсВ \е Вауе 4еуе]оред. 

Маепта! апа Мефосд$. Ап АК аррйсайоп 1$ а зоЙ\аге аррИсаНоп а{ ищеогаез 1е(а| у15иа1, аа до ап оег {урез 
оР сощепе ш® а геа1-\о4 епугоптепе. ТВе зоЙ\аге диаШу ап4 ре{огтапсе аге {е таш свагасег1$Ис$ оЁ Фе ар- 
рИсайоп, уасВ аге Кеу асюг$ Юг АК аррПсайоп$. ТВе апа[уз1$ оЁ зслепийЯс рарегз, Чоситеп5 ап4 запдаг$ таде и 
роз$1Ые ю д&егите свагас{ет15Нс$ Ша аге Ще 11051 $1е1Исап дааШу ш1саог$ Базе оп \еП-эточи4е4 изег$” пее 45 
ап4 аетап($. 

Кези5. 'ТВе сгцепа \ууе Вауе деуе]оре4 Гог еуаайпе аррПсаНоп$ у аиотешще4 геаШу епаШе 4еуе]орег$ №0 стежме 
ег о\уп зоЙ\аге ргод4ис{5 ш $асез, Базе оп %ер-Бу-%ер гедишетеп{$ ог Фет, еуашанпе фе деуе]ортепй ргосез$5 
Бу сБагацегзИс$. ТЬ1$ арргоасЬ \Ш аПо\ уоц 0 сгеме ыеВ-дааШу зой\аге ргодис{$ изшо запдаг1те, то4еги 
еуе!ортепе 10015. 

Обси5яоп апа Сопсшяоп5. ш аа4\оп, деуеорегз \Ш Вауе а д&аПе4 ипдегапЧ те оЁ еасВ асе оЁ стеайпо Фе 
аррПсайоп ап4 Фе песеззагу деуе]ортеп 10015 апа {есВпо]о21ез ф0 оМаш Ше Ь12Ве${ ацау гези. Тваё м Ш олуе ап 
орробипйу тю есле оп зресИс 4еуортепе 10015, те о4$, ппо4е!5 апа {есбпо]ослез Беюге Уагипе \уогК оп а 
рго]ес+. Аз а гезо №, ® \ Ш ргом14е Ше Нпа| №121 -диаПу зоЁ\аге ргодисё у соо ежепз16Ииу ап сотрйапсе ум 
фе то4еги гедитетеп ое 41а! шдизгу пагКей. 


Кеунога$: ааотеще4 геаШу, зоЁ\аге 4еуе]ортепЕ Ки, пау1сайоп апа {тасКте, сощеп тапазетепь иза у. 
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Ноп5. Адуапсеа Епошеегис Везеагсй, 2020, уо|. 20, по. 4, рр. 414—421. 6@рз://401.0ге/10.23947/2687-1653-2020-20- 
4-414-42] 





© АБуает М. К., АБПаютота А. М№., ве4атеюта /. 5., 2020 


Критерии оценивания приложений дополненной реальности 


М.Р. Абляев, А. Н. Аблякимова, 3. С. Сейдаметова 
ГБОУВО РК «Крымский инженерно-педагогический университет имени Февзи Якубова», (г. Симферополь, Республика Крым) 


Введение. Дополненная реальность (АК) является быстроразвивающейся отраслью, которая значительно про- 
двинулась в вопросах взаимодействия, навигации и отслеживания объектов. В настоящее время существует 
множество направлений АК, которые применяются в различных областях (образование, развлечения, бизнес, 
медицина и т.д.). Однако отсутствуют чёткие критерии оценки фреймворков АК-приложений для поддержки 
разработки таких приложений под конкретные задачи. В статье предложен подход, позволяющий количествен- 
но оценить некоторые характеристики АВ-приложений. В данной статье сделана попытка разработки критериев 
оценки приложений дополненной реальности. Предложены критерии выбора измерений пространства разра- 
ботки, такие как стандарты для АК, инструменты для разработки АК, навигация и отслеживание объектов, 


Абляев М. Р. Критерии оценивания приложений дополненной реальности 





управление контентом, удобство использования (юзабилити). Предоставлены анализ и оценки АК-приложений 
по каждой предложенной характеристике. 

Материалы и методы. АВ-приложение — это программное приложение, которое интегрирует цифровой визу- 
альный, аудио-контент и другие типы контента в реальную среду. Качество и производительность программно- 
го обеспечения являются основными характеристиками приложения, которые являются критически важными 
факторами для приложений АК. 

Анализ научных работ, документов и стандартов позволил определить характеристики, которые являются 
наиболее значимыми качественными показателями, основанными на обоснованных потребностях и требовани- 
ях пользователей. 

Результаты исследования. Разработанные нами критерии оценки приложений с дополненной реальностью поз- 
воляют разработчикам создавать свои собственные программные продукты поэтапно, основываясь на пошаго- 
вых требованиях к ним, оценивая процесс разработки по характеристикам. Такой подход позволит создавать 
высококачественные программные продукты, используя стандартизированные, современные средства разра- 
ботки. 

Обсуждение и заключения. Разработчики будут иметь детальное понимание на каждом этапе создания прило- 
жения, что позволит выбирать необходимые инструменты и технологии разработки для получения результата 
наивысшего качества. Это также даст возможность определиться с конкретными инструментами, методами, 
моделями и технологиями разработки до начала работы над проектом. В результате будет получен конечный 
высококачественный программный продукт с хорошей расширяемостью и соответствием современным требо- 
ваниям рынка цифровой индустрии. 


Ключевые слова: дополненная реальность, комплект разработки программного обеспечения, навигация и от- 
слеживание, управление контентом, юзабилити. 
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шеодисНоп. ТБе деуеортепЕ о? Фе то4еги тагКей о? аисотещеа геаШу (АК) 1есбпо]осу соптЬщез о Фе 
етегоепсе оЁа 1агоег питбег ог АК-аррИсаНоп$ у14еу изе4 ш шдизелез шсТадта® Беа саге, раб с заеу, саз ап 
01|, {ош5т ап тагКейпе, емецалатеп{ ап4 асадепла. ш соппесноп ул Ше шсгеазе4 ищегез{ 1 15$ {еспио]осу, Ве 
деуе]ортепЕ о? ГипсНопа| сара вез оГ АК рго]ес{5 1$ а|$о0 сагле4 оц ууВлеВ ити ез ап шсгеазе ш Фе пее4$ о? 
АК-есрпо]осу из 1ш Ше 11031 Чтуегзе зесфогз оР тодеги зослейу. 

ОпюЮнипа у, зоте $1211Ясап{ азрес{5 оГ Фе АК аррИсанопз$ деуеюортеп ап4 паретещайноп апа АВ $ег- 
У1се5, оцеп 1епоге4 Гог Фе Д4ез1еп зпарИсИу ап пиретепипо зрее4, аге е сотрПапсе узи зисВ зузетз ю геа| 
сопЧ!юоп$ ап фе еуашайние ипдег геа| орегайие сопд1опз$. То деуеор зассез$@ ап меШу еЁРйсет АВ зу$ет5$ 
Ба сап Бе адоре4 ш еуегудау зсепа11о$, азег аззеззтепй ап Гееаск аге уегу пароцаиа [1]. 

Апстеще4 геаШу 1$ а уегу уоци® шдизу, ап4 ете аге зи] по сепегаПу ассереа запдага$ юг деуеорше 
АВ аррИсанопз. АШоиэВ бе таш р!аогиз ог сгеание апотеше геаШу аррИсанопз Вауе 5ееп дейпеа: АВКи', 
АВСоге” апа Опиу” (тозНу изе) — еуеп по\и АВ аррИсаНопз деуеюре4 изше ШФезе ю015 аге ауаЙаЫе ог 1аппсВ 
оту опа птпеа фуре оЁ 41 а| ап тоб!е 4еу1сез. 

Тре зпиаНоп 15 ехасЙу Ше зате \и \еагае 4еусез, Ноо[Гепз ап Маелс Г[.еар 1аз5ез — еасВ тапиасвиг- 
ег оНег$ {5 омут ип1аче зоЙ\аге юг сгеайпе аиотеще4 геаШу ргосгаплз. Ехрей$ олуе зеуега| уеагз Беоге ше шдиз- 
ту сопзо4а{е$ ап соттоп зап4ага$ аПо\ такте е 4еуе!ортепЕ о? аиотеще4 геаШу аррПсайопз ассезз1е ап 
иптуегза| ог а] АВ деуеюорегс. 

Кар14 ргозтез1оп оЁ Ме АВ Не гедигез еНесйуе ап4 уа|й4ае тео4$ о{ 4ез1еп еуашайоп о Бе 4еуе]- 
оре4. ЕаЙиге ю сопз14ег Фе изаб Ииу оЁ АК аррИсанопз$ диние Фе 4ез1еп ргосез$ \уШ тез ш ап шсгеазе ш чзег 
еггог5 ап4 асс14ет, пит те ег 131 ог Пе {1есппо]оэу ап ипдегитите изег регсерйоп$ оЁ Фе {есбпо]осу, Юг Бо 
АК апа Ушба! КеаШу (УК) 1есбпо]оглез [2]. 

Тре аи ог оЁ Ще рарегз [1], [2], [3] ргомае4 ап оуегуле\м оЁ Фе ппропапё 4ез1ет® ап пар!етепиио Геа- 
{игез оГ АВ аррИсайоп$ ап4 ргорозе еогейса| еуашайноп оЁ АК зузетз ап Нате\муогк$ @гоиеВ фе запдага1та- 
оп азрес{5. Епа$[еу апа офегз [4] дезстеа ришстр[ез о? аез1еп Веиизис$ Юг АВ Юг шшШ-аппепз1опа! ааотещеа 
епугоптеп 5. Зоте ехатр/ез оРзег ехрепепсе еуааНоп$ \еге ргезеще4 ш е рарегз [5], [6], [7], [8], [9], [10]. 

Г. Г. СабБаг4 апд .. Е. Э\ап [11] ргорозеа а Озаб ку Епошеегие (ОЕ) юг Аисотещеа КеаШу арргоасВ ай 


' Ацпотешеа ВеаШу — Арр!е Оеуеорег. Вирз://4еуе]орег.арр!е.сот/аиотлегие-геаШу/ 
* АВСоге — Сооб]е Оеуе]орегз. В рз://Чеуеорегз.ооз]е.сот/аг 
°Опйу Ююга|. БИрз://лайу.сото/ 


Информатика, вычислительная техника и управление 


мч 
= 
л 


ЬИр://уезш-допза.га 


416 


А4уапсеа Епатеетиз Кезеагсй 2020. Т. 20, № 4. С. 414—421. 155М 2687-1653 





115е{$ Цегайуе[у а зетез оЁ изег-Базе за ез шю а Над 1юопа1 изабШу-епошеетие Пе сусе. Зеуега! изабиу {е${- 
шо тео4$ оЁ Фе АК аррПсайоп (заБ]есйуе шеазигетеп изшо Битап регсерноп, об]есиуе теазиге Нот обзегуа- 
бой, еуааноп Бу ехрей ФгоиеЙ соотуе уа@гочэв, Веиизис еуашайоп, 1аб обзегуайоп, даезНоппаге) \уеге 
езсгфеа Бу Ргапою апа оегз [12]. Магил$ ап оегз [13] ргезещеа ргасиса[ изе оР Фе иза у шефо4$ Гог еуа- 
айпс ап АК сЫ@агеп’$ БооК \иВ шире тефо4$. Оег азрес$ оГ АВ {есвпо]осу {0015 ап АК аррПсайопз еуа[ча{- 
шо \еге ргезеще4 ш е рарегз [14] апа [15]. 

АКоч ап офегз [16] адар{е4 ап ех15Нпе 4еер 1еагита® агсбиесваге ю зо]уе Фе дщесйоп ргоМет ш АК ар- 
рисайоп изо сатега-базе4 гасКт». Офег тео4$ Ваё Вер №0 ппргоуе АК аррПсайоп ап4 © шсгетепе Ве ргодис- 
ПуЦу ш тапийсвте \уеге аезсгеа ш [17], [18], [19]. 

Треге аге зеуега| зигуеу рарегз оп АВ деуеортепь Би{ попе 1$ дефдсаже4а ю Мое Аиеотещеа ВеаШу. 
Ниапс ап4 оегз [20] ргезепе Ше гези65 ог Ше [а{е31 {есБпо]о21ез ап те о4$ зигуеу аё пиргоуе$ гап-Ите регЮг- 
тапсе апа епегоу еРйслепсу Юг Ще ргасйса| пиретещайоп оЁ тобШе АК арр|Псайопз. 

\е Вауе ргезещеа АК аррйсайоп ““ГИзиий апНег” (“Мазлс [е{ег$”) ап зресла1 Геафигез ог 4ез1еп апа 4е- 
уе|орше шобШе АК аррИсайоп г епВапсше еа]у Щегасу $К1$ ш ог рарегз [21], [22]. 

эоЁууаге фиаШу ап4 ре Тогтапсе. Ап АК аррПсайоп 15 а зоЙ\аге аррПсаНоп а пмеотаез а1епа1 у1зиа1, 
ай41о ап оег {фурез оГ сощеп ш а геа|-\/от14 епутоптет. ТБе зоЙ\аге ацаШу ап ре огтапсе аге фе таш 
свагасет15$Яс$ оЁ Фе аррИсаНоп, у шсВ аге Кеу сюгз Юг АК аррйсаНоп$. зоЙ\аге деуеюортеиа 15$ а сошр[ех ап 
па[а-Рсее4 ргосез$, ш \мсьЬ а 1агое питабег о? зрес1а1${5 оР уапоц$ агеаз оРГехрегизе ап4 уаг1о$ $КШ ]еуе]5 рагис- 
рае. ш аа оп, тапу фесбитса1, {есбпо]021са1, ап тапасема] 1554ез ищете ш Ше аррПсайоп 4еуе]ортеп рго- 
сез5. ТВе зиссез$ оР Пе рго]ес{ ап е ацаШу ое деуеоре4 ргодис* 4ереп4 оп феи адеацае шуо]уетепе. 

Тре апа[у$15 о рарегз, Чоситеп (е.о., [2], [3], [10], [16], [22]) ап апдаг4$ (е.с., ЗО-ТЕС ТГС 1 ЪС 24, 
18521-1', 130 9000) зВо\з фа Фе юПо\ие сВагасег$Нс$ аге Фе Кеу апаШу ш@саюгз Базе оп \уе-огони4ед 
изег5’ пеед$ ап 4етап4$: 

— ша4деацае РапсНопие оЁ Фе зоЙ\аге рго4исЕ; 

— шза стей ищегасйоп оЁ Фе ргодис{ мл оег зоЁ\аге, Баг4умаге, @есоттитсайоп5; 

— Та|огез ое зо\аге аррИсайоп дигио Ше ицепде4 пе; 

— Фе $1о\ме4 до\уп ите о? Фе зоЙ\аге рго4ис( ап Фе 4е]ау ш Ше ргезещайноп оЁ пцегте(ае ап ошрие 
шРюппайоп; 

— шсотшр!ае 41$р1ау оЁ погтайоп; 

— шсоп$1${епсу оЁ 5юге4 ааа ап шогтайоп ещегеа Бу Фе орегафог; 

— [055 орт@еуапсе ое шЮгтанНоп; 

— у1о|айоп ое сопН4епнаШу о? шЮгтайоп. 

1 аа оп © зисВ “ргитагу” ацаШу даа сое диесЙу Кот Фе сопзитег, деуе]орегз изе “имегпа” ш91- 
сабог$ ю еуаае Ше рагатаег$ оЁ Ве сиггепй ргоесё: 

— Фе Ппез оРсо4де ш Фе зап даг тоде; 

— Фе питьег ор ааес{е4 еггог$ рег 1000 Ппез оЁ со4е; 

— ргоэгат сотр]ехйу рагатеегз; 

— Ме ргоБабИиу оЁ оссиитепсе о? 5реслйс еггогз; 

— Те рго]ес{ сотр]ехйу апа со5{ оРа со4е ипи; 

— рисе ога “тап-топ”; 

— ЗайзИса| сВагафет1$Ис$ о ргосез5ез (ехресаНоп, уапапсе, сотгеайоп Гпсйоп, ес.) ап4 оег езитае4 
рагатетегз. 

Аистетщеа ВеаШу Арр|сайоп Еуашайоп СгКема. Теп ехре$ у тоге Фап 3 уеаг$’ ехрепепсе ш АВ 
аррйсайоп деуеортеп Вер и$ ю стежме стиепа Гог еуаайпо АК аррз. \\е а1у14е4 Фе таш стиепа Гог аззеззте е 
ацаШу оГаиотеще4 геаШу аррИсаНоп$ шю 5 эгопрз (Е1> 1): 

— АК арр 4ез1еп ап4 Ап 4ез1еп; 

— Огаршс ргоэтапил1е; 

— АК аррИсайоп$ ргоэтаплил1е; 

— АррПсайоп ргойПипс ап4 ориплтайоп; 

— РибИ$шие аррИсайоп$ (Би|9). 

\/Веп 4еуеорше ап АК аррйсапоп аеядеп, \е гесотлтепа @е ЮПоулих сВагасцег1$Ис$ Каф звошШ@ Ъе соп- 


' Кип С. Аиотещед ВеаШу Сопбпиит Сопсер апа ВеЕегепсе Моде! — Рай 1: АВС ВеЁегепсе Моде! (\Могк Цет Ргороза!), 1$О-ТЕС ТС 1 $С 
24, 18521-1. 2012. 
2 150 9000 ЕатЙу ОпаШу Мапазетепе. В рз://\у\у\\.150.0г9/150-9001-диаПу-папазетлеп. Вит! 
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514егеа: 

1. Тре загсе. Но\ Фе аррПсайоп сап адарЕ ю уапочз зигРасез. Ш уои изе Нате\могК$, уой сап зе Фе 
Бат за асе д&есйоп Рлиспоп. Еог ехатр/е, {Ве арр|сайоп сап гесозт1те Ше Йоог, \аП$ апа обес. 

2. эБше. Еуаайоп оЁ Шилптшайоп 1$ уегу ппропап Гог Ше геа т оЁ об]ес5. Ц 15 аау1зае 1ю изе упатс 
Попних у звадо\уз ш геа| те. 

3. орасе пеедеа. Озегз сап ехрепепсе АВ ш Штее Ч1Шегеп: 517е3: фа е зсайе, гоот зсае, ап ореп епушгоп- 
тег. ТВе изег зВоц!4 а[\мауз Вауе епочэй зрасе ю зпар/у еп]оу Фе ехрепепсе. Тиз, И 1$ песеззагу ю {фшК оуег {15 
БеГоге фе изег аз изше фе аррПсайоп. Еог ехатшр/е, Ш ап арр|сайоп гедштгез ап ореп епупоптей, уои пее4 ю 
шЮюги Фе аи 1епсе ш адуапсе Беоге Шеу ай изше Фе рго4ис{. 

4. зшфе-изег ог ти-изег ргодис*. [ уоц аге деуорте а ти@-изег име асе, уоц звошШ4 4ез1еп @е ргод- 
ис \Пеп еуегуопе 15 шуо[уе4 ш Фе ргосез$. [Ц 15 ипроцап{ ю стеже а зепзе оР ая епсе соппесйоп уу Ше рго4исй. 
Еог ехатр]е, 1Ё уоц аге деуеорше ап АК сате, уоц сап ргом1ае а тар а зВо\з Фе [осайоп оР изег$ апа ргоу14е$ 
геа]-{те заби5 поиЯсаНопз. 





Е!о. 1. Еуааноп стгцепа ог АК аррИсаНопз$ 


Тре таш стиепа аге о оЩе4 аз Поз (ТаЫе 1). 
ТаЫе 1 
АКВ Арр Оез12п апа Ап Оез1еп еуаайоп ро 


| Анатоль (ого нае челов Гарт [2 
О ЕО МОЕ 
И О ООО ОО 
ЕО ОО 
ЕЕ ООО ПООИ 
ЕО ОО ОО 


Рипс1р]ез оЁ могкше уу зВадегз ап таета|$, геп4етис Ееабагез 


мы | -|-|М 


Ееаагез Юг сизютитилие то4е]$ апа фехагез, аз \е| аз тайепа|$ Гог ехрог {0 Фе гате 


епотте 
Стеже ЦУ-5сап ш тодеги едЦог$ 
Техцлте п1оде]5 ш тодеги е 4 Ког$ 


В1еошФ ап4 апипайпе по4е]$ ш то4еги е4йог$ 
Сизют7е апипаноп ш сате епете$ 
Стезме ап сизюпитте звадегз, пабета1$ ш зиауайоп ргоэтаил$ ап сате епеттез 


Тре па] уегзлоп оЁ Та е 1 \аз ргорозеа Бу Фе аи®огз, еп И \уаз 415сиззе4 ми {еп ехремз ш АК 4еуе]- 
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ортеп: изше Бгалазюгииас {естаче$. 

То зоуе Фе рго[еилз оЁ а15р1ауше сошр]ех огар1с об]есб, и 1$ песеззагу №0 4еуеор еНесиуе те о4$ ап 
а|>огибил$ Юг ргосеззше огарб1с шРюгтайоп а Фе %асез оЁ шриф епсо41те, гапзКогтайоп ап ппаге ЮгтаНоп. АП 
ОР 15 таКез пр а $е{ оЁ Ъаз1с сошрщег отарШс$ фазК$, зо \е Вауе 14епиНе4 фе таш сотропеп пееде4 ю сгеме а 
Ь1>1-аца\у этарб1с сотропепё (Тае 2). 

Отарс ргостатилит? епа ]ез №0 сгеме у1зиа| еНес{5 гоией Фе 4еуе]ортеп: о{ зва4ег$ изше фе Кисйоп$ 
ОР этаршс Пбгамез, сазюпитлие Ше тепдегио оЁ фе деуортеп епутоптепь, изше еНесйуе розёргосеззште Пбгаг- 
1е$, сазюпилте звадо\/з, ап тоге. АП езе сотропеп таКе И роз е ю оМаш №16 -диаШу птасез Бу арр!уше 
уаг10и$ еНес{5, аз \еП аз пиргоуе фе орипитайоп оГарр|сайоп ре{огтапсе. 

ТаЫе 2 
Отаршс ргостатиитите еуалайоп ро 


Г ше оне веент — [о 
о бжиуам иена ие ион [о 
о або овен [о 
Ре обомбение проно Беван [о 
ОЕ ООО ОЕ 


42| 42| || м 


| РомгехоРттрететите Иена вла вар оп во ней ват ое а атеетрвНти [9 
—______  Собаио чане ар дула ела и волне атоме [о 
КО ЕО ОЕ 
—_________ Съомиоромрюбозн ай ИмГпивеарание [о 


Еуааноп рош 6 ргеземе ш ТаЫе 2 ууеге ргорозеа Бу Ше аа ог$ ап4 1аег Феу \уеге 415сиззе4 ми {еп 


| № м 





| 


ехре{$ ш АК 4еуе|ортепи. 
Рипис Фе 4еуе!ортепЕ оЁ Фе зой\аге рай оЁ е то БШе арр|Псайоп \И аиотеще4 геаШу, И 1$ песеззагу 
{0 ргес1з@1у а&щегите Ше Рапдатета| деуе]ортеп: 10015, ап4а, ассог4те]у, Ве ргостатитисе |апеиарез, ргоотатил? 
епугоптеп$, раНегиз, агспиесате, ес. ТЬ1$ 15 песеззагу Гог а с1еагег ап тоге соогдшае4 {еат \огк, ипаег{ап те 
Фе гедлитетеп ап фазК$, сотргецепзоп оЁ Ше 1021с ап4 зедаепсе оР арр|йсаноп деуеортепь, опешщаноп ш Фе 
рго]есф, гедистпе т15Кз, ап ицегасйоп о{ еасВ {еат тетбег у еасВ офег. То астеуе 1$, ме Вауе 14епйНе4 Фе 
РоПо\лие рот Юг Фе ргоогатиие рВазе, уу1сВ зо Ъе ЮПо\е4 ш ог4ег фо оБа а розтйуе гезий о{ Фе 4еуе]- 
ортепЁ ргосез$ (ТаЫе 3). 
ТаЫе 3 
АК АррИсайопз$ Ргоотатилие еуаайоп ро 





Тре еуааНоп ро ргезещеа ш ТаЫе 3 \еге ргорозеа Бу Фе аи®огз ап Бгалазюгтей у {еп ехрем$ ш 
АКВ 4еуеортепе. 
\/Веп 4еуе!орше арр|саНоп$ УИ апотеще4 геаШу, И 1$ пироцап © сопз14ег Фе Геагез оЁ аррИсайоп 
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орИпитаНоп ог РС ап тоБ!е 4еу1сез, аз уе аз Ше агспКесаге оР тобШе деу1сез ш Ше сощехЕ оЁ аррИсайоп орН- 
п7айоп. @ 1$ гесаг4, \ме Вауе 14епи Вей Фе рот Фа{ \ме аду1зе ю ааБеге ю дигие Фе 4еуортеп ш ог4ег © 
астеуе тахилит орйпл17айоп оЁ АВ арр|сайоп$ (ТаЫе 4). 
Тье рго ег 001 ргоу14е$ зрес1Яс даа оп гате ре{огтапсе ап4 Тас|иаез 15 орипитайоп ргосез$. ТВе рго- 
ег ргоу14ез Кате-Бу-ЕКате тейлс$ уу уве уоц сап тоге еазПу 14еп Фу рго ет агеаз. 
ТаЫе 4 
АррИсайоп РгоНПп® апд ОриптайНоп еуаайоп ро 


—__ ОртраотоГАВ аррйсавоп реноме |3 
Ооо о 
буи о 


[—— оримравот опеке вла ыпонаь отеле рано |2 
—_ Орванот ого пай аррйеаноп росе | 


ТаЫе 4 уаз сгеме4 Бу Ше аи®ог$ ап еп Ё \аз 415сиззе ум {еп ехре {5 ш АК деуортете. 





АЁег сотр]ейоп оГаП $асез оГ Фе аррПсаНоп 4еуе!ортепь, # 15 ппроцап фю рибШ$В е аррИсайоп. ТВе 
раб ПсаНоп гергезет5 Фе г@еазе оР Фе арр|сайоп оп апу р|аогт уПеге Фе сизютег (еп4-изег) сап еазПу 4о\уоа4 
Фе Нпа| уегзлоп о? ргодис®, хе асаиаииеа уу 1, сеё аП {Пе песеззагу дЧоситетщаноп$, феси1са1 зиррой ап4 Геедбаск 
Кот Фе деуе]орег. Еасв ое р!аНогил$ риёз Ююгу\аг4 15 зрес1с гедигетеп 5 0 Фе ри 15511 аррИсайоп, св аге 
песеззагу Юг соггес{ 415р]ауте Фе аррПсайоп ш Фе р1аМогил’$ тагКеф, оМаште аП шРЮгтайоп абойй Фе аррйса- 
Ноп’5 орегайоп, епзагше еп4-изег зесигиу, рготойи® @е арр|сайоп, ап4 тоге. [шп 1$ гегага, фе раб Псайоп 15 опе 
ог Ше п105 ипропапё ап сгасла| асе оЁ фе 4еуеюортепе. У!е Вауе 14епиНе4 зеуега| тала р]аогилз$, Беоге изо 
\сь и 1$ песеззагу ю гатШагте уопгзе ул аП Фе доситешайоп ап4 аззет у Еагез Юг Фе арргориае р1а+- 
Гогт (ТаЫе 5). 

ТаЫе 5 
Риб1$Шиз аррПсайоп$ (Би 9) еуаайоп ро 


Риб 15$ Шие арр|сайоп$ (Би! 9) 9 | 
Ееафиге$ ор Баас ап аррИсайоп Юг \\!1140\%5 


Ееафиге$ оЁ бе Би! 4 аррИсайоп Юг Апдго19. 
Ееафиге$ оЁ бе Би! 4 аррйсайоп Юг ОЗ Х/ 105. 


Еуашаноп ро! ргезеще ш ТаЫе 5 \еге ргорозе4 Бу Ше ехремз ш АК 4еуе|ортепе. 





Сопси$10п$. Мо\адауз, Ше аиотешщеа геаШу 15 опе ое п1о5 шпоуайуе ап4 а пе\ 41еа1 тепа ш Фе 4е- 
у@оршсе аррПсайопз$ Юг аП туре оЁ деу1сез. ТВе АК {есбпо]оэу орепз$ а пе\ Вот17оп ап4 1$ сое фо её тоге роршаг 
ш Фе Югезееа е Клбиге. 

Тре стцепа \ме Вауе д4еуеоре4 юг еуаайи® аррПсаНоп$ УИ аиотеще4 геаШу епа е деуеорегз ю сгеме 
Фе о\п зоЙ\аге ргодис5 ш засе$, Базе оп %ер-Бу-%“ер гедитетеп$ Гог Фет. Т1$ м аПо\ уоц 0 сгезже 121 - 
ацаШу зоЙ\аге рго4ис{5 изше фе запага1те4, то4деги 4еуеортеп {юо[5. 

1 ааЧюоп, деу@орегз у Вауе а деае4 ипдегубап те о{ еасВ асе оЁ сгеайпе Фе аррПсайоп ап ®е 
песеззагу 4еуеортеп: {1001[5 ап4 1еспиооз1ез 1ю оМаш Ше В121е${ ачаШу геза. 'ТБаё Ш олуе ап оррогаииу ®ю 4е- 
с14е оп зресйс 4еуеортеп: 10015, тео4$, по4е|$ ап4 1есрпооглез БеГоге загипе \уотК оп а рго]есй; апа, ш Фе 
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оиоШу уепФу Ще ргорозе4 сгцепа, аа4юопа] {езНпе у Бе гедаше4, уПеге тоге зоЁ\аге Не ехреи$ зВоч4 Бе 
шуо[уе4. Епгаета® ита-рацу ехрег5 \Ш аз5е5$ Фе зицаб у оЁ Фе ргорозе4 стцепта. \Уе Бейеуе ай 1$ арргоасВ 
\Ш Бесоте рай о? Фе АК аррПсайоп деуеортеп: ргосез$. 
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Введение. Пакеты верификации проектируют и анализируют корректность параллельных и распределенных 
систем в рамках различных классов темпоральных логик линейного и ветвящегося времени. В статье рассмат- 
ривается полиномиально реализуемый класс ЛоГ-формул, интерпретируемый на многосортных моделях с 
иерархическими надстройками. Структура надстройки описывается произвольной контекстно-свободной (КС) 
грамматикой. Предикаты и функции сигнатуры модели интерпретируются на исходном КС-списке, достраива- 
емом в процессе интерпретации. 

Материалы и методы. Для теорий из ДоТ-квазитождеств, обладающих свойствами нётеровости и конфлюент- 
ности, строится константная модель. Рассматриваются формулы многосортного языка исчисления предикатов 
(ИП) 1-го порядка с переменными сорта «список», интерпретируемые на моделях с иерархизированной 
надстройкой. Теория интерпретируется на деревьях вывода грамматики, описывающих поведение специфици- 
руемой системы. Правила КС-грамматики иерархизируют пространство действий моделируемой системы. От- 
мечено, что для моделирования систем реального времени недостаточно выразительных возможностей ЛоТ- 
формул. Поэтому для спецификации используются выражения с неограниченным квантором всеобщности У, 
известные как ПТ-формулы. 

Результаты исследования. В качестве примера приводится логическая спецификация автоматизированного 
комплекса, который состоит из манипулятора, обрабатывающего детали. Положение позиций фиксируется дат- 
чиками. Описывается цикл работы манипулятора. Спецификация его функционирования состоит в иерархиза- 
ции действий правилами КС-грамматики и их описании формулами ИП 1-го порядка с учетом значений време- 
НИ. 

Обсуждение и заключения. В статье показано, что класс рассмотренных формул можно использовать для моде- 
лирования систем реального времени. Приводится пример логической спецификации управляющего устройства 
поведением манипулятора. 


Ключевые слова: логическая спецификация, модель теории, реагирующая система, КС-грамматика, формула 
ИП 1-го порядка. 
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угастаге 1$ 4езсгеа Бу ап агогагу сощех{-Нее (СЁ) огататаг. ТВе ргеФсайез$ ап РапсНоп$ оЁ фе тоде! $1епабиге 
аге ицегргае4 оп Фе отетпа| СЕ-1$6 \Шсь 15$ сотр! ае4 4игте Фе ищегргеайоп ргосез$. 

Маепта[5 апа Мефо4$. А соп$апе тоде| 15 сопзгисе4 Гог Феопез Нот ЛоТ-ацазидепийез ут Моеепап апа 
сопНчепсе ргорегиез. У!е сопз14ег юга оР Фе ши!@-зоцеЯ Нг$-ог4ег ргефсае са1сиз$ (РС) 1апеиаге ми 
уама ]ез оГ Ше “15” зо ищегргае4 оп то4е|$ уу а Шегагс1те4 зизрепз1оп. ТВе Феогу 15$ пиегргае ш {егил$ ог 
отаплтлаг шЕегепсе еез дезстте фе Бебамтог оЁ Фе зрес1Не4 зует. Тве СЕ-сотаплтлаг тез №егагсЬ1те Фе асйоп 
зрасе ог ше тодеа зует. И 1$ пое4 фаЁ Ше ехргезууе сара Шйез ог А, Т-Ююгич[а$ аге шзиЕЯслет ог тодейп® 
геа]-{ите зузетл$. ТБегеЮюге, ехргеззлоп$ УИ ипбоипае4 итуегза! дцапийег \, Кпо\ги аз РТ Юпиша$з, аге изе4 Гог 
Фе зрес1ЙсаНоп. 

Кезий5. ТБе 1оглса| зрес1Ясайоп оЁ ап ащотае сотр! ех \/№1сВ соп$1${5 оРа \огкртесе таптрщаюг 1$ о1уеп аз ап 
ехатр/е. ТВе |осайоп оЁ Ше роз1оп$ 15 Нхеа Бу зепзог$. ТВе орегайпе суфе оЁ фе таптриаог 15 4езсте4. ТВе 
зрес1Нсайоп оР 1$ орегайоп соп51545 ш фе ШегагсЬ17тайоп оЁ асйоп$ Бу Фе пез оЁ Фе СЕ-отаттаг ап4 Фе 
дезсириоп Бу Фе ЁВгз{-ог4ег РТ-Ююгиач]а$ фаК® шю ассоций Фе ите уашез. 

Оксиуятоп апа Сопс[ияопб. ТЩе рарег зВо\з фай Ше с1а5$ оЁ е сопз14егеа Гогич]а$ сап Бе изе4 {№0 то4е] геа]-ите 
5узеп15. Ап ехатр!е оЁ Ше 1о21са| зрес1Исаноп оРа татрщаюг Бевау1ог сопго] 4е\1се 1$ слуеп. 


Кеунога$5: 1ослса[ зреслйсаНоп, еогу то4е|, теасйуе зузет, СЕ-отаплтаг, Яг${-ог4ег РС-ЮгиШа. 


Еог сйапоп: У. М. ОТа$ЮКоуа, К. 5. Когоуша. Ро]упопимаПу сотриаЫМе »—зрес1ИсаНопз$ оЁ шегагс1са| ппо4е]$ оЁ 
геасип® зузеп5. Адуапсеа Епошеегие Кезеагсв, 2020, уо1. 20, по. 4, р. 422—429. В@рз://Ао1.0го/10.23947/2687- 
1653-2020-20-4-422-429 


Введение. Математически обоснованные практически значимые методы верификации сложных про- 
граммных и технических систем базируются на аппарате математической логики [1-4]. Известна методика 
применения этого подхода для различных видов реагирующих систем реального времени (протоколов связи, 
систем управления, встроенных бортовых систем космической техники и др.). Данная методика позволяет ве- 
рифицировать системы то4е| свесКла® [5—7]. Многочисленные пакеты верификации поддерживают проектиро- 
вание и анализ корректности параллельных и распределенных систем в рамках различных классов темпораль- 
ных логик линейного и ветвящегося времени: ГТГ, СТГ, ТСТЕ и др. [8]. 

Для моделирования времени в этих системах используется стандартная модель временного автомата. 
Это конечный автомат, снабженный переменными специального вида — локальными часами. Количественный 
анализ временных характеристик системы затруднен сложными экспоненциальными алгоритмами построения 
часовых поясов как классов эквивалентности [9]. Поэтому для упрощения анализа необходимо разработать бо- 
лее выразительный практически значимый языка спецификации. 

Предлагается использовать для моделирования язык У-спецификаций, выделенный в концепции семан- 
тического программирования, которая основана на теоретико-модельном подходе [10]. В этом случае для по- 
строения формальной модели анализируемой системы можно использовать язык исчисления предикатов 1-го 
порядка, расширенный аксиомами для операций над списочными структурами ”””. 

Материалы и методы. В работе используется терминология статей [11, 12]. Пусть М — многосортная 
модель сигнатуры < =< [, С,Е, Ё>. Здесь [Г — множество сортов, включая сорт «список» (15. С, Е, В — 


множества констант, функций и предикатов соответственно. Все символы сигнатуры имеют тип. Если {ЕЁ яв- 


ляется и-местной функцией, п > 0, то ее тип равен <, ... 1, Г>, где й, ... и. 1=Г причем й, ... й — типы аргу- 
ментов, { — тип значения функции. Аналогично и-местный предикат гЕР имеет тип < й, ... >. Носитель моде- 


ли М — индексированное семейство множеств (,=С, ,/Е/, где С, — множество констант сорта }; 


О: ЖехО, = О Ре Ч. ХьхО. , 


1 
Для модели М над множеством констант С формируется списочная надстройка 2. (С) из иерархизи- 
рованных КС-списков, структура которых задается КС-грамматикой С =(М№,Т,Р). Здесь №, Т — множества 


нетерминальных и терминальных символов. Множество О.(С) определяется как наименьшее множество всех 
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списков <1 ...., „>, формируемых для каждого правила А -> Х\...Х, ЕР, и>1 следующим образом: Е — произ- 
вольная константа из С’, ‚ если Х; ЕТ; в противном случае для Х; ЕМ элемент Е, — произвольный список сорта 


ре 


Ао-формулы определяются традиционным способом как формулы сигнатуры с с использованием всех 


логических связок (-, Л, \, —) и ограниченных кванторов Ухеф, херь \Ухеь 3Зхс. Здесь х — перемен- 


ная произвольного сорта; Е — терм сорта Д5Е , не содержащий х; е — отношение принадлежности списку; < — 
отношение вложенности для списков, определяемое как < 0 ,..., 0, > < <а4.,....@, >, Т2П. 


Ниже будем использовать лишь ограниченные кванторы вида Ух Е у, Ух Е у, где переменная у сорта 
[51, отношение Е — транзитивное замыкание отношения Е. Обозначим через Х индексированную последова- 


тельность переменных х;, через = — отношение принадлежности или его транзитивное замыкание. 

Правила КС-грамматики иерархизируют пространство действий моделируемой системы. Из соображе- 
ний эффективности вычислений выделяется класс АоТ-формул с «древесным» префиксом. Введем для элемен- 
тов списков отношение = — «левее», а именно, для списка < .-: ©,В.... > считаем © - В. 

Определение. Ло-формула вида 


(Ууем)... (Уз Е») (п: < П)... (и, < Ь) ФГ), т> 1, р> 0 
называется Ло1-формулой, если п;,[; Е (9, Т), 1 < /< р; для всех переменных префикса выполняется условие: 
ТЕТ ,1 << т или 7141 = Фь, К < Е. Если Т; +1 = Ук, ТО 4 к 3 Тки ик = Т, для всех К < [. 

Легко показать, что префикс АоТ-формулы в силу ограничений на переменные можно представить в 
виде дерева с корнем 71, вершинами "%,, т; и дугами, идущими из вершины 7; в вершину 17;, 1 << т. 

Выразительных возможностей Ло/-формул недостаточно для моделирования систем реального време- 
ни, циклически функционирующих неограниченно долго. Будем использовать для спецификации ПТ-формулы 
с неограниченным квантором всеобщности У. 

Определение. Формула, полученная из ЛАо1-формулы Ф навешиванием неограниченного квантора все- 
общности Уп Ф(и), называется ПТ-формулой. 

Модель М определяется теорией из квазитождеств вида: 

(Ум, Е?))... (Ут, Е т) (< Ц)...(пь < 5 )(ф(®, т) > ф(®,Т)). 
Здесь формула ф (\/) — это конъюнкция атомных формул (или их отрицаний) вида г, т. =то, (Ех), ЕР,ТЕ В, 
Т1, 1› — термы сигнатуры о. 

Алгоритм построения модели [11 реализует правило вывода то4и$ ропепз (если фи Ф-—»у, то чи). Вход- 

ные данные для интерпретатора — множество начальных значений функций и предикатов вида: 
50 = Ф(2), (с) = сир ЕР, ДЕЁ}, 
где с — множество констант из и элементов, п > 1. 

Аксиомы (ах) обрабатываются в определенном порядке, сначала с позитивными вхождениями преди- 
катов в левую и правую часть ах, потом с негативными вхождениями до получения неподвижных точек вычис- 
лений. Область определения функций и предикатов, входящих в правую часть ах, расширяется при истинности 
ее левой части. Это объясняется тем, что интерпретатор задает новые значения функций и предикатов так, что- 
бы правая часть ах была тоже истинной. Пусть состояние 5„ анализируемой системы на -м шаге вычисления 
содержит значения всех предикатов и функций сигнатуры модели, а функция т Р(5) > Р(5) в (терминоло- 
гии [3] — преобразователь предикатов) отражает изменение состояния при переходе интерпретатора /иё от и-го 
шага вычисления к п+1. Интерпретатор /и{ строит для монотонного преобразователя т на Р(5) наименьшую 
неподвижную точку 7. 

т(7) =Ч; т ($5), где 10 (7) = 7, < *1(7) = т((7)). 

Формально функции } Е РЁ и предикаты р Е Ё интерпретируются на КС-списке (п), представляющем 
дерево вывода {"(п) в грамматике С, где п — шаг работы /мг. Из-за трудности представления компактной фор- 
мы записи алгоритма интерпретации на элементах КС-списка приведем сначала словесное пояснение алгорит- 
ма, ориентируясь на дерево {"(п). Входные данные для алгоритма построения модели (18): исходное дерево 
вывода {7(0) в грамматике С, расширяемое в процессе построения модели М и Рас! =5°. КС-грамматика ис- 
пользуется в процессе построения модели следующим образом. Во-первых, правила Р иерархизируют про- 
странство действий и состояний анализируемой системы. Будем считать, что имена действий, представленные в 
сигнатуре модели предикатами, и имена соответствующих нетерминальных символов грамматики совпадают. 
Во-вторых, символы из алфавита Г грамматики однозначно определяют сорта всех элементов универса модели, 
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включая списки, которым ставится в соответствие сорт, определяемый меткой корня соответствующего дерева. 
Сорта будем обозначать мнемонично начальными строчными символами имен нетерминальных и терминаль- 
ных символов грамматики с добавлением в конце символа 5 (507). Основное достоинство КС-грамматик состо- 
ит в возможности использования эффективных синтаксически ориентированных (СО) методов для анализа кор- 
ректности (верификации) модели, развитых в теории синтаксического анализа языков программирования. 

Интерпретатор начинает работу с просмотра дерева {7(0) от корня сверху вниз, слева направо. Пре- 
фикс всех аксиом удовлетворяет ограничениям До,ТГ-формул. Сорта префикса определяются символами КС- 
грамматики С. Интерпретатор выбирает в качестве констант области истинности предикатов, входящих в акси- 
ому ах Е ТИ, константы, сопоставляемые вершинам куста дерева, просматривая его сверху вниз, слева направо. 
Причем корень куста помечен именем соответствующего предиката. Для отражения зависимости моделируе- 
мой технической системы от последовательности входных сигналов необходимо достраивать исходное дерево 
7(0). С этой целью последовательность правил рт*(ах) Е Р\ приписывается к выводу дерева, полученного на 
предыдущем шаге работы алгоритма. Причем константы из области истинности предиката т используются в 
качестве терминальных символов, подчиненных в дереве вершине, помеченной нетерминальным символом 7. 

Опишем алгоритм интерпретации [т более формально, не детализируя процедуру Соп(О, ТИ) — полу- 
чения всех логических следствий из множества формул О на основе аксиом теории ТА. ТЛАу = 5%. Теория 
ТАроз = ТА включает только позитивные вхождения предикатов. ТАиед Е ТА включает негативное вхождение 
предикатов в правой части аксиом. Обозначим ох, нед — деревья выводов, генерируемые в процессе интер- 
претации. 
9:=9; 

О’:=ТИо; 
\ПШе 0 = 0'4до 
Чроз = ©; 
О в - 0’; 
У\ВШе Оо; == Оро; 40 
роз у ре 

О раь == Сот(Оъо«, ТПроз) 
епа \Ю|е 
гебагт (Ороз , Ироз) ; 
пед - 0; 
О Е = Ороз; 
У\ВШе Опед 3 Опед 90 
пед = Опед; 

О пед т Сопт(Олед, ГТтеа) 
епа \№|е 
гебагт (Опед, тед)} 

а = Чроз; 
0’ = Опед 

епа \№|е 
геаги (Опед» пед) 

Верификация модели ЛМ состоит в проверке свойств, которым должна удовлетворять анализируемая 
система. Будем выражать эти свойства произвольными ДоГ-формулами. Используя СО-методы проверки фор- 
мул, можно строить доказательство аналогично [13]. 

Теорема. Произвольная ДоТ-формула с т-ограниченными кванторами всеобщности проверяется на КС- 
списке мощности п за время О(п"*Т). 

Мощность списка &[ равна мощности множества {5 |5 Е Ц}. 

Оценка является верхней, и ее можно понизить до линейной, если проверять формулы с использовани- 
ем специфичных СО-методов обработки языков. 

Результаты исследования. Приведем логическую спецификацию автоматизированного комплекса, со- 
стоящего из манипулятора, обслуживающего технологическую линию (1) с двумя позициями: загрузки и раз- 
грузки деталей (4 и и/А соответственно) [14]. Датчики фиксируют положение позиций. Манипулятор функцио- 
нирует циклически, начиная с позиции загрузки. 
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ЦИКЛ 

1. В начальном положении /[4 для загрузки детали манипулятор поднимает электропривод за 4 сек. 
Сжимает автоматизированные клещи и берет заготовку (2 сек), опускает электропривод (4 сек) и передвигается 
вправо к автомату до срабатывания датчика положения И. 

2. Чтобы установить заготовку на автомате на позиции И, манипулятор поднимает электропривод, раз- 
жимает автоматизированные клещи (2 сек), опускает электропривод. Далее манипулятор ожидает 4 мин, после 
чего повторяет те же процедуры, что и на позиции /4. Затем манипулятор двигается влево к позиции разгрузки 
до срабатывания концевого выключателя и/4. 

3. За 8 сек деталь разгружается на конвейер. Манипулятор передвигается влево до фиксации датчика [4 
на позицию загрузки. Далее процесс работы комплекса циклически повторяется. 

Спецификация системы состоит из нескольких уровней. Поведение манипулятора определяется сигна- 
лами датчиков, фиксирующих его положение: [/4, и, И (> 4, —ша, —И, отрицание указывает на отсутствие со- 
ответствующего сигнала). Эта последовательность сигналов представлена кортежем тс=<х, у, 2...>, где х = А 
(—/4) у =иМ@ (— и), 2=И (— Ц). Она генерируется конечным автоматом с начальным состоянием х. Этот авто- 
мат строится по праволинейной грамматике с правилами: 

51 4 51 | 4$ в 
5#1—> 1 512 | 912 
5$12—> и! 51 | — ша 54. 

Обозначим через 0)51е множество списков, составленных из цепочек символов, порожденных этой 
грамматикой. 

Внешнее дискретное время (переменная и в логической спецификации) определяется количеством пе- 
реходов в автомате. Для описания второго уровня функционирования манипулятора используется КС- 
грамматика, которая указывает последовательность действий (Орег) манипулятора: 

— /,, [а — загрузка заготовки манипулятором в положении [4 и И соответственно; 

— От, Оша — разгрузка детали в позиции и/4 и И; 

— Моуег, Моуе[ — движение манипулятора направо и налево; 

— 1[,5юр, А5юр, Ибюр — остановка манипулятора в соответствующем положении; 

— Вхр — ожидание; 

— Сг — выход из строя управляющего манипулятором устройства. 

На состояния манипулятора (символ Ро5) влияют его действия. В указанном примере состояние харак- 
теризуется непрерывным временем Тйтес и дискретным Тйпе4, задаваемым натуральным числом. Значением 
сорта Типес являются сегменты вида < 1, Ь>, 4, # — константы, причем < заменяется на ( или [ в зависимо- 
сти от того, включена левая граница в сегмент времени или нет, аналогично для >. 

При спецификации поведения манипулятора абстрагируются от значения времени перемещения мани- 
пулятора из одного положения в другое (определяется сигналами с датчиков положения — входными к устрой- 
ству управления манипулятора). Сигналы, которые подаются на исполнительные механизмы манипулятора для 
осуществления передвижений и работы автоматизированных клещей, являются выходными. 

Ниже представлены правила грамматики С. 
‚ а! {Орег} *. 
| Орег->> Г. Га| Оп! Ота| Моуе!| Моуе!| Г5ор | Аяюр| (Лор | Ехр| Ст. 
. [-— 5. 
‚„ Га— 591. 
„ И — 51а. 
‚ Иша > 5е. 
.‚ Моуег — 51ие. 
‚ Моуе[ —> 5кие. 
‚ [5юр —> 5!е. 
10. А5юр — эюие. 
11. Обюр->> зе. 
12. Ехр —> ие. 
13. Си >> 5юе. 
14. кие Ттес Тте4| Тте4 Ттеа. 
15. Тйпес-—> Титеа | (Трпеа, Гтеа) | [Глтеа, Гите4) | (Трпеа, Гте4] | [Глтеа, Гите4]. 
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Птеа — класс лексем, значениями которых являются натуральные числа, вычисляемые во время ин- 
терпретации теории ТИ. 

В теории Тй переменные в формулах обозначаются мнемонично в соответствии с их сортом: р (зе) = 
р (51е), р (орег) = р (Орег), р (п) = р (Ё =р ([Гтеа,, р (сё) = р (Гтес). Предикаты Га, Т/, Оша определены на 
множестве Типеа. Га (п) истинен, если манипулятор находится на позиции загрузки. Аналогично для Т/ (и) — 
на позиции обрабатывающего автомата, О (п) — на позиции разгрузки. Перечислим области определения 
остальных предикатов: [,50р, А5юр, Цб5юр, Моуег, Моуер, Су < ПтеахТГитеа; Г, Га, ЦИш, Иша, ЕхрхТитес х 
Гтеа. В формулах используются стандартные функции йеаа (<х,,...х,>) = х,, @Щ< х,..х, >) =<х,,..х, > и 
функция Мс: Титеа- Ос. Здесь Ос — множество списков из сигналов датчиков. 

В начальный момент времени [ = 0, и = | и на 1-м шаге вычисления выполняется предикат [.50р(0,1); 


* 


Мс (Т) = тс, где тсеПс. В аксиомах 1-11 переменные &, и и сЁ связаны ограниченным квантором \ зе е орег, 
\ Би, СЕ е зи. В аксиомах 12—17 переменная ип связана неограниченным квантором \; 5 = <<<<0,1>>>> — 
начальное значение списка, на котором интерпретируется теория. Составляющие его списки в порядке глубины 


вложенности имеют сорта р (К), р (Оре’), р (Ё510р), р (0) =р (Г) =р (Птеа). Списку 5, соответствует дерево 1, 
(рис. 1). 





Ттеа. 0 


Рис. 1. Дерево Г, , соответствующее списку 5, 





В аксиомах теории справа в квадратных скобках приводится последовательность правил грамматики С, 
достраивающих дерево Г. 
Аксиомы теории 
1. [5юр(@, п) л Га (п) >> МЬЕ+Т), п) ^ Мс (и + Г) = @ Ме (п)) [1; 2.1; 3; 14.1; 15.3]. 
2. Ксь п) ^ ТЦп + Г) -> Моуег (с [2], п + Г) л Ар (с! [2], п + 1) [1;2.5; 7; 14.2;1; 2.8; 10; 14.2]. 
3. Аяюр ( п)-> Иша 6 +7], п) [1; 2.4; 6; 14.1; 15.2]. 
4. Оша (сь п) —> Ехр ((сЁ [2], сё [2] + 180), п) [1; 2.10; 11; 14.1; 15.2] 
5. Ехр (сь п) —> Га [сё [2], с [2] +3), п) ^л Мс (п + Г] =@ Ме (п)) [1; 2.2; 4; 14.1; 15.2]. 
6. [а (с п) ^ Ц (и + Г) > Мое! (сЕ [2], п + Гл Ир (сЁ [2], п + 1) [1; 2.6; 8; 14.2; 1; 2.9; 11; 14.2]. 
7. (15ор (6 п) > Ош (ЕТ), п) ^ Мс (п + ГП = @ЩМС (п) [1; 2.3; 5; 14.1; 15.2]. 
8. Ци! (с, п) л [А (п + Т) >> Моуе! (сё [2], п + Г) л 15юр (© [2], п + Г) [1; 2.6; 8; 14.2; 1; 2.7; 9; 14.2]. 
9. Ош (сь п) л —[а (и + 1) -> Си (сЕ [2], п + 1) [1; 2.11; 13; 14.2]. 
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10. Га (с п) л — Ипа (п + 1)- Си (сЕ [2], п + 1) [1; 2.11; 13; 14.2]. 
11. 6 (6 п) л — 11 (п- 1) —> Си (© [2], п + 1) 1; 2.11; 13; 14.2]. 

12. реаа (Мс (п)) =1а > Га (п). 

13. йеаа (Мс (и)) = — 4 —> — Га (м). 

14. реаа (Мс (п)) =И-—> Т[ (п). 

15. йеаа (Мс (п)) = —И — — Т/ (и). 

16. леаа (Мс (п)) = ша > От (п) 

17. йеаа (Мс (п)) = ша > — ИМ (п). 

Теория Тй обладает свойством нётеровости, т. к. изменение переменной под квантором У ограничено К 
— количеством элементов в исходном списке тс. При этом йеаа (Мс (К-Т)) не определено, т. к. Мс (А+Г=< >. 
Отметим, что цепочки —/А и др. в правой части аксиом 12-17 имеют сорт 517те, и — рассматривается не как 
логическая операция, а как символ. 

Для начального значения функции Мс (1) = <Я, И, ша, 14, — И, шА > получаем множество следствий: 
[510р(0,1), ГА (1), Ё ([0,7), Т), Ме (2) = < #, и, 4, — И, ша>, Моуег (71,2), Азор (7,2), Иша ([7, 14], 2), Ехр([14, 
194], 2), Га([194, 197), 2), Ме (3) = <и4, 14, — И, и 4>, Моуе! (197, 3), И15юр (197, 3), Чт ([197, 204], 3), Мс (4) = 
<1а, — И, и4 >, Моуе! (204, 4), [.51ор (204, 4), Г. ([204, 211), 4), Мс (5) = < — И, иа>, Си (211, 5). 

Полученное множество следствий иерархизируется в соответствии с выводом в грамматике С, полу- 


ченным в результате правил, приписанных к интерпретируемым аксиомам. В соответствии с ними к дереву Г 0 


справа от узла, помеченного символом Орег, добавляются еще 12 вершин, помеченных этим же символом и 
связанных ребрами с корнем. К новым вершинам подвешиваются поддеревья с корнями, помеченными симво- 
лами [4, Моуег, Азор и т. д. с их состояниями и константами сорта р (Титеа), полученными в результате интер- 
претации. 

Обсуждения и заключения. На построенной модели можно проверять истинность произвольных Ло/1- 
формул. Например, формализуем утверждение: «Если манипулятор стоял на позиции загрузки в момент # на 
шаге и цикла его функционирования, то на шаге п + 2 через 197 сек он начинает разгрузку в течение 7 сек». 


Приведенная ниже формула проверяется на заданном списке Орег сорта р (Оре’): 
(У5ище е орег) (УЕ Зще) (Упеза) (Ё51ор(Е, п) —> Ци ([&+ 197, #+ 204], п +2). 
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Введение. Системы электронного документооборота (СЭД) используются для хранения, обработки и передачи 
больших объемов информации. Автоматизация данных процессов является актуальной проблемой, требующей 
комплексного решения. Ее решение позволит снизить временные и материальные затраты на проектирование и 
осуществить переход к более совершенным, адаптивным СЭД. Статья посвящена разработке новых методов 
автоматизации процесса распределения информации в СЭД. 

Целью работы является повышение точности распределения информации в СЭД за счет перехода от 
аналитических или алгоритмических решений к использованию новых методов, функционирующих на основе 
технологий машинного обучения. Использование нейронных сетей для достижения поставленной цели также 
позволит повысить эффективность разработки программного обеспечения за счет автоматизации процессов 
анализа и обработки информации. 

Материалы и методы. Предложен новый метод автоматизированного распределения информации на основе 
технологий машинного обучения, включающий математическое описание правил распределения информации. 
Сформулированный перечень условий распределения информации позволяет реализовать программное 
обеспечение на основе нейронных сетей для решения задачи автоматического распределения данных в СЭД. 
Результаты исследования. Метод автоматизированного распределения информации апробирован на примере 
предметной области СЭД при решении задачи анализа правильности информации, введенной пользователем. В 
ходе экспериментальных исследований получено, что предложенный метод, функционирующий на основе 
технологий машинного обучения, обеспечивает лучшую точность (на 8 % выше) и более эффективен (в 
соответствии с метриками Джилба и цикломатической сложности) 

Обсуждение и заключения. Полученные результаты подтверждают эффективность и точность предложенного 
метода. Представленные результаты могут использоваться для автоматизации процессов распределения и 
проверки информации в адаптивных СЭД, а также других информационных системах. На основе 
разработанного метода также возможно решение связанных задач: поиск дубликатов и похожих документов, 
классификация и размещение по категориям файлов. 


Ключевые слова: системы электронного документооборота, распределение информации, автоматизация, 
адаптивность, машинное обучение. 
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Введение. Системы электронного документооборота (СЭД) широко используются для хранения, 
обработки и передачи информации различного типа: документов, электронных таблиц, графической и 
аудиоинформации, технической и бухгалтерской документации и так далее. Основная функциональность по 
работе с документами за годы развития СЭД сформирована. Дальнейшее развитие направлено на 
автоматизацию человеческой деятельности, повышение гибкости и надежности, реализацию интеллектуальных 
модулей по поддержке принятия решений [1-3]. 

Однако подобная гибкость и адаптивность СЭД приводит к дополнительным временным и 
материальным затратам, увеличивает сложность проектирования и модернизации системы. Поэтому 
актуальным вопросом при разработке адаптивных информационных систем, в том числе СЭД, является 
автоматизация проектирования системы, что включает множество отдельных задач по автоматизации 
процессов анализа, обработки и передачи документов, решение которых в комплексе позволяет снизить 
нагрузку на разработчиков в ходе реализации функций адаптивности в СЭД. 

В рамках исследования будет рассмотрен вопрос классификации информации и ее последующего 
автоматического распределения по заданным категориям в СЭД. Данная подзадача является одной из наиболее 
распространенных при организации электронных хранилищ и архивов, картотек, заполнении форм документов. 
Под распределением понимается размещение данных и файлов на заданных позициях: по категориям, 
носителям, каталогам и так далее [4]. Классические алгоритмы не всегда способны определить ошибку в 
размещении информации без привлечения эксперта или модератора. 

Рассматривая подходы к решению задачи распределения информации, можно выделить подходы, 
основанные на применении машинного обучения'. Использование технологий искусственного интеллекта для 
решения задач классификации доказало свою эффективность в ходе многочисленных исследований и 
экспериментов [5—8]. При решении поставленной задачи и реализации метода распределения информации для 
автоматизации процесса классификации данных будут использованы технологии машинного обучения. 

На основе анализа процесса движения информации можно выделить следующие актуальные задачи по 
классификации и распределению данных в адаптивных СЭД: 

— классификация и определение соответствия введенных данных той категории, в которой их 
расположил пользователь [9]; 


ГОтадеу! $. А загуеу оп даа тише с1аззсайоп аеотйВииз / $. Отладеуь К. 5. 7. Магзейпе // 2017 Пцегпанопа! Сопёегепсе оп 51е па! Ргосеззте ап Сопиташсайоп (1С$РС): ТЕЕЕ, 2017. 
Р. 264-268. 
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— классификация данных по категориям” и, в случае ошибки, перемещение их в правильные категории 
либо выдача пользователю предупреждения [10]; 
— определение дублирования данных по признакам [11]. 

Для решения поставленных задач предлагается разработать метод автоматизированного распределения 
информации в адаптивных СЭД, обобщающий существующие подходы к классификации информации и 
основанный на применении методов машинного обучения для автоматизации процессов обработки данных. 

Материалы и методы. Формализуем основные этапы метода автоматизированного распределения 
данных в адаптивных СЭД. 

В основе метода — формирование и обучение нейронной сети для классификации информации. 
Поэтому на первом этапе реализации метода необходимо подготовить множество информационных объектов 
Х ={х,..ху} для обучения и тестирования нейронной сети. Объект х, е Х может быть представлен текстовой 
или численной информацией, введенной пользователем в поля форм, либо файлами, загруженными через 
интерфейсе СЭД [12]. Процесс подготовки данных может включать нормализацию, токенизацию, 
лемматизацию, извлечение свойств и атрибутов файлов. Для собранных данных задается множество категорий 
У={у.....Уи!, то есть существует непрерывное отображение Х —>У. Аппроксимируем это отображение 
нейронной сетью №№: 

ММ(Х) =У. 

При достаточном объеме тренировочных данных можно обеспечить необходимую точность 
классификации. В соответствии с исследованиями, представленными в работе [13], на каждый выходной 
признак необходимо предоставить не менее 50 экземпляров тренировочных данных. Тогда для Ух, Е Х 
получим М№ММ№(х,) = у, , где 1=1..М. 

На второй этапе метода осуществляется распределение информации. Пусть в СЭД задан набор 
информационных объектов Х = {х,....х,} и множество соответствующих им категорий У =\{),....уи}. Для 
распределения информации необходимо осуществить проверку следующих условий: 

— если Ух, Е Х, ММ№(х,) = у, и объект х, помещен пользователем в категорию у,, тогда положение х, 
остается без изменений; 

— если 3х еХ, №ММ(х) = у, (где у, — категория вредоносных объектов), то объект х, удаляется из 
СЭД, администратору системы отправляется предупреждение, пользователь, добавивший объект х , 
блокируется; 

— если хе Х, М(х,) =у, , но объект х, введен пользователем в категорию у, (у, =у, ), то 
необходимо осуществить перераспределение информации в СЭД. Рассмотрим все возможные варианты 
перераспределения: 

1. Категория у, свободная (пустая), тогда объект х, переносится из категории у, в у,. 


2. Категория у, занята некоторым объектом х, (таким, что М№ММ№(х,) = у,), тогда объект х, переносится 


в буфер, пользователю выдается предупреждение о неправильном вводе х, и дублировании информации. 


3. Категория у, занята некоторым объектом х, (таким, что №М№(х, ) = у, ), тогда объект х, переносится 
в буфер, а объект х, переносится из категории у, в у, ‚ пользователю выдается предупреждение о 
неправильном вводе х, . 


Рассмотренные в рамках метода этапы и условия распределения информации включают основные 
сценарии добавления информации в СЭД. За счет применения технологий машинного обучения возможно 
автоматизировать распределение информации по заданным категориям. 

Результаты исследования. Для реализации и интеграции метода в адаптивную СЭД будем 
использовать микросервисный подход. Нейронная сеть реализуется посредством библиотеки Кегаз (Ру®оп), 
после этого импортируется в микросервис, который может быть реализован на основе любого фреймворка, 
например, Е1азК. Данные между СЭД и микросервисом передаются по протоколу НТТР в формате ЛЗОМ, что 
обеспечивает совместимость с любой реализацией СЭД [14]. 

В качестве примера задачи распределения данных будем использовать форму карточки документа в 
СЭД с 5 полями: наименование документа; автор документа; контактная информация (адрес) автора; дата 
создания; описание документа. В ходе экспериментального исследования будет моделироваться ввод данных 
пользователей в форму и заполнение информации в указанные 5 полей. Нейронная сеть будет 
классифицировать введенную информацию и проверять ее на соответствие заданной категории. 


? Попова Е. С. Использование искусственных нейронных сетей для решения задачи классификации текста / Е. С. Попова, В. Г. Спицын, Ю. А. Иванова // Труды 29-й междунар. конф. по 
комп. графики и зрению «Графикон». Брянск, 2019. С. 270-273. 
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Для формирования набора исходных данных используем генератор на основе открытых баз русских 
имен, городов и стран, а также последовательностей из русских слов и цифр. Это позволит получить 


конструкции, соответствующие реальным данным, но не хранящие персональные данные реальных людей. 
Текстовые данные будут обработаны с помощью токенизатора для перевода их в численный формат. 


Максимальная длина предложения ограничена 20 словами. Сгенерируем тренировочную выборку. Для каждой 
категории сформируем по 10 000 элементов для каждой категории. Также добавим 50 000 элементов с 
неправильными и ошибочными данными. Таким образом получим массив из 100 тысяч элементов. 


Зависимости между входными признаками и категориями представлены на тепловой карте на рис. 1. 


Карта позволяет визуализировать матрицу корреляций и провести наглядный анализ данных для определения 


того, как параметры влияют друг на друга и на выходные переменные [15]. 
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Процесс обучения сети отражен на рис. 2. Итоговая точность нейронной сети на контрольном наборе 
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Рис. 1. Тепловая карта исходных данных 


после 5 эпох составила 97,8 %. 
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Рис. 2. Процесс обучения нейронной сети классификации 
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Количество ошибок не превысило 7 % для категории «Автор» и 1 % для категории «Описание». В 
остальных категориях некорректные данные были распознаны в 100 % случаев. 
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Рис. 3. Тестирование нейронной сети на ввод некорректных данных 
Второй эксперимент (рис. 4) заключается в вводе данных из / категорий в 1 (7=1). Формируется 


набор из 100 элементов для Г каждой категории, состоящий из элементов других категорий (по 25 каждой 
категории). График отражает количество тестов, в которых нейронная сеть неправильно распознала данные 
других категорий, как соответствующие текущей категории. Во втором эксперименте был допущен 1 % в 
категории «Наименование», остальные категории отработаны без ошибок. 

Далее сравним точность предложенного метода (далее обозначен как «нейросетевой метод») с 
классическим решением задачи классификации и распределения информации (обозначен как «классический 
метод»). Для классического метода получены следующие результаты: 16 % ошибок в категории 
«Наименование» и 25 % — в «Описание», что существенно хуже показателей нейросетевого метода (рис. 5). В 


среднем использование нейросетевого метода обеспечивает повышение точности на 8 %. 
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Рис. 4. Тестирование нейронной сети на соответствие категории 
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Рис. 5. Сравнение точности нейросетевого и классического методов 


Далее осуществим сравнение методов по метрикам сложности реализации программного обеспечения. 
В качестве метрик используются цикломатическая сложность, выраженная числом МсСабе (общая и 
усредненная) [16], а также сложность по метрике Джилба [17]. Итоговые результаты представлены в таблице 1. 


Таблица 1 
Сравнение классического и нейросетевого методов 


Метрика Классический метод Нейросетевой метод 


Цикломатическая сложность (средняя) А (3,71) А (1,37) 
Метрика Джжилба 





Таким образом, сложность нейросетевого метода по метрикам цикломатической сложности и Джилба 
ниже. Точность разработанного метода выше почти на 8 %. Стоит заметить, что сложность реализации 
классического и нейросетевого метода классификации и распределения данных относительно невелики, так что 
их можно считать сравнимыми. Однако, полученные результаты по точности классификации подтверждают 
большую эффективность нейросетевого подхода. 

Обсуждение и заключения. В работе поставлена задача автоматизации процессов классификации и 
распределения данных в адаптивных СЭД. Предложенный метод автоматизированного распределения данных 
включает формализацию процесса классификации и распределения данных, использование технологий 
машинного обучения для автоматизации решения задачи. Сформулированный перечень условий распределения 
информации позволяет реализовать программное обеспечение на основе нейронных сетей, решающее задачу 
автоматического распределения данных. 

Для апробации разработанного метода проведены экспериментальные исследования на базе 
стенерированных данных о документах в СЭД. Точность обученной нейронной сети составила около 98 %. 
Дополнительные испытания показали, что нейронная сеть может обнаруживать некорректно распределенные 
данные почти в 100 % случаев, причем, при наихудших условиях ошибка не превысила 7 %. Таким образом 
подтверждена эффективность и точность предлагаемого метода. Разработанный метод, по сравнению с 
классической реализацией на основе алгоритмического обеспечения, показывает следующий положительный 
эффект: повышение средней точности на 8 %, снижение сложности реализации. 

Представленные результаты могут использоваться для автоматизации процессов распределения и 
проверки информации в адаптивных СЭД, а также других информационных системах. Кроме того, на основе 
предложенного метода возможно решение связанных задач: поиск дубликатов и похожих документов, 
классификация и размещение по категориям файлов. 
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Введение. Работа посвящена численному решению задач динамики фитопланктона в прибрежных системах 
на примере Азовского моря. При этом учитывается трансформация форм фосфора, азота и кремния, а 
также движение водной среды, распределение температур и соленостей по акватории моря. Речной сток, 
меняясь по объему и по химическому составу, существенно влияет на изменчивость гидрофизических и 
биогеохимических параметров процессов, происходящих в прибрежной системе. Этим обусловлена 
необходимость статистической обработки данных многолетних наблюдений за характеристиками речного 
стока. 

Материалы и методы. Математическая модель биогеохимических циклов основана на системе 
нестационарных уравнений конвекции — диффузии — реакции параболического типа с нелинейными 
функциями источников и младшими производными, к которой добавлены соответствующие начальные и 
граничные условия. В ходе статистического анализа рядов многолетних наблюдений за стоками рек 
найдены значения следующих показателей: коэффициенты асимметрии, эксцесса, дисперсия и стандартное 
отклонение, коэффициент вариации, коэффициент автокорреляции, отношение Неймана и критерий 
Андерсона. 

Результаты исследования. Статистический анализ рядов многолетних наблюдений за гидрохимическими 
показателями реки Дон свидетельствует о неоднородности натурных данных. Это связано со 
стохастичностью поступлений биогенных веществ и объема пресноводного стока в море в результате 
воздействия природных и антропогенных факторов. Натурные данные целесообразно соотносить с 
сезонными изменениями температуры водной среды. В работе представлены результаты вычислительного 
эксперимента по моделированию динамики фитопланктонных популяций в летний период, когда 
температуры благоприятны для их размножения и роста. Предложенная математическая модель учитывает 
пространственно-неоднородное распределение и трансформации форм фосфора, азота и кремния, а также 
изменения солености, температуры и движения водной среды. 

Обсуждение и заключения. Рассмотрена многовидовая математическая модель динамики фитопланктонных 
популяций с учетом трансформации форм фосфора, азота и кремния в прибрежных системах. Проведен 
анализ данных натурных наблюдений, для чего рассчитаны его основные статистические параметры. В 
итоге сделан вывод о существенной изменчивости данных многолетних наблюдений. Это обусловлено 
двумя причинами. Первая: случайный характер поступления биогенных веществ и объемов стока рек в 
результате воздействия антропогенных факторов. Вторая: чередование в последние 12-15 лет 
относительно многоводных и маловодных периодов по пресному стоку. Гидрологический режим меняется 
главным образом за счет сокращения среднегодового пресноводного стока Дона и частично Кубани. Эта 
тенденция, по-видимому, будет усиливаться из-за климатических изменений, а также при дальнейшем 
регулировании стока Дона после ввода Багаевского гидроузла. Численные эксперименты на основе 
натурных данных подтвердили предсказательную ценность разработанных моделей и программ. С их 
помощью можно прогнозировать изменение состава и численности (концентраций) в Азовском море 
основных планктонных популяций, которые определяют, с одной стороны, кормовую базу, а с другой — 
состояние водной среды в условиях продолжающегося осолонения моря. 


Ключевые слова: биогеохимические циклы, фитопланктонная популяция, биогенное вещество, химико- 
биологический источник, уравнение конвекции — диффузии — реакции, натурные данные. 


Информатика, вычислительная техника и управление 


> 
9% 
1 


ЬИр://уезш-допза.га 


м 


А4уапсеа Епотеетпи? Кезеагсй 2020. Т. 20, №4. С. 437—445. 155М 2687-1653 





Для иитирования: Моделирование биогеохимических процессов в Азовском море с использованием 
статистически обработанных данных о речном стоке/ А. И. Сухинов, Ю. В. Белова, А. В. Никитина, 
А. М. Атаян// Адуапсеа  Епошеегие — Кезеагсй. — 2020. — Т.20 №4. — С. 437- 


445. № рз://401.0г2/10.23947/2687-1653-2020-20-4-437-445 


Финансирование: работа выполнена при финансовой поддержке РФФИ, проект № 20-01-00421. 


© Сухинов А. И., Белова Ю. В., Никитина А. В., Атаян А. М., 2020 


Мод4ейте ЫМосеосвеписа| ргосез$е$ ш {Пе Ал7оу ЪЗеа изто убай$ЯсаПу ргосе$$е4 даа оп ттуег 
Поуу 


А. 1. биквтоу", У. У. Веоуа', А. У. МИЧ@та”, А. М. А'ауап' 


'Роп Эа е ТесЬшса1 ОшуегзИу (В озюу-оп-Ооп, Визап Еедегайоп) 
* Зирегсотриегз ап Меигосотриег$ Везеагсь Сещег (Тазапгое, Визчап Еедегайоп) 


птодисйоп. ТЬ1$ \о!К 1$ аппте4 а зо1уше фе рго ет о? рвуюр!апКюп дупапс$ ш Фе соаа| епугоптет$ ие 
Фе ехатр[е о Фе Атоу Зеа. Т№$ фаКез шю ассоииЕ Фе тапзогтайНоп о Гогил$ оЕ рвозрВоги$, пигосеп ап $Шсоп, 
аз уе аз Ше адиайс тедта тоНоп, Фе 415и1БиНоп о{ {етрегаагез ап заПи1ез оуег Ше зеа агеа. Клуег Ноу, 
уагуше ш уоате ап4 сВеписа| сотрозШоп, аНес з1етИсап у Фе уапаб у оЁ Бу4дгорВуз1са1 апа Мосеосвеписа1 
рагатеег$ оЁ Фе ргосеззез оссиггие ш Фе соаз{а1 епугойптеп. ТЬ1$ ехр|аштз Фе пееа Гог %Зайзиса| ргосеззше оф бе 
айа Кот 1опе-{егт обзегуанопз$ оуег Ше пуег Но\ сВагас{е!15Яс$. 

Маетай5 апа Мешоа$. Те тафетайса| то4е| ог Мосеосвеписа! субез 15 Базе оп а зу$фет оЁ поп-$Зайопагу 
едчайоп$ ог фе сопуесНоп-Ч1а$10оп-теасНоп оЁ рагабойс буре м попПпеаг РаисНоп$ оЁ зоигсез ап4 1о\уег-ог4дег 
депуайуез, ю \сЬ е соггезропд4те ша ап Боипдагу сопд 11015 аге аа4ед. шп Фе сопгзе о? ЗайзИса| апа1у$1$ оЁ 
Фе зетез ог 1опе-шеги обзегуайопз$ оуег пуег Йо\з, Фе уашез оЁ Фе ЮПоулие ш саг \еге Ючп4: $Ке\упез5 
сое слеп, 4еогее оЁ Кибо$1$, уапапсе ап ап аг4 демайноп, сое Исеп" оЁ уапайоп, амюсогге]айоп соеЁЯслеть, 
Меитапп гано, ап Апдегзоп сгцепоп. 

Кезий5. Тбе уайзИса| апа[у$15 оЁ Фе зетез оЁ 1оп2-{егт обзегуаноп$ оуег фе пу4госвепса| ша1саогз ог Фе Ооп 
пуег $102е55 Ваегорепейу оЁ Фе Не!4 д=а. ТЬ1$ 15 дие ю Фе зюсвазисйу оГ пиблет шри апа Фе уоште о? 
Кезбуужег Но\/ ю Ше зеа аз а гези оЁ пабига| ап ап@горогетс Ёасюгз. Е1е!4 даа зВоч]4 Бе соттёаеа \м зеазопа] 
свапоез ш Фе адиайс епутоптей фетрегате. ТЬ1$ рарег ргезеп5 Фе гези65 оРа сотриаНопа] ехрегитеп: © то4е| 
Фе Дупапл1с$ оР рвуюр!апКюп рорчайоп$ ш заттег зеазоп, \УПеп 1етрегаигез аге Рауога Ш е Юг фе гергодисйоп 
ап ото. ТБе ргорозе4 тафетайса| ое] сопз14Аегз фе зранаПу швотозепеочц$ 415б1БиНоп ап нап$огтаноп 
ог Ююгил$ о рвозрВоги$, ийгосеп, ап4 $Исоп, аз \еЙ аз свапоез ш заПипу, {етрегаиге, ап тойоп ог Ше адчаайс 
епугоптейе. 

О5си5яопт апа Сопсизоп5. Тве тизреслез та етайса| то4е] ог Ве дупапс$ ог рвуюр!апКюп роршаНопз 15$ 
соп514еге4 У ассоипё г Фе НапзЮгтайоп оЁР Юпиз$ оЁР рбозрвога$, пигозеп, ап $Исоп ш Фе соаза 
епугоптеп 5. ТБе апа[уз15$ оР даа Нот Не обзегуанопз, ог \мсЬ {$ та]ог Зайзиса! рагатаегз аге са]саеч, 15 
сагией ое. Аз а гези\, { 1$ сопса4е4 а даа оЁ Фе 1оп2егил обзегуайопз$ аге з1отиЙсайЙу уапа е. ТЬ1$ 1$ 4че © 
{\о геазоп$. Кап4от паге оЁ Ве 1приё оЁ пиблет5 ап4 Фе уо[ите оЁ пуег Но\м аз а гезий ог ап горозешс Гасогз 1$ 
Фе Нг5{ геазоп. ТВе зесопа геазоп шс[а4ез Ве аЦегпайоп оЁте[ануеу 61>1-уажег апа 1о\-\ужег рето4$ Гог НезВ Йом 
оуег Ше 1а3{ 12-15 уеагз. Тве пудгооз1са[ тесипе 15 свапоше таш]у дие ю Фе гедисйоп оЁ Фе ауегазе аппиа| 
Кезбуажег Но\у оЁ Фе Поп ап4 рагИу о{ Фе КорБап. Т№$ иепд 15 ПКеу ю шсгеазе дие ю сПтае сВапзез, аз \еП аз 
У Рег геошаноп ог Фе Поп пуег Но\у айег Ше Васаеузку пудго@есилсе шуаПаноп зай-ир. Митепса]| 
ехрегитеп!{ Базе оп фе Не]4 аа сопйгте4 Фе ргефсиуе уаПАКу оЁ Фе деуе]оре4 то4е[5 ап4 ргоэгапаз. ТВеу сап 
Бе изе4 №ю ргефсе свапое ш Фе сотшроз!юоп ап абипдапсе (сопсепайоп$) ш Фе Атоу зеа соге рапКюшс 
рорч1аНопз$, ”мсВ дейпе, оп Ще опе Вапа, Ююо4 гезотгсез, апа, оп Ше офег Вапа, Фе аацайс епутоптепе 1 {егил$ оЁ 
Фе опоошс зеа зайтитайоп. 


Кеумога$: Блозеосвеписа| сусез, руюрапКюп роршайНоп, Мосеплс заб$апсе, свеписа1-Б1о]ое1са| зопгсе, 


сопуесйоп-ЧИРа1оп-геасНоп еацайоп, Не]4 даа. 


Еог сйаноп: А. Г. ЗаКыпоту, У. У. Веота, А. У. МПайпа, е{ а1. Мо4дейи» Блосеосвеписа| ргосеззез ш Фе Атоу Зеа 
изшо ЗайзисаПу ргосеззе4 аа оп пуег Йо\. Адуапсеа Епошеегио КезеагсВ, 2020, уо]. 20, по. 4, рр. 437— 


445. ВИрз://401.0г2/10.23947/2687-1653-2020-20-4-437-445 


Сухинов А. И. и др. Моделирование биогеохимических процессов в Азовском море с использованием статистически обработанных данных 





Еипатэ трюгтапоп: Ще гезеагсй 15 4опе у Фе Йпапсла1 заррой Нот КЕЕГ (ргоуесй по. 20-01-00421). 


Введение. Азовское море — крупная прибрежная система эстуарного типа. Это самое мелкое море в 
мире. Оно почти равномерно прогревается летом (с перепадами температуры в среднем не более 4°С). При 
этом характеризуется большим перепадом солености — от 0 % до 12-15 %, т.к. стоки рек обеспечивают 
приток пресной воды, соизмеримый с общим объемом вод моря, а из Черного моря в районе Керченского 
пролива поступает соленая черноморская вода. Речной сток существенно влияет на биохимический состав 
водоема [1]. Представляется актуальным математическое моделирование биогеохимических процессов, 
позволяющее проводить диагностические и прогностические расчеты динамики экосистемы моря. Изменчивый 
по объему и по гидрохимическому составу речной сток заметно влияет на параметры гидрофизических и 
биологических процессов, происходящих в прибрежной системе. Поэтому целесообразно провести 
статистический анализ данных многолетних наблюдений, в частности гидрохимических показателей стока рек, 
впадающих в Азовское море, и на основе статистически обработанных входных данных выполнить 
прогнозирование биогеохимических процессов. 

В области исследования гидродинамики и прогнозирования морских систем следует отметить работы 
Марчука Г. И. [2], Матишова Г. Г., Сухинова А. И. [3], Бердникова С. В.., Тютюнова Ю. В. [4], 
Якушева Е. В. [5], Ильичева В. Г. и др. В статье представлены результаты комплексирования математической 
модели биогеохимических циклов с моделью гидродинамики Азовского моря [6-8]. Это позволяет повысить 
точность моделирования и учесть такие факторы, как гидродинамические процессы в прибрежных системах, 
неоднородное распределение температур, соленостей и биогенных веществ, влияющих на развитие 
фитопланктонных популяций, переход биогенов из одной формы в другую [9]. Следует заметить, что 
численные модели пространственно-трехмерных гидрофизических процессов в прибрежных системах являются 
предметом отдельного исследования авторского коллектива. Они дают возможность учесть динамически 
изменяющуюся геометрию дна и береговой линии, ветровое напряжение на свободной поверхности и ее 
возвышение, трение о дно, силу Кориолиса, турбулентный обмен, испарение, речные стоки, отклонение 
значений давления в водной среде от гидростатического приближения и др. В представленной работе входные 
данные (распределение трехмерного вектора скоростей, а также соленостей и температур) являются 
результатами численных расчетов на основе гидрофизической модели [10]. 


Материалы и методы. Для описания модели формулируется начально-краевая задача для системы 


уравнений параболического типа с младшими производными и нелинейными функциями правых частей: 
ое о 94 1, 
ее. т, р = 7(Когаа4;) + Кч,, (1) 
где 4; — концентрация [-й компоненты [мг/л]; (ЕМ, М = {ЁРу, Р›, Ез, РО4, РОР, РОР, МОз, №О>, МНа, 51; 


2 
{и,у, и} — компоненты вектора скорости водного потока [м/с]; А — коэффициент турбулентного обмена [м’/с]; 


Ка; — функция — источник биогенных веществ [мг/(л°с)]. 
В уравнении (1) индекс 1 указывает на вид субстанции (таблица 1). 
Таблица 1 
Биогенные вещества в модели динамики фитопланктона 


РО. 


РО 
ОР 
0: 





Химико-биологические источники описаны следующими уравнениями (1 Е 11, ул 3} ‚где 1 — это СЯТ, 


2— АК-А,3 — 5с,а СИТ, АЕ-А,5с — символические обозначения видов планктона): 


Кв; = Ск, (1 — Кев)ЧЕ; — Кер@ЧЕ; — КвЕЧв,» [=1Ъ3, 


3 
Кроь = 2.5КЕРЧЕ — К рр4рор — К рм@ров› 
1= 
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3 
Кроь = 2.5 КЕЧЕ + К ррЧ4рор — К риЧроь› 
3 
Кро, = 2.52 [К — | Ч + К руЧроь + К риЧрор › 


— 1 3 


Кин, = ие. а 1 з но, Чно ны, Че +2. 9х К +Кир | ЧЕ — К Ч кн, › 


1 . 


1 


= 


ое: 2. 
Е О Е 
>. а ыы | | Уи (4, , Ч хо, , Ч хн. ) Чо, ы Ч хо 


2 (4, › Чл, ›Ч и} 4 
т — у С те 1 | ПА = 
г. _ | } Ру "а ) Чо, ) Чл. ) Ч хо, и Ч лс 
К, — 5» (К, — 1) Ч», и 5% крЧь - 


Здесь К,, — удельная скорость дыхания фитопланктона; К,, — удельная скорость отмирания 


4, + К Чи, — К ›34мо, > 





Ч, + К›з4», › 


фитопланктона; К,, — удельная скорость экскреции фитопланктона; К,„, — удельная скорость автолиза РОР; 


К», — коэффициент фосфатофикации РОР; К», — коэффициент фосфатофикации ДОР; К, — удельная 


скорость окисления аммония до нитритов в процессе нитрификации; К›, — удельная скорость окисления 
нитритов до нитратов в процессе нитрификации; 5,, 5,, 5; —_ нормировочные коэффициенты между 


содержанием №, Р, эт в органическом веществе [11-12]. 
Скорость роста фитопланктона определяется выражениями 


С, № а Л: (Т) № (5) тт 41 (4„», | р (4. Чо, › Ч хн. )} 2 
С № Кн Л: (т). (5) тт т (4„„ а А - (4) > 


где Кл — максимальная удельная скорость роста фитопланктона. 


Зависимости температуры и солености: 
Г — Ре 5 — О ори 


ор ор 


где Г„, 5,„ — температура и соленость, оптимальные для данного вида фитопланктона; «>0, р>0 — 


ор: 
коэффициенты ширины интервала толерантности фитопланктона к температуре и солености соответственно. 








Рис. 1. Модельная схема биогеохимической трансформации форм фосфора, азота и кремния 


Ниже представлены функции, описывающие содержание биогенов. 


Ч ро. 





‚где К„, — константа полунасыщения фосфатами. 


Для фосфора У, (4„,)= Е 


Ч 
Ч; 5 
Для азота Х, (4 ‚Ч, Ч, }= Л (Ч, ›Ч,ь Чи, ) + ГО (а. ) › 


(1) — (Чо, че Ч хо, ) ехр(-К „,; Чун. ) (2) Е 
но 


Для кремния /, (аз) = ‚ где К, — константа полунасыщения кремнием. 
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где К мо — Константа полунасыщения нитратами; К мн, —` Константа полунасыщения аммонием; Ре 


коэффициент ингибирования аммония. 

Предположим, что коэффициенты, входящие в выражения для функций источников, положительны и 
не зависят от времени [. 

Для системы (1) ставится начально-краевая задача в цилиндрической области С. Пусть граница У 


цилиндрической области С’ является кусочно-гладкой поверхностью и Х=»Х,„ \/», о , где Х„ — поверхность 


дна водоема, », — невозмущенная поверхность водной среды, с — боковая (цилиндрическая) поверхность. 


Пусть и, — нормальная по отношению к » составляющая вектора скорости водного потока, п — 
вектор внешней нормали к »%. Допустим для концентраций 4; на боковой границе: 








4, =0, на с, если и, <0,1ЕМ; (2) 
О0. 
С 0, на с, если и, 20,1еМ; (3) 
[6] 
О4. 
— =0, на », — поверхность водоема, 1{ЕМ; (4) 
р 
да. да. 
= 81:9; 1Е Е» Яо, Е}, — 8> 4; › 
02 02 
ге {РО4, РОР, РОР, МОз, №О›, МНа4, 5} на дне У». (5) 


Здесь &,,, &., — неотрицательные постоянные; &,,;, 1 {Р\, Р›, Рз}; учитывают опускание водорослей на дно и 


их затопление; =,,, {= {РО4, РОР, РОР, МОз, МО», МНа4, 51} учитывают поглощение питательных веществ 


донными отложениями. 
Для системы уравнений необходимо в любой момент времени задавать вектор скоростей водного 
потока, поле солености и температуры, а также начальные значения функций 4, : 


Ч, (х, у,2,0) = Ч; (х, у,2), (х, у,2) ЕС, 1=0,1ЕМ, (6) 


У (х, у, 2,0) = \ (х,у,2), Т(х,у,2,0) = (х,у,2) Ь 5 (х,у,2,0) = 5) ) 
Статистическая обработка данных многолетних наблюдений о стоках рек в Азовское море. 
Значительный относительно объемов моря речной сток существенно влияет на биологические и 
гидрофизические процессы, происходящие в Азовском море [13]. Со стоками рек в водоем поступает большое 
количество биогенных веществ, в том числе азота, фосфора и кремния — основных питательных элементов для 
фитопланктона. В ХХ веке основная часть притока вод в Азовское море приходится на сток Дона — 


63 % (рис. 2—4)". 
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Рис. 2. Ряды многолетних наблюдений за стоком реки Дон (1993—2012 гг.): концентрация азота (№-МН4) 
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Рис. 3. Ряды многолетних наблюдений за стоком реки Дон (1993—2012 гг.): концентрация фосфора (РО4) 


' Экологический атлас. Черное и Азовское моря / НК «Роснефть» ; Арктический научный центр ; НИР. М., 2019. 464 с. 
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Рис. 4. Ряды многолетних наблюдений за стоком реки Дон (1993—2012 гг.): концентрация кремния (5104) 


Таблица 2 
Результаты расчета статистических параметров данных натурных наблюдений 





Исследование рядов многолетних наблюдений за стоком Дона позволяет сделать ряд выводов. 

—[ Рассматриваемые концентрации питательных веществ обладают как положительной, так и 
отрицательной асимметрией. 

— Случайные величины для азота и кремния смещены относительно центра распределения, о чем 
свидетельствует высокое значение коэффициента асимметрии. 

— Большие значения дисперсий и среднеквадратических отклонений получены для всех биогенов. 

— Коэффициенты автокорреляции малы, следовательно, для рядов натурных данных характерна 
сильная нелинейная тенденция. 

— Вариация во всех рядах больше 20 %, следовательно, ряды сильно вариативны. 

— Для азота очевидно наличие автокорреляционной связи по критерию Андерсона при количестве 
значений в выборке 20, так как коэффициент автокорреляции превышает 0,299 при уровне значимости 5%. 

— Для азота очевидно наличие автокорреляции остатков по критерию Неймана при уровне значимости 
5% для 20 наблюдений, т. к. отношение Неймана меньше 1,2. 

— Фосфор и кремний не демонстрируют автокорреляционные связи, гипотезу об автокорреляции 
остатков отвергаем. 

В результате статистического анализа натурных данных [14] можно сделать вывод об их сильной 
изменчивости. Это обусловлено стохастичностью поступления биогенных веществ со стоком Дона и 
существенно меняющимся объемом стока под воздействием природных и антропогенных факторов. Для 
использования натурных данных в модели (1)-(6) целесообразно учитывать сезонные изменения. Далее при 
моделировании будем рассматривать летний период. 

Результаты численных экспериментов. Проведено численное моделирование решения задачи 
динамики фитопланктонных популяций в летнее время с учетом трансформации форм фосфора, азота и 
кремния на примере Азовского моря. Моделируемая область соответствует физическим размерам Азовского 
моря: длина — 355 км, ширина — 233 км, шаг по пространству в горизонтальных направлениях — 1000 м. На 
рис. 5 приведен спутниковый снимок Азовского моря, подтверждающий соответствие результатов работы 
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натурным данным. На снимке отчетливо видны распределения зеленых и сине-зеленых водорослей в области 
Таганрогского залива и диатомовых водорослей — в центральной части моря. 





Рис. 5. Спутниковый снимок Азовского моря, сделанный с помощью спектрорадиометра с умеренным разрешением 
(МОГ) спутником МАЗА Адра 31 июля 2004 года 


В результате вычислительного эксперимента получены сеточные распределения концентраций 
основных популяций фитопланктона и питательных веществ в Азовском море (рис. 6). Период расчета — 30 
суток. Этого достаточно для установления стационарных режимов в задачах динамики фитопланктона. 
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Рис. 6. Распределения концентраций: зеленой водоросли Сй[огеЙа уш[ватх (а); сине-зеленой водоросли Арйаттотепоп 
Ло5-адиае (6); диатомовой водоросли эсе/еюпета созшиит (в); фосфатов (г); нитратов (0); 
растворенного неорганического кремния (е) 


Рисунки отражают динамику фитопланктонных популяций, циклы фосфора, азота и кремния. В 
процессе экскреции и отмирания фитопланктон выделяет фосфор в растворенной и взвешенной органических 
формах, затем в процессе фосфатофикации они переходят в неорганическую форму — фосфаты, которые 
потребляет фитопланктон. Также описан цикл азота: в процессе жизнедеятельности фитопланктон выделяет 
азот в органическом виде, который разлагается до аммиака. Аммиак в процессе нитрификации окисляется до 
нитритов, а затем до нитратов. Стоит отметить, что фитопланктон потребляет все три формы азота. Отмечено 
потребление и выделение кремния диатомовыми водорослями. Сравнение с результатами моделирования для 
многоводных периодов показывает, что в последние маловодные годы ареалы обитания зеленых и сине- 
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зеленых водорослей в области Таганрогского залива существенно (на многие километры) сместились в 
восточном направлении, ближе к Дону — источнику пресной воды. 

Обсуждение и заключения. В работе приведена многовидовая математическая модель динамики 
фитопланктонных популяций, трансформаций форм биогенных веществ — соединений фосфора, азота, 
кремния. В модели учтено: 

— влияние солености и температуры на развитие трех основных видов фитопланктона (зеленых, сине- 
зеленых и диатомовых водорослей); 

— поглощение фитопланктоном фосфатов и форм азота, 

— переход форм фосфора и азота из одной в другую, 

— поглощение кремния диатомовыми водорослями, 

— адвективное и микротурбулентное движение водной среды, 

— стоки и источники на границе. 

Для анализа натурных данных разработана и адаптирована методика статистического анализа 
многолетних рядов наблюдений концентраций биогенных веществ (фосфора, азота, их соединений и др.), 
которые поступают в море со стоком Дона. Статистический анализ данных многолетних наблюдений, в 
частности гидрохимических показателей стока Дона, стал базой для прогнозирования биогеохимических 
процессов с учетом движения водной среды, распределения температур и солености. Результаты численных 
экспериментов согласуются с данными космического зондирования Азовского моря, что является 
подтверждением предсказательной ценности применяемых моделей и методов их численной реализации. 
Сравнение распределений популяций зеленых и сине-зеленых водорослей в Таганрогском заливе для 
многоводных и маловодных периодов показывает, что ареалы их обитания в последние маловодные годы 
существенно (на многие километры) сместились в восточном направлении, ближе к Дону — источнику пресной 
ВОДЫ. 
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